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Introduction 

Satellite systems to date have been mainly scientific in nature. Only a few 
systems have been of direct use to the public such as for telephone or television 
transmission. Space enterprises have remained a mystery to the general public and 
beyond the reach of the small business community. The result is a less than 
supportive public when it comes to space activities. 

The purpose of the ISAT-1 program is to develop a small and relatively 
inexpensive satellite that will serve the State of Iowa, primarily for educational 
purposes. It will provide products, services, and activities that will be educational, 
practical, and useful for a large number for people. The emphasis, is on public 
awareness, "space literacy", and routine practical applications rather than high 
technology. 

The initial conceptual design phase was complete when the current team 
took over the project. Some areas of the conceptual design were taken a little 
farther, but for the most part this team started at the detailed design stage. This 
team was split into eight groups: 

• Integration 

• Structures 

• Payload 

• Power 

• Thermal 

• Attitude Determination and Control 

• Ground System and Operations 

• Communications 

Each groups work is expanded upon in the following report. Each chapter 
relates the group work to date and any considerations that the group felt were 
necessary for future work. 




Chapter 1 
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Chapter 1 Integration 
1.1 Introduction 

The integration team's responsibility was to combine all incoming designs 
from other teams and combine them into one feasible design. Determine packing 
considerations in the satellite and construct at least one full scale model of the 
satellite, for packing considerations. Determine a mass and power budget for the 
satellite and select a launch vehicle for the satellite. 


1.2 Launch Vehicle 

The Delta II was selected as the launch vehicle. Several options were 
investigated including the Space Shuttle, Ariane and Atlas. The Shuttle was not 
chosen because it doesn't reach a high enough altitude. The Atlas was ruled out 
because it has no designated secondary payload accommodations and the primary 
payload option was prohibitively expensive. The Ariane was not used because the 
launch apparatus and configuration were not compatible with the proposed ISAT- 
1 design. The only drawback to using the Delta II was that the largest secondary 
payload envelope offered was not large enough to accommodate a satellite the size 
of the initially proposed ISAT-1. The size of IS AT- 1 was therefore reduced to fit 
into the available envelope. The current launch date is set for December 1995. 
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1.3 Mass Budget The current mass budget is reflected in Table 1.1. 


# 

Component 


1 

Structure 

3.6 

1 

G. G. Boom 

7 

2 

Communications 

3 

2 

Flight Computer 

3 

2 

Large Experiment Module 

3 

2 

Medium Experiment Module 

1.1 

6 

Small Experiment Module 

0.9 

1 

Pitch Torqrod 

0.1 

1 

Roll Torqrod 

0.1 

1 

Yaw Torqrod 

0.45 

2 

G. G. Boom Tanks 

0.1 

1 

Earth Radio Frequency Experiment 

5.4 

2 

CCD Camera 

0.39 

1 

Telescope 

3 

6 

Battery 

0.6 

6 

Battery Mount 

0 

2 

Magnetometer 

0.22 


Table 1.1 Mass Budget 










































Chapter 1 Integration 


4 


1.4 Power Budget The current power budget is reflected in Table 1 .2 


# 

Component 


1 

Structure 

0 

1 

G. G. Boom 

0 

2 

Communications 

7 

2 

Flight Computer 

7 

2 

Large Experiment Module 

4 

2 

Medium Experiment Module 

1.5 

6 

Small Experiment Module 

1.1 

1 

Pitch Torqrod 

0.25 

1 

Roll Torqrod 

0.25 

1 

Yaw Torqrod 

2 

2 

G. G. Boom Tanks 

0 

1 

Earth Radio Frequency Experiment 

7 

2 

CCD Camera 

5 

1 

Telescope 

0 

6 

Battery 

0 

6 

Battery Mount 

0 

2 

Magnetometer 

0.04 


Table 1.2 Power Budget 
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1.5 Packaging 

For the integration and packaging of all the components of the bus 
subsystems and payloads, the major factor in choosing the components was their 
size. Also playing a major part in the choice of components was cost. The 
following is a summary of the bus subsystems to be included in the satellite: 

a. (2) magnetometers, by ITHACO, chosen because of their size. 

(Figure A. 1.1) 

b. (3) Torqrods, by ITHACO, chosen because this was the only place they 
could be found. There are two 1 a*m 2 torqrods, and one 10 a*m 2 
torqrod. (Figures A. 1 .2 and A. 1 .3) 

c. (6) batteries, by Eagle Picher, which are nickel hydride inside, chosen 
because of durability. (Figure A. 1 .4) 

d. (2) flight computers, which will be made by Iowa State University. 

(Figure A. 1.5) 

e. (2) communication units, made by Cynetics. They were chosen because of 
their relatively small cost. (Figure A. 1 .6) 

f. (1) telescope made by Iowa Scientific Optical, chosen because they could 
build it to specifications. The CCD camera to go with it is built by Canon. 
(Figure A. 1.7) 


The payloads, up to this point are specified below: 

a. Earth Radio Frequency Experiment (ERFE). (Figure A. 1 .8) 

b. Two large experiment modules, one which will hold the seed experiment. 
The other is uncommitted. (Figure A. 1 .9) 

c. Two medium experiment modules, both which are uncommitted. 

(Figure A. 1.10) 

d. Six small experiment modules, all of which are uncommitted except one, 
which will hold the electronics needed for the robot arm. (Figure A. 1.11) 
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The experiment modules will be constructed of 6061-T6 aluminum. The 
satellite structure will be graphite epoxy. The solar cells are gallium arsenide 
germanium. The packaging arrangement for the satellite is seen in Figure A. 1.12. 
There are three shelves to which the various components can be attached. The 
satellite is packed in this configuration because it's the only way everything would 
fit. 

1.6 ISAT-1 Model 

Part of the integration team's job was to build a full scale mock up of the 
proposed satellite. This project was started just after mid-term spring semester 
1993. As of this date 95% of the payload and 90% of the bus components are 
complete, just needing to be painted and mounted on the appropriate shelves 

During the planning stages of the model construction, size of all internal 
components needed to be considered along with the size of the launch vehicle s 
cargo area. The cargo area defined the maximum dimensions of the satellite and 
the size of the internal components dictated where they would be placed, besides 
restrictions due to the component itself, such as the telescope needing to be 
earthward looking, or the magnetometers needing to be mutually perpendicular. 
The packing arrangement finally arrived as in Figures A. 1.12 and A. 1.13. This 
helped during the finalization of the packing arrangement. 

The model is being constructed out of extruded polystyrene and modelers 
plywood. The polystyrene is being used for the solid models of the bus and 
payload components, along with portions of the satellite outer surface. Modeler's 
plywood is being used where load considerations would make the use of 
polystyrene questionable; such as the shelves, mountings for the bus and payload 
components, and outer skins. Due to the difficulty in manufacturing cylindrical 
parts out of polystyrene, the bodies of the batteries were cut from plumbing pipe of 
an outer diameter matching the actual batteries (2 1/2 inches). 
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Model construction is expected to be complete by the end of June 1993. 
This will include two full scale mock-ups, complete with internal bus and payload, 
and full scale models of all internal components for individual display. 


1.7 Future work 

Further work for the satellite configuration will happen after every design 
iteration. The integration team will continue to repackage until the final design. 
The integration team will also have to repackage probably to account for wiring. 


1.8 References 

Delta II Launch Planners Manual, Secondary Payload. 


Chapter 2 
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Chapter 2 Structures 
2.1 Introduction 

The structure of the spacecraft provides the necessary physical support for 
all subsystems. The structure must not only contain the payload and bus systems 
while on-orbit, but must withstand launch conditions and ensure that all of the 
payload and bus systems will arrive on-orbit in an operational condition. 

Previous work on the ISAT structure included the determination of the 
shape of the exterior. This was selected as hexagonal based on solar cell power 
calculations. In addition, the maximum dimensions of the exterior were 
determined based on the launch vehicle envelope chosen (Delta II SPE1). A 
previous analysis of a metal truss structure had been conducted by the Aerospace 
Engineering 464 class of Spring 1992. In addition, the internal configuration was 
provided by the Integration Group. 


2.2 Structure Performance Criteria 

Several criteria for the performance of the ISAT structure were developed. 
These included the requirement of supporting the internal payload arrangement as 
specified by the Integration Group, supporting the McDonnell-Douglas specified 
SPE1 launch adapter, supporting attitude control mechanisms, and using as little 
internal volume as possible. 

The internal configuration of ISAT-1 is shown in Figure 2.1. This 
arrangement shows the payload and bus systems arranged on three internal shelves. 
The payload consists of a series of differing size boxes, for which the mass has 
been estimated. 

The location and arrangement of attitude control mechanisms are also 
shown in Figure 2.1. Most notably, the structure must be capable of transmitting 
control torques from the three torque-rods and from the gravity-gradient boom 
attached at shelf 2. 
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Figure 2.1 IS AT Internal Configuration 


The McDonnell-Douglas launch adapter for SPE1 payloads is shown in 
Figure 2.2. This adapter attaches to the payload with six (6) bolts as shown. This 
attachment configuration must be supported by one side panel. Figure 2.3 shows 
the geometry of the launch adapter attachment points with respect to the ISAT 
panel design. 
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2.3 Structure Concepts 

Five concepts for the internal structure of ISAT were considered. They 
included Central Support, Modified Central Support, Full Monocoque, Semi- 
Monocoque, and Truss. These five concepts are shown in Figure 2.4. These 
structural concepts were considered prior to determination of the internal payload 
arrangement. 



Figure 2.4 Internal Structure Concepts 


The Central Support concept would place a large structural member along 
the central, longitudinal axis of the satellite. Payloads would then be attached to 
this structural member, with the skin panels being supported by additional braces. 

The Modified Central Support concept would place the large structural 
member along one side of the hexagon structure to accommodate the launch 
adapter attachment for the Delta II SPE1. Again the payloads would be attached 
to the major structural member and the skin panels supported by additional braces. 
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The Full Monocoque concept would consist of thick skin panels that would 
be designed to support the payloads either by direct attachment or by attachment 
to cross shelves. The panels would have to be designed thick enough to transmit 
the launch adapter loadings as well. 

The Semi-Monocoque concept would consist of skin panels attached to 
longitudinal stringers. The design would reduce the panel thickness while 
transferring more structural material into the comers of the hexagon. Again, one 
panel would have to be designed to accommodate launch loads. 

The Truss concept would consist of thin skin panels attached to a 
framework of structural elements. This framework would support the internal 
arrangement of payloads. 


2.4 Structure Selection Criteria 

One of the major considerations in the final selection of a structural 
concept was the amount of internal volume that would be required. It was felt that 
the semi-monocoque structure would best fulfill these criteria, as it would make 
the best utilization of previous unused space in the comers of the hexagon shape. 

Another consideration was the ability of the structure to be adapted to the 
unique launch interface. The placement of the launch interface on the side of the 
satellite, require that non-uniform loading be a consideration. It was felt that the 
semi-monocoque structure would again adequately meet this consideration. The 
semi-monocoque structure can be adapted to the launch adapter by varying the 
stringer and panel thickness on the side of the launch adapter. 

Fabrication and assembly of both the structural components and the 
integration of the payload were also considered. The semi-monocoque concept 
provides for greater flexibility in fabrication and assembly sequences than the other 
concepts. It should be possible to design the structure such that all segments of 
the payload are accessible even with the launch adapter in place. 
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2.5 Preliminary Structure Sizing Analysis 

A preliminary analysis of the structure was conducted for sizing purposes. 
For this case, each payload was considered to be of uniform mass, and its mass 
would act at its center of mass on each shelf. These loads were then transferred to 
the six stringers for each shelf. The details of the this analysis are presented below 
in Table 2. 1 . The mass budget may be seen in Table A.2. 1 . 


Bottom 

0.0 

Shelf 1 

92.4 Newtons 

Shelf 2 

138.5 Newtons 

Shelf 3 

156.0 Newtons 

Top 

to be determined 


Table 2.1 Static Shelf Load Estimates 

These loads were then used to conduct a stress evaluation of the stringers 
and the shear stress in the panels. This analysis was based on the asymmetric beam 
bending equation. 

The program calculated the axial stress in the stringers and the shear stress 
in the panels due to the maximum shelf loading case. The code may be found in 
Appendix B.2. Results indicated that for the static case, a panel thickness of 0.040 
inches and a stringer with flanges 15% of the panel width and a thickness of 0.080 
inches should be adequate. 


2.6 Preliminary Finite Element Analysis 

The results from the preliminary sizing analysis were used to construct a 
finite-element model. This model used the panel and stringer thicknesses 
mentioned above. The element plot for the internal configuration is shown in 
Figure 2.5. The element plot for the external panel configuration is shown in 
Figure 2.6, and the figure for the external stringer configuration is shown in Figure 
2.7. An additional element plot may be found in Figure A.2.1. The model was 
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fixed in rotation and translation at the four launch adapter attach points. The 
model was loaded with the estimated shelf loads. 



Figure 2.5 ANSYS Internal Element Plot 



Figure 2.6 ANSYS External Element Plot - Panels 
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Figure 2.7 ANSYS External Element Plot -- Stringers 


2.7 Structural Elements 

2.7.1 Gravity-Gradient Boom Design 

The structures group also performed the design of the gravity gradient 
boom. For additional information on the Gravity-Gradient Boom see Appendix 

C.2. 

2.7.1. 1 Boom Section Selection 

The closed-section design was selected for further development. This 
section was chosen for its torsional properties and its uniform bending stiffness. 
The simplicity of the deployment mechanism was also a consideration. The current 
assumption is that the magnetic portion of the satellite attitude control system will 
be able to invert the satellite with the boom deployed, so the boom will only be 
deployed once. 
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Figure 2.8 and Figure 2.9 show the details of the joint design. Figure 2.8 is 
a cut-away of the outer wall of a section. Figure 2.9 is a view of an inner section, 
with the outer section not shown. This design provides for a forward stop collar 
located on the inner wall of the outer boom section, and a series of forward 
segmented stops located on the outer wall of the inner boom section. A 
segmented stop collar is also located on the inner wall of the outer section. A rear 
stop collar is also located on the inner section. When pressurized, the boom 
extends until the forward stop collar on the outer section comes into contact with 
the segmented collar on the inner section. The segments on the inner section slide 
through the gaps in the collar on the outer section, and are then rotated to lock the 
boom in its deployed state. The joint design also provides for the placement of 
three gaskets: one each at the base of the inner section, the top and base of the 
rear stop collar to provide a gas seal for the pressurized deployment of the boom. 



Figure 2.8 Joint Detail A 




Chapter 2 Structures 


18 



Figure 2.9 Joint Detail B 


2.7. 1.2 Boom Sizing 

The operational task of a gravity-gradient boom is to modify the inertia 
properties of a spacecraft to increase its stability. The degree of stability increase 
desired is the driving factor in sizing the boom's length and tip mass to a particular 
satellite. Iowa State University is currently undertaking an effort to design, 
construct and launch a small satellite, ISAT-1 (Iowa Satellite One). The 
preliminary design of this satellite was used to size the boom length and mass for 
this project. 

Current efforts of the Iowa Satellite Project provided much of the needed 
information about the satellite attitude dynamics. The attitude determination and 
control group provided the mass moment inertias of the spacecraft body. These 
inertia calculations assumed uniform mass distribution within the satellite. A 
design sizing code with boom length and tip mass as control variables was written 
to determine an optimum length and tip mass for the satellite. The code allows the 
user to select the desired stability in terms of the stability criterion 0 X . The user 
also selects the desired range of boom lengths and tip masses. The code then 
iterates through these two variables, calculating the new mass moment of inertia 
I x , and uses the new inertia to calculate the value of 0 X . The calculated value of 
0 X is subtracted from the target value, and the absolute value of the difference is 
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written to a data file along with the boom length and mass. This information can 
then be plotted on a contour plot. Figure 2.10 shows the contour plot for a 
stability criterion (0 X ) of 0.98. Several candidate designs meet the desired 
stability. The final boom size was selected to be 3 meters in length with a 6 
kilogram tip mass. This represents almost a fifteen fold increase in I x . 



Stability considerations determined the boom length and tip mass. The 
diameter and thickness of each section would be determined according to the loads 
acting upon each section Loads acting on the gravity-gradient boom were 
modeled as the Earth's gravitational force acting on the tip mass. This in turn 
caused forces and moments to be exerted on the boom. A section sizing code was 
written to size the boom sections based on the known loads, the material 
properties of graphite-epoxy composites, and a selected factor of safety. It was 
not possible to accurately determine dynamic loads on the gravity-gradient boom 
due to insufficient source data from the ISAT project, so the factor of safety used 
in the sizing analysis was increased. 

The code was set up to allow the user to choose the desired tip mass and 
boom length for section sizing analysis. The forces and moment on the boom were 
then determined based on the user's input of the tip mass and the maximum 
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deflection of the satellite from the vertical axis. Once the user had completed the 
inputs of number of sections, minimum section thickness, gap between sections, 
and initial base section outer diameter, the code would size the boom sections. 
This entailed determining the necessary outer diameters, section thicknesses, and 
inner diameters corresponding to each boom section. If the section thickness 
needed was found by the code to be smaller than the minimum input section 
thickness, the code would substitute in the minimum section thickness and 
continue the analysis. This measure was included in order to allow for more 
reasonable composite layup thicknesses in the fabrication process. 

Results of section sizing yielded necessary section thicknesses for a number 
of base section outer diameters. A base section having an outer diameter of 
0.0381 meters was chosen based on the limited volume fraction of ISAT-1 that 
would be required. Table 2.2 shows the section sizes in terms of outer diameters, 
thicknesses, and corresponding inner diameters. It should be noted that in all cases 
the thicknesses used are the minimum section thickness, and the gap between 
sections is constant for ease of fabrication. 


Boom 

Section 

(#) 

Outer 

Diameter 

(m) 

■KIM 

Inner 

Diameter 

(m) 

1 

0.0381 

0.00102 

0.0371 

2 

0.0361 

0.00102 

0.0351 

3 

0.0340 

0.00102 

0.0330 

4 

0.0320 

0.00102 

0.0310 

5 

0.0300 

0.00102 

0.0290 

6 

0.0279 

0.00102 

0.0269 


Table 2.2 Boom Section Sizes 


2.7.1.3 Fabrication Methods 

Fabrication methods were developed to produce the necessary sections to 
conform to the design. This included the construction of suitable mandrels and 
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development of layup procedures. This also included developing a fabrication 
technique for the section joints. 

The mandrels used for curing of the boom sections consisted of two 
sections of steel pipe, each 0.61 meters in length. The two pipes had inside 
diameters of 0.0343 meters and 0.0394 meters. These pipes were split along their 
length at the diameter. These were complemented by three rubber hose sections of 
0.76 meters in length. The hoses had outside diameters 0.0262 meters, 0.0315 
meters, and 0.0361 meters. This provided the capability to fabricate boom 
sections of varying diameter and thickness. The hoses were fitted internally with 
steel rods along their entire length. Both ends of the hoses were sealed, one 
containing a valve, to allow pressurization during the curing process. 

Each section was fabricated from three plys of woven [0°,90°] graphite- 
epoxy composite. The composite material was cut to a length of 0.56 meters and a 
width equal to the circumference of the desired section. The three plys were then 
debulked with the edges staggered prior to being placed in the mandrel. The 
debulked composite was then placed in the mandrel and wrapped around the hose. 
Finally, the pipe section was clamped around it. Bleeder cloth was inserted 
between the composite material and the outer pipe as necessary to form each 
section. 


With the mandrel clamped together, the inner hose was pressurized to 
275,790.3 Pa (40 psig), and the mandrel was placed in an oven to cure for 3 1/2 
hours at 176° C. Upon completion of the cure, the mandrel was disassembled, and 
the completed boom section was removed. 

To avoid the added complexity and difficulties involved with co-curing the 
joints along with the respective boom sections, joints from separate cure cycles 
were epoxied in at a later time. A suitable gas seal was added to each section, and 
the boom was assembled. 
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2.7.1. 4 Testing 

Once the boom was complete, testing was conducted to ensure that the 
mechanical performance of the gravity-gradient boom design was satisfactory. 
The first test to be run was an analytical test involving ANSYS finite-element 
analysis. The boom was modeled on ANSYS using plate elements. A static 
analysis was performed. Figure 2.11 gives a representation of the deflection 
resulting in a maximum load condition applied to the boom and its tip mass. The 
boom was found to have a maximum deflection at the tip of 0.05 meters (1.969 
inches). 



Figure 2.11 ANSYS Model 


2.7.2 Shelving 

2.7.2. 1 Design Considerations 

Each shelf must be capable of supporting all loads and payloads. In 
addition the shelves must accommodate the payload boxes, including any unique 
interfaces (radiator, tape, electrical/power connections, etc.). Some payloads may 


Chapter 2 Structures 


23 


also require vibration isolation. The internal volume of the spacecraft and the 
payload requirements limit the allowable deflection for each shelf. Other concerns 
include thermal conductivity, vibration characteristics, cost, weight, ease of 
fabrication, and ease of access. 


1 . 1.12 Design 

Three general types of shelving are available: solid aluminum, or graphite 
epoxy with cut out sections, stringers for supporting payloads with skin as needed, 
thin honeycombed construction. 


1 . 1.12 Assumptions 

Preliminary modeling and analysis of the shelves assume that the payload 
containers carry no load. The payloads are modeled as point masses on 
appropriate elements of the shelf. 


2.7. 2.4 Results 

The preliminary shelf modeling was included in the ANSYS finite element 
model. Appendix B.2. Detailed analysis of individual shelf behavior has not been 
completed. The ISAT structural mass estimate may be seen in Table A.2.2. 


2.7.3 Launch Adapter 

The launch adapter ring provided by McDonnell-Douglas for the SPE1 
payload fairing on its Delta II rocket is shown in Figure 2.2. This adapter attached 
to the launch vehicle in the payload envelope shown in Figure 2.12. The ISAT 
structure must be designed to accommodate this launch adapter and fit within the 
designated payload envelope. 
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Figure 2.12 Delta II SPE1 Orientation 


2.7.3. 1 Payload Envelope Considerations 

The ISAT main structure has been sized to fit within the allowable payload 
envelope. This results in a regular hexagon with a side length of 0. 170 meters, and 
a height of 0.643 meters. This size is shown in Figure 2.1. The allowable and 
negotiable space must also accommodate the gravity-gradient boom that extends 
above the satellite 0.2134 meters without the tip mass. 


2.7.3.2 Launch Adapter Considerations 

The launch adapter design is fixed by McDonnell-Douglas Space Systems. 
The adapter consists of a ring that is clamped to the launch vehicle. Six 
attachment bolts are provided to attach the launch adapter to the payload. The 
ISAT structure must be designed to accommodate the given interface and be able 
to jettison the launch adapter when orbit is achieved. The size of the ISAT panel is 
fixed by the size of the launch envelope. The width of the launch adapter exceeds 
that of one side panel. Therefore, two of the attachment bolts will not directly 
contact any side panel. The geometry of the launch adapter and the ISAT side 
panel is shown in Figure 2.3. 
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Several designs have been considered to meet this challenge. One 
possibility is the reinforcement of the side panel facing the launch adapter by 
addition of stringer, increasing the size of existing stringers and/or increasing panel 
thickness. This design may make it possible to not use the two attachment bolts 
not directly in contact with the panel. Further study of this design should .require 
consultation with McDonnell-Douglas. 


2.8 Recommendations for Future Work 

The dynamics loads during the launch phase are of primary concern. In 
addition, on-orbit attitude dynamics such as tumbling, de-spin, and inversion need 
to be characterized and modeled. 

A detailed design of the launch adapter attachment to the satellite needs to 
be conducted to determine the transfer of loads from the attachment points to the 
remainder of the structure. McDonnell-Douglas needs to be consulted concerning 
the possibility of using only four of the six available launch adapter attachment 
bolts. 


A more detailed study of the trade-offs involved with panel thickness and 
stringer area needs to be conducted. 

The boom attachment mechanism and location of support hardware need to 
be designed. 

Additional research is necessary to understand the effects of the space 
environment on composite materials. 

Additional structural requirements, such as deflection limits, thermal 
considerations, payload interfaces, and payload geometry, need to be compiled 
from other groups. 

Further work is needed to understand the application of graphite-epoxy 
composites to gravity-gradient booms. The feasibility of the application has been 
demonstrated. However additional development work and testing are needed to 
provide the validation necessary to prepare hardware for use in a space application. 
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Chapter 3 Payload 

3.1 Earth Radio Frequency Experiment 

3.1.1 Introduction 

The Earth Radio Frequency (ERF) Experiment will measure the intensity 
and spectrum of terrestrial communication signal "leakage" through the Earth's 
ionosphere over the frequency range of 1.5 - 34.5 MHz. This experiment will 
collect data that, when analyzed on the ground, will help determine the spectrum, 
intensity, and temporal characteristics of this signal leakage. This data will be used 
to determine whether low frequency radio astronomy is feasible from Earth orbit. 

The ERF Experiment consists of a swept LO heterodyne receiver with 
analog-to-digital conversion, some basic signal processing, and data storage 
elements, Figure 3.1. 



voltages from bus 


Figure 3.1 ERFE Block Diagram 
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3.1.2 Mechanical Requirements 

The ERF receiver shall be designed and built to withstand the hostile 
environment or space and in such a manner that it will successfully accomplish its 
mission. The following paragraphs provide detailed mechanical design 
requirements to meet these goals. 


3. 1.2.1 Physical Dimensions 

The experiment shall be contained in an enclosure measuring 6" wide by 8" 
long by 4" high. 


3.1.2.2 Enclosure 

The enclosure containing the ERF experiment shall be made of 0.25" thick 
machined aluminum stock to provide adequate radiation shielding for the RF, 
analog, and digital circuitry inside. The enclosure shall be constructed in two 
sections. When the receiver is operational, the circuitry must be sealed inside the 
enclosure using non-outgassing materials. 


3. 1.2.3 Connectors 
3. 1.2.3. 1 Power 


The connector providing the supply voltages to the receiver from the 
power bus shall be of (TBD) type, with the following pin output: 


1) +5V 

2) +15V 

3) -15V 
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3.1.2.3.2 Antenna/RF Input 

The antenna/RF input connector shall be a 50 O, type (TBD) connector. 


3.1.2.3.3 Digital Data Input/Output 

The connector(s) for the receiver's digital command inputs and data 
outputs shall be the (TBD) pin D sub miniature type, with built-in filtering. 


3.1.3 Environmental Requirements 
3.1.3.1 Temperature Range 

The receiver must meet performance requirements over the following 
temperature ranges: 

-30° to +70° C operational 

-40° to +80° C non operational 

The receiver must be capable of functioning within design specifications 
after withstanding temperatures in the non operational range. 


3.1.3.2 Vibration 

(to be determined) 


3. 1.3.3 Radiation 


(to be determined) 
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3.1.4 Power Consumption 

The receiver shall have a maximum power consumption of 7W while 
operational. Its standby power consumption shall be (TBD) W. 


3.1.4.1 Power Conditioning 

The receiver should be capable of conditioning the power received from 
the spacecraft. The available input voltages shall be : 

+ 5 V ± (TBD) % digital circuitry voltage 
± 1 5 V ± (TBD)% analog/RF circuitry voltage 

3.1.4.2 Power Stability 

The power conditioning circuitry in the receiver shall provide a stable 
source of power to the instrument. The voltages shall not vary by more than 
(TBD)% over a period of 12 hours. 


3.1.5 Future Work 

Additional work is required in the following areas: 

1) Connector pinout compatibility 

2) Digital command and output data format 

3) Proof-of-concept RF circuit 
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3.2 Ground Transponders 
3.2.1 Introduction 

The power/transponder will provide a ground based 
measurement/datalogging capacity to ISAT-1. These transponders will have the 
capability of communicating with ISAT-1 to upload previously logged data. All 
remote ground transponders will be based upon the general power/transponder 
unit. 


3.2.2 Power/Transponder Unit 

The power/transponder unit will allow a common component to all 
configurations of the remote ground transponder. It will provide communications 
capability that is independent of the data to be uploaded to ISAT-1 and will 
provide the power required for communications. The environment is assumed to 
be harsh with a possible temperature range of -35°C to 50°C (-31°F to 122°F). 

The purpose of the ground transponder will be to uplink data gathered from 
ground sensors when the satellite comes in the coverage region. 


3.2.2. 1 Purpose (assuming FDMA) 

a. A receiver installed in the transponder would be constantly checking for a 
pilot signal (constant frequency tone) from the satellite. 

b. When it receives the signal, the transmitter would be switched ON and 
data coming from a data storage bank would be uplinked to the satellite. 

c. This would continue until the pilot signal level falls below certain 
threshold, indicating that the satellite has moved away from the coverage zone. 
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3.2.2.2 Data Rate / Bandwidth 


Average number of passes per day 

= 4 

Average contact time per pass 

= 8 min. 

Approximate number of sensors per transponder 

- 10 

Number of quantization levels 

= 1024 

Number of data bits per sample 

= 10 

Time between two samples 

= 15 min. 

Total data generated per day 

= 96 k bits 

Total data to be uplinked per pass 

- 24k bits 

Uplink data rate (assuming FDMA) 

= 24k/480 = 50 bps 

Bandwidth requirement per transponder 

= 50 Hz 

Frequency spacing between two transponders 

= 10 Hz 

Bandwidth requirement for 200 transponders 

= 12 k HZ 


3.2.2.3 Frequency Selection 

Selection of frequency is based upon the following considerations: 

a. Antenna size is smaller and thus easily manageable at higher frequencies. 

b. Free space loss is lesser at lower frequencies. 

c. Doppler frequency shift is higher at higher frequencies. 

d. Operating frequency of COMM- 1 and COMM-2 should be close enough 
to share a common antenna for on-board communication equipment. 

e. Ultimate frequency selection is subject to the availability of a particular 
frequency band and the FCC location. 


3.2.2.4 Reliability / Redundancy 

Since this is an experimental communication system and also the type of data 
is such that certain level of transmission errors can be tolerated, the maximum Bit- 
Error-Rate is 10' 5 which is lower than that of COMM-1 which was 10' 6 . 
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Also the transponder would remain on ground and hence the equipment 
redundancy is not critical. However since these transponders would operate as 
unmanned up-link stations under normal weather conditions, they are expected to 
have a mean-time between-failure of at least few months. 


3.2.2.5 Weight /Size 

Again since this is not an on-board equipment, the weight/size considerations 
are not very critical. The transponder size is expected to be around 10 lbs while the 
antenna weight/size is largely dependent upon frequency of operation with a lower 
limit of 10 lbs and an upper limit of 50 lbs. 


3.2.2.6 Antenna Design 

In order to minimize system complexity a fixed antenna with relatively wide 
beam width is considered suitable. The antenna polarization would be circular. 

Antenna type: Two Yagi Arrays arranged perpendicular to each other. 


Antenna Elements: Central supporting boom 

Reflectors 

Dipoles 

Directors 

Antenna Beam Width: About 100 degrees. 


1 

1 for each array 
1 for each array 
3 for each array 


3.2.2. 7 Power Consumption 

It is assumed that 1 10 V supply would be available for transponder's power 
requirements so power efficiency of the system is not very critical. 


Transmitter 

Receiver 


8-10 watts (only during transmission) 
2-4 watts (continuous) 
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3.2.2.8 Cost 

Ranging between $1000 to $2000. 


3.2.2.9 Suppliers 

a. Cynetics Corporation 
P.O. Box 2422 
3824 Jet Drive 
Rapid City, SD 57709 


b. Campbell Scientific 
P.O. Box 551 
815 W. 1800 N. 

Logan, Utah 84321 

(801) 750-9555 (Dave Meek) 


3.2.3 Weather Stations 

A primary application of the ground transponder is the automation of 
weather/meteorological data gathering. There are presently several companies that 
distribute weather and stock information to farmers around Iowa. The weather 
information provided by these companies are very general and cover very large 
areas in the map grid. More accurate weather information can be obtained at a 
higher cost than the already expensive standard weather package. 

One hundred fully equipped weather stations will provide automated data 
gathering around the state of Iowa. A Campbell Scientific Universal 10 meter 
tower with a wide range of sensors will provide the measurement/datalogging 
instrumentation. A standard IBM 9-pin serial connector will be used to interface 
with the common power/transponder unit. 
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These weather stations will provide the farming industry with weather data 
essential for crop management at a very low cost, if not free. The weather stations 
will also provide useful meteorological information to the Iowa State University 
meteorological center. 


3.2.3. 1 Weather Sensors 

A sorting process was made from the "Price List" Campbell Scientific 
catalog in order to determine which sensors were most adequate for farming 
purposes and which were most adequate for meteorological studies. 


3.2.3. 1.1 Farming Sensors 

Interviews were conducted with a couple of farmers and with an applications 
engineer at Campbell Scientific in order to determine which sensors would provide 
the most useful information for farming purposes. The following Campbell 
Scientific sensors were selected as "farming" sensors: 

1) Air temperature and relative humidity sensor. 

Model: HMP35C 

Description: Vaisala Temp. & RH Probe 

2) Leaf wetness sensor. 

Model: 237L 

Description: Wetness sensing grid 

3) Precipitation Sensor 
Model: TE525 

Description: Texas Electronics Tipping Bucket Rain Gage 

4) Soil moisture sensor 
Model: 227-L 

Description : Soil moisture block 
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5) Soil temperature sensor 
Model: 107B-L 

Description: Temperature probe 

6) Solar radiation sensor 
Model: LI200S-U 

Description: LI-COR Silicon Pyranometer 
Note: Most useful solar radiation sensor to farmers 

7) Wind speed and direction sensor 
Model: 05103-U 

Description: RM Young wind monitor 
Note: Most robust of all the wind sensors. 


3.2.3. 1.2 Meteorological Sensors 

A professor of the ISU Agricultural Department, Dr. Taylor, was 
interviewed to acquire information on the most useful sensors in meteorological 
studies. A recommendation was made to add the following sensors to the weather 
station in order to make it useful for meteorological studies: 


1) Barometric Pressure 
Model: PTA427 

Description: Vaisala barometric pressure sensor 

2) Snow Depth 
Model: UDG01 

Description: CSC ultrasonic snow depth sensor 
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3.2.3.2 Datalogger and Other Accessories 

Additional peripheral devices and extraneous hardware are needed to make 
the weather station operational. The following is a list of some of the accessories 
needed for the weather station: 

1) Ten meter tower, mounting base, grounding kit, anchors, and guy kit 

2) CR10 datalogger, extended temperature test, and keyboard / display 

3) Instrument enclosure and kit. 

4) Power supply and mounts. 

5) Datalogger software, interface and cables. 


3.2.4 Price List 


liUSSItW 



Model 

Description 

Price 

HMP35C 

Vaisala Temp. & RH Probe W/5'LD 

470.00 

UT12VA 

Solar radiation shield for HMP35C 

210.00 

05103-U 

RM Young wind monitor 

663.00 

019ALU 

Aluminum crossarm sensor mount 

70.00 

LI200S-U 

LI-COR Pyranometer W/34'LD 

220.00 

LI2003S 

Pyranometer base leveling fixture 

44.00 

025 

Pyranometer crossarm stand 

33.00 

237-L 

Wetness sensing grid W/50'LD 

60.50 

227-L 

Delmhorst soilmoist blk W/25'LD 

54.75 

TE525 

TE rain gage (0.01 "/tip) W/25'LD 

235.00 

107B-L 

Water/soil temperature probe W/50'LD 

67.00 

PTA427 

Vaisala barometric pressure sensor 

1,045.00 

UDG01 

CSC ultrasonic snow depth sensor 

800.00 


HMNhB 

UT930 

Universal 1 0M tower & adjustable mast 

550.00 

B18 

10M tower concrete mounting base 

80.00 

UTGND 

Grounding kit for 1 0M tower 

30.00 
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UTGUY 

Guy kit for 10M tower 

160.00 

UTEYE 

Eyebolt anchors for 1 0M tower 

30.00 



llllilliliill: 

CR10 

Measurement and control module W/CR10WP 

1,090.00 

XT 

Extended temperature test -55 to +85°C 

110.00 

CR10KD 

Keyboard/display 

275.00 

1 

12 x 14" Instrument enclosure 

180.00 

TN 

BKT kit for UT930 

25.00 




PC208 

Datalogger support software 

220.00 

SC32A 

Optically isolated RS232 interface 

145.00 

7026 

9-pin computer to 25-pin cable 

25.00 

ma&xst ^ '>■ • ^ - " ■ • - 


Subtotal 

ns** 

Less 5% educational discount 

-344.61 

Freight & insurance 

-160.00 

Total 

| 


Table 3.1 Weather Station Cost Estimate 


3.2.5 Future Work 

One fully equipped weather station should be purchased from Campbell 
scientific and tested in order to determine how to program the CR10 datalogger 
and to set data retrieval rates. The weather station can be easily connected to an 
IBM compatible computer through the 7026 9-pin cable. Tests on the ground 
transponder should also be made in order to determine the efficiency of the system 
as automated data retrieval. 


3.2.6 References 


(1) Campbell Scientific US. Price List. December 1991. 
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3.3 Primary CCD Cameras 


One of the most common uses of a low Earth orbiting satellite is as a 
photographic platform. Nothing seems to show off the accomplishment of a 
successful launch into orbit as well as a picture of the planet as seen from space. 
This will be the primary purpose of the Charge Coupled Device (CCD) cameras 
carried on board. They will give direct, unarguable evidence that ISAT-1 has 
achieved orbit. 

Prior to this semester, a trade study was performed to determine which of 
several models of CCD cameras would be best suited for the project. It 
recommended the purchase of the Panasonic GP-KR202 full motion model. While 
no doubt this is a high quality camera, the author does not believe it is the best 
camera to fly on ISAT-1 . Instead it will be shown that the higher resolution Canon 
Ci-20 still video system may be a better choice. 

A Charge Coupled Device (CCD) camera is, put most simply, a photon 
counting digital camera. The photographic surface is made up of thousands of 
picture elements (pixels), each of which is an individual photo diode. As photons 
of light strike each pixel, electrons are knocked from the surface and passed to a 
series of dynodes each of which amplifies the signal passed to the next. Finally the 
signal is collected and assigned a value based on the number and the quantum 
energy (or wavelength) of electrons collected. 


3.3.1 Camera Recommendations 

Most CCD arrays consist of over 250,000 pixels. As was mentioned 
above, each pixel is assigned of value based on the photons it receives. For a grey 
scale camera there are usually 8 bits per pixel (or 2^ = 256 levels). Therefore one 
grey scale picture will require ~ 0.25 MBytes of storage. For a full color picture 
24 MBytes are required, (8 for each of red, blue and green), therefore tripling the 
required storage. Fortunately, digital information has the quality that it can be 
easily manipulated. If the detector can only discern 300 individual colors from its 
palette of 16.8 million colors, then the MBytes computer can compress the data 
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file by using only the color values present. Other techniques such as pixel 
averaging can also reduce storage space required, but they degrade image 
resolution. 

Even with data compression techniques, it is unlikely that the 9600 MBytes 
downlink rate will be enough to transmit more than one or two pictures per pass 
making full motion video transmission impossible even by storing frames for later 
transmission. This is the primary reason for selecting a still video system rather 
than the motion video camera. If the latter were chosen, it would require some 
type of frame grabbing hardware and software to pick out an individual frame. 
The still camera, however, outputs one frame of digital information directly into 
the flight computer for either on board processing or downlinking. 

With the decision made to use a still camera, the possibility opens up of 
using a hi resolution color camera that will create enormous data files for each 
picture taken and then downlinking that image information in segments over 
several passes, (another benefit of digital technology). 

The camera chosen, the Canon Ci-20, was the only industrial grade, still 
video model for which information was available. Table 3.2. 


Number of pixels 

380,000 pixels (774Hx488V) 


20-lux 

Power Supply 

9V, DC 

Power Consumption 

5W Max 

Operating Temperature 

-10 to 45 C (14 to 113 F) 

Dimensions (WxHxD) 

71x58x90 mm 

Weight 

Approx. 390 g (13.7 oz) 

Minimum exposure time 

1 ms 


Table 3.2 Technical Specifications for Canon Ci-20 


It should be noted that, while the 5W power consumption is more than the 
Panasonic model, the camera need only be turned for a few minutes around the 
time of imaging. 
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It is therefore the author's recommendation that ISAT-1 carry two Canon 
Ci-20 still video color CCD cameras. The two cameras will share the satellite's 4 
inch Maksutov-Cassegrain telescope. In addition to the telescope, one camera will 
be fitted with additional optics to increase the surface resolution. The optical 
system is being designed by Gary Cameron of Iowa Scientific Optical. 


3.3.2 Data and Resolution 

Some preliminary calculations have been made for data storage 
requirements and image downlink time as well as first estimates of surface 
resolution without the use of the telescope. 

First of all, assuming a full 24 bit per pixel image with the Ci-20's 
774Hx488V pixel array, one image will require 1.1 MByte of storage. At a 
downlink rate of 9600 bits per second, this will take just under 16 minutes which is 
equivalent to 3-4 satellite passes. 

Next, it was assumed that the camera had an angular resolution of 10 
micro-radians per pixel (commonly found in text book research but not listed in 
technical specifications) and the satellite at an orbit altitude of 600 km. Doing this, 
the surface resolution is calculated to be 6 m without use of the telescope. 

There is some concern as to whether the velocity of the satellite will have a 
large impact on image quality. Again assuming a circular orbit with an altitude of 
600 km, the orbit velocity is determined to be 7.56 km/s. Also, it is assumed that 
the minimum exposure time of 1 ms is used and the satellite is continuously 
pointing toward the center of the Earth. With this, the satellite covers 6.9 m of 
ground distance during the exposure time which is about the same as the per pixel 
resolution. Further calculations will need to made in designing the telescope. 
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3.3.3 Future Work 

First and foremost, a price needs to be found for the Canon model. If it is 
deemed too costly, it may be necessary to research other manufacturers. It is still 
the author's recommendation that a still video system be used. Other specifications 
which need to be determined are the angular resolution in micro-radians per pixel 
(from which surface resolution and field of view can be calculated. These will help 
determine the image quality and assist the design of the telescope. 


3.4 GPS Unit 

One of the hottest topics in astrodynamics and navigation is the 
NAVSTAR / Global Positioning System (GPS). The system consists of a 
constellation of 24 satellites arranged in six orbital planes. Using coded signals 
from any four satellites, an on board GPS receiver can either calculate its position 
on orbit at a precise time or downlink the information for computation on the 
ground. Knowing the precise position of the satellite containing the receiver at 
various times provides an efficient means of orbit determination which can be 
monitored by the ISAT Company rather than having an outside source track the 
satellite. 

The author consulted Professor Peter MacDoran from the University of 
Colorado at Boulder. Mr. MacDoran holds a number of patents in small GPS 
receiver technology. He was kind enough to send some background information 
since the author was too ill equipped to even ask the technical questions required 
in choosing GPS receivers and antennas. 

It was learned, however, that GPS units can be made as small as 0.75 
grams with a power consumption of 0.3 W. These, though, are the type used on 
meteorological radiosonde balloons and therefore are probably not suited to the 
rigors of space flight. The antennas are small domes measuring 2-3 inches in 
diameter. Mr. MacDoran also gave his assurance that the data stream from the 
GPS receiver to the ground station is minimal. 
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For purposes of national security, errors are built into the GPS satellite 
signal. Therefore, it is only possible to obtain an accuracy on the order of tens of 
meters which is probably accurate enough for ISAT’s purpose. However, using 
differential techniques (DGPS) it is possible to be accurate to within centimeters. 

However, just as with the personal computer industry, big names mean big 
prices. Fortunately, as with computers, there are smaller companies which 
produce the same high quality hardware at more reasonable prices. Large 
companies such as Motorola currently manufacture flight ready hardware. 
Professor MacDoran has expressed some interest helping the ISAT project 
through the small business he operates in GPS technology. GPS is certainly the 
wave of the future and it is the author's strongest recommendation that this study 
be pursued further. 


3.5 Satellite Laser Ranging 

Satellite Laser Ranging (SLR), is a relatively new means of tracking earth 
orbiting satellites. It consists of two main components. First, the satellite must 
have some sort of reflector on it, either a mirror or a comer cube lens. Second is a 
ground based laser unit and signal receiver. If it is known approximately where 
and when the satellite is going to pass over, a laser beam can scan a portion of the 
sky. As the satellite passes through the scan region, the laser light reflects from 
the mirrors and returns to a point near the source. By measuring the time required 
for the light to travel to the satellite and back, its range can be found from which 
the orbit may be determined. The author is by no means an expert and suggests 
that future workers do some research to better understand methodology. 

Very little had been done prior to this semester in determining the 
feasibility of SLR on a satellite as small as ISAT-1. In fact all that had been done 
was obtaining a catalog from Rolyn Optics. Dr. Jeyapalan of the Civil Engineering 
Department here at Iowa State University was of the opinion that higher quality 
reflectors would be necessary than those in the catalog. 
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The advantage of SLR over GPS tracking lies in the non-mechanical nature 
of the on board equipment. If the electronics in a GPS unit fail, there is no 
recourse but to have the military track the satellite. With SLR, however, there are 
no electronics which may fail and absolutely no power consumption on the 
satellite. The drawback is cost. Constructing a laser and operating it continuously 
will cost several million dollars which the Iowa Satellite Company can not yet 
afford. Another possibility is to have some other group such as the University of 
Texas at Austin, continuously track the satellite. This will still cost several 
hundreds of thousands of dollars. The third and most financially attractive option 
is to have it ranged only on occasion. This can be done at an annual cost -10-30 
thousand dollars. 

Clearly, much more work needs to be done in determining the feasibility of 
this segment of the payload. The first thing which needs to be determined is 
whether the surface area on the bottom of the satellite, (around the telescope 
opening), is enough to even reflect a discernible signal. Next, it would be 
beneficial to know how large of a laser would be required to obtain a measurable 
signal. 


3.6 Micro-Meteor Detection 
3.6.1 Introduction 

In the space environment where the satellite will be orbiting, there will be 
much space debris. Because of the size of the satellite, micrometeorites, meteors 
less than 10 mm in length, will be the main concern. In the earth's orbit right now, 
there are at least 3.5 million particles of this size. A system needs to be devised to 
detect these micrometeor collisions and locate where on the satellite the particle 
actually hit. From this an analysis of the possible damage to systems can be 
ascertained. For additional infomation on space debris see Appendix C.3. 
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3.6.2 Design 

A catalog from Kistler Piezo-Instrumentation was obtained previously, and it 
was this catalog where information on accelerometers was found. In order to 
properly size the accelerometer, the g-force that the accelerometer would see in 
one microsecond was needed. This is a standard way for comparing 
accelerometers. By applying the conservation of mass and momentum equations 
to the satellite and a particle, heading towards each other with equal velocities, the 
one microsecond impulse was calculated to be about 20,000 g's. This narrowed 
the process to choosing a high shock accelerometer. Taking into consideration the 
temperature range involved during orbit as well as size and weight, as preliminary 
accelerometer was chosen. In addition, charge amplifiers are needed to run the 
accelerometers. The smallest charge amplifier compatible was chosen. 

Information about the accelerometer, a Kistler high impedance 8044 
accelerometer. Table 3.3; and the charge amplifier, Kistler Model 5037A, Table 
3.4. 


Size 

0.74"length x 0.425" diameter (cylinder) 

Mass 

7 grams 

Power 

see Charge Amplifier (Table 3.4) 

Cost 

530.00 each (24 required) 


Table 3.3 Kistler 8044 High Impedance Accelerometer 


Size 

4.7" wide x 1.3" high x 2.9" deep 

Mass 

358 grams 

Power 

0.06 watts 

Cost 

$800.00 (4 required) 


Table 3.4 Kistler 5037A Multi-Channel Charge Amplifier 
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The charge amplifier only has six inputs on the three channel model, so for 
the current configuration of 24 accelerometers, four charge amplifiers would be 
needed for a total power cost of 0.24 Watts. 

3.6.3 Future Work 

This has become a preliminary design. Kistler has just come out with a new 
line of accelerometers and accessories that are designed for the space environment. 
The people at Kistler strongly recommend looking into these products before using 
the previously stated ones, and I must agree. The information has not come yet, 
but they should be smaller and stronger than those in the current catalog. In 
addition, a configuration for the accelerometers needs to be developed to properly 
locate the impacts of the micrometeorites. Minimal redundancy should be 
considered to lower the cost and size constraints of this subsystem. 


3.6.4 Appendix / Theory 

In order to locate the collision, a series of equations must be solved. A 
summation of forces in each direction is first performed, as well as a summation of 
moments in each direction. These equations need to be related to the center of 
gravity, a known location. From this information and the fact that the force times 
the distance equals the moment, the impact location may be found. The equations 
needed are listed below. 

Z F * - m c* U x Z F y = m Z F * = m * u z Et l n 51 


M 

2 

X 

II 

I w - (I -I ) w w = F d 

A XX X V A yy A zz / T ’y z X X 

Eqn. 5.2 

Z My = 

Iw - (I -I )ww = F d 

yy y V zz xx / z x y y 

Eqn. 5.3 

Z = 

Iw - (I -I ) w w = F d 

A zz z V A xx A yy * x y z ^ z 

Eqn. 5.4 

KNOWN: 

meg, I xx, Iyy, Izz, Wx, Wy, Wz 


MEASURED: 

Ux, Uy, Uz 


UNKNOWN: 

dx, dy, dz 
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3.7 Life Sciences 

3.7.1 Introduction 

The life sciences project was desired to provide students ranging from 
kindergarten through twelfth grade the opportunity to use an interactive 
experiment. 


3.7.2 Criteria 

The experiment needs to be simple. There are two reasons for this, the 
first being the age level of the students and the second being the limited space 
available. Dr. Misra of the I.S.U. Seed Sciences Department has been given a 
copy of all data sent to ISAT dealing with the life sciences project. 


3.7.3 Research 


Many government agencies were contacted to determine if there was data 
available to compare with any experimental data we might receive. The 
application of any experiments that might be adaptable for our satellite were also 
examined. 

After checking with Goddard, Ames, and Johnson Space Centers (all of 
which deal more with the human size life sciences), the name of Dr. Knott from 
Kennedy Space Center provided some useful information. Dr. Knott sent 
summaries of the experiments that he thought might be adaptable to our 
specifications, but on review, these seemed slightly more difficult than originally 
thought. Dr. Knott also gave the names of two others that were of great help. 

Dr. Brown from Kennedy is a plant physiologist had suggested the use of 
increased levels of biological activity to determine germination or plant growth. 
He also suggested the possibility of attaching a transducer to the seed to measure 
pressure changes in the hull of the seed. 
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Dr. Guikema from Kansas State University was more than helpful. He sent 
a package of detailed information which K-State has collected from their 
experiments. He suggested using one of two methods in which to get the nutrients 
to the seeds. 

1 . a rotating well using two cylinders with off-set wells which allow 
the mixing of seeds and nutrients when aligned. 

2. a syringe with multiple chambers separated by a membrane which is 
punctured when the plunger is depressed, thus mixing the seeds 
and nutrients, Figures A. 3.1 and A. 3. 2. 

After studying the two methods for seed/nutrient mixing, it has been 
decided to attempt to use the syringe with two chambers. This was decided for 
multiple reasons, two being that the rotating well was to big to fit the allowable 
space and the problem of not easily determining the germination or growth of the 
seed. 


The syringe is produced by BioServe Technologies in Boulder, Colorado. 
When contacted. Dr. Stodieck sent detailed design information describing the 
syringe. Dr. Stodieck mentioned that they could possibly produce the syringe to 
our specifications for the cost of the materials, and the sharing of the data 
collected. 


3.7.4 Suggestions 

By utilizing the two large payload boxes, two mini-micro-cams (Figure 
A. 1.14), and smaller syringes, it should be possible to have many individual 
experiments. Two suggestions to depress the plunger are using compressed gas or 
a spring. 

For the compressed gas, the exact amount to obtain the right pressure must 
be calculated as well as the release mechanism. Also the necessity of a light source 
will reduce the "experiment" space. 
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The spring idea would need some release mechanism as well, but that could 
possibly come from a glow-plug. This glow-plug could also be the light source. 

A second experiment could be run using brine shrimp. When the eggs are 
dried, they are capable of lasting indefinitely. Then when a saline solution is 
administered, the eggs hatch. The shrimp would have an advantage over the plants 
if visual data is sent back, because they would swim around for a couple of days to 
a week unlike the plants which would simply germinate and die in a matter of a 
couple of days. 


3.7.5 Future Work 

Dr. Misra has been given a copy of what has been done and what has been 
suggested and he sounds very enthusiastic about the project. He is going to look 
into the research that has already been done to make sure that we can accomplish 
it also. 


A working model of the syringe needs to be obtained and specific 
dimensions need to be determined. Dr. Stodieck said he would try to answer any 
questions about the syringe or is construction. Dr. Guikema also would be 
interested in any questions or comments about either the experiments or the 
hardware we have chosen. 

If it is found that these decisions will not work, the next step would be to 
try to contact Dr. Alstin of Parks Seed Co. in South Carolina. Dr. Alstin helped 
with LDEF and may possibly have another experiment which could be used on 
ISAT. Numerous attempts have been made to contact Dr. Alstin, but he is a busy 
man and very difficult to catch. 
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3.7.6 Phone Numbers 


Dr. Misra 
Dr. Guikema 
Dr. Stodieck 
Dr. Knott 
Dr. Brown 
Dr. Alstin 


ISU Seed Sciences 
Kansas State University 
BioServe Technologies 
Kennedy Space Center 
Kennedy Space Center 
Parks Seed Co. 


(515) 294-6821 
(913) 532-6011 
(303)492-4010 
(407) 853-5142 
(407) 853-3165 
(803) 941-4445 


Chapter 4 
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Chapter 4 Power Systems 
4.1 Introduction 

Designing the power system for a large satellite is a relatively easy task, 
unfortunately since our satellite is a small one the power system design is not as 
simple. 


With a limited amount of surface area on which to attach solar cells the 
power which can be generated is limited. With this low power generation, battery 
design is a very crucial part of the overall process. 

There were many problems encountered in the design process of this 
system. The rest of this section describes the steps and problems encountered in 
designing the power system for ISAT-1 . 


4.2 Solar Array Design 

Power analysis for the satellite consisted of first defining the power needed 
by all of the subsystems onboard. Then with a given solar array size the minimum 
efficiency of the cells was determined. However this efficiency was out of the 
range of current technology, so a different approach was used. 

First, with a given outer surface area it was decided that the solar cells 
must cover all of the available area on the satellite. This ruled out earlier concepts 
where some of the surface area of the satellite was to be used for a video display 
or sign of some sort. Additional information may be found in Appendix C.4. 

The solar cells decided on are to be Gallium Arsenide. These cells are 18% 
efficient at beginning of life and greater than 1 5% efficient at end of life. The cells 
are to be arranged on panels, with each panel consisting of 32 cells in the 
arrangement shown in Figure 4.1. The satellite will have four panels per side on 
five of the sides and 2 panels on the sixth side. These two configurations can be 
seen in Figures 4.2 and 4.3 respectively. 
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Panel Configuration 
32 Cells / Panel 
Volatage 32V 

Current 225mA 

Power 6W 
Cell: 

GaA 

> 15 % 

20 x 40 mm 
series 


Figure 4.1 Solar Panel Configuration 


Main Panel 
Configuration 
4 Panels / Side 



Figure 4.2 Main Solar Array 
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Bracket Panel 
Configuration 
2 Panels / Side 



? 

■ 


Figure 4.3 Bracket Solar Array 


4.3 Power Regulation and Control 

The power system for most large satellites runs by drawing power from the 
solar cells most of the time and only draining the batteries when the amount of 
power consumption is high. This configuration will not work for a small satellite 
since the solar array on a small satellite puts out much less power. So the 
approach used in small satellites is to draw power from the batteries all the time 
and when the solar array can, it will charge the batteries. This is the setup of the 
system chosen for ISAT-1. 

The output from the solar cells is monitored and controlled by an electronic 
device that finds the maximum power point for the array and operates the array at 
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that voltage and current. This power is then used to charge the batteries. After 
the batteries there is a power regulation unit. This unit has 12, +5, -5 volt outputs 
and a ground. After the power regulation unit there is a power distribution unit. 
This unit will be controlled by the flight computer and will control all the power to 
all the subsystems. It will have current meters on all the outputs to keep track of 
how much power each subsystem is using. This will allow the flight computer and 
ground systems to keep track of the power being used and avoid discharging the 
batteries past their limit. 


4.4 Batteries 

For a small satellite the selection of batteries is crucial to the performance 
of the power system. The batteries need to be able to go to low depths of 
discharge without damaging the battery. They must also be light weight and 
reliable. There are two kinds of batteries currently being used in satellite 
applications Nickel-Cadmium and Nickel-Hydrogen batteries. Typical NiCd 
batteries have an energy density of less than 20 Whr/kg and allow a moderate 
depth of discharge. NiH2 batteries have an energy density of 28 Whr/kg and allow 
for deeper depths of discharge. 

The batteries will be drawn upon for power at all times so the battery must 
be able to go to deep depths of discharge. This consideration and the need to 
keep weight a minimum necessitated the use of the NiH2 batteries, even though 
they are more costly and less common than NiCd batteries (O'Donnell, p. 46). 


4.5 Final Configuration 

The proposed final configuration is as follows: 

Power outputs at 12V, -5V, +5V 
Maximum power per orbit of 46.8 Watts 
Minimum power per orbit of 33. 1 Watts 
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The amount of power required currently by the subsystems can be seen in 
Table 4.1. From this it can easily be seen that more power is needed than can be 
generated if all of the systems are to be running at one time. Fortunately this is not 
the case. With proper power management there should be enough power to power 
all the subsystems. 


# 

Component 

Power (W) 

1 

Structure 

0 

1 

G. G. Boom 

0 

2 

Communications 

7 

2 

Flight Computer 

7 

2 

Large Experiment Module 

4 

2 

Medium Experiment Module 

1.5 

6 

Small Experiment Module 

1.1 

1 

Pitch Torqrod 

0.25 

1 

Roll Torqrod 

0.25 

1 

Yaw Torqrod 

2 

2 

G. G. Boom Tanks 

0 

1 

Earth Radio Frequency Experiment 

7 

2 

CCD Camera 

5 

1 

Telescope 

0 

6 

Battery 

0 

6 

Battery Mount 

0 

2 

Magnetometer 

0.04 


Table 4.1 Power Budget 
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4.6 Future Work 

To get more power it may be possible to place solar cells under the 
mounting bracket. This will give an increase of 4 watts in the low power end of 
the orbit as seen in Figure 4.5. The original power generated without placing solar 
cells under the mounting bracket may be seen in Figure 4.4. This addition of solar 
cells and careful power management should allow for all the subsystems to have 
enough power. 



Figure 4.4 Solar Array Power 
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Inclination Angle (deg) 


Incresed Normal 

Figure 4.5 Improved Solar Array Power 


Things that still need to be examined include the cost of all the components 
of the system as well as deciding if off the shelf equipment will be used for the 
power regulator and controlers or if these components will need to be built by the 
ISAT team. 


4.7 References 

O'Donnell, Patricia M. "Space Batteries for Mobile Battlefield Power 
Applications" EEEE AES System Magazine, December 1991. 
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Chapter 5 Thermal Considerations 
5.1 Introduction 

The purpose of a thermal-control subsystem is to maintain all the sections 
and other subsystems of the spacecraft within their specific temperature limits. 
This includes all mission phases from takeoff to the death of the satellite. A 
satellite experiences two major environments direct solar heating and solar eclipse 
while it is in orbit around the Earth. The thermal environment changes 
dramatically for these environments and the thermal subsystem must be able to 
compensate for these changes. Other heat inputs to the spacecraft environment 
include Earth reflected solar albedo, Earth infrared, and electrical components on 
board the satellite. For the case of the ISAT-1 all of these have been taken into 
consideration since they will all affect the thermal environment of the satellite. All 
of these inputs vary highly with time and it is necessary for the thermal-control 
system to be able to maintain a working thermal environment for the satellite and 
all of its components. Although the thermal-control subsystem is one of the most 
important systems it usually only makes up about two to five percent of the total 
cost and about the same percentage of the dry weight. 

The work that had already been done, was just some basic ideas about 
what was needed for a thermal -control system for the ISAT-1. Some things 
considered were such ideas about special coatings for the satellite's skin, heaters, 
multilayer blankets, and a space radiator that radiate waste heat into deep space. 
The bulk of the work on the thermal subsystem has been done by the Mechanical 
Engineering Department at the University of Iowa. This coordination of efforts 
has turned out to be very educational for both parties of interest. The students at 
the University of Iowa who were taking a thermal design class got to actually 
work on something that is going to be reality and it also gave the students in 
Aerospace Engineering at Iowa State University a chance to learn a little more 
about thermal aspects in design. 
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5.2 Background 

The spacecraft environment for the satellite was modeled as four types of 
radiation. These radiations included solar, solar albedo. Earth infrared, and 
satellite emission and their values are given in Table 5.1. The mean temperature of 
space was modeled as 4 K. These radiations all affect the thermal environment for 
ISAT-1 and the thermal-control system to be designed will have to consider all the 
types and compensate for them. The spacecraft's environment is a harsh one with 
very transient temperatures and a thermal-control system is the only thing 
protecting the delicate payloads from this harsh environment. 


Spacecraft Environment 

-Radiation Solar = 1430 W/m A 2 

- Tspace = 4 K 

Albedo = 0.3(Solar) 


Earth Infrared = 230 W/m A 2 


Satellite Emission = eoAT A 4 


Table 5.1 Spacecraft Environment 


5.3 Current Work 


As a first approximation some assumptions were made to make the first 
estimates of the thermal environment an easier task. The satellite will be spinning 
about the centroidal axis so it was modeled as a cylinder. The total internal power 
generation was estimated to be around 40 watts. The mass of the bus was 
estimated at 15 kg of aluminum. The transient response was modeled by a 
mathematical expression relating all radiations presented earlier. A FORTRAN 
program was then written using these assumptions which were used to specify the 
thermal environment more accurately and account for the transient response the 
satellite will encounter. A finite difference method was also added to this program 
to give even more accurate information about the thermal details. The solar 
radiation for the program was modeled as a step function having its maximum 
value stated above when in direct contact and going to zero in eclipse. The albedo 
radiation had to be modified a little to account for the dusk and dawn periods that 
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would affect its intensity. It was modeled as constant during daylight, a ramp 
function during dusk and dawn and zero at eclipse. 

Through the use of this program the skin temperature transient response 
was determined. It was shown that after about five orbits from the initial orbit the 
fluctuations in the response died out. This left the transient response the satellite 
would encounter for the rest of its life. High skin temperatures were on the order 
of about 330 K and low skin temperatures were seen to be about 280 K after the 
fluctuations died out. A sensitivity analysis was also performed using the program. 
The length of the satellite was varied while holding all other properties constant. 
The results showed that for increases in length the mean temperature of the skin 
went up fairly dramatically. 

Once a nominal orbit environment was established it was necessary to 
upgrade the math model to account for the internal analysis of the satellite. In this 
upgraded model the satellite was broken up into six sections consisting of top, 
bottom, sides, and three levels of shelving. Thermal resistance values between 
these sections were set up in order to isolate certain sections from others. 
Temperature tolerances between these sections were estimated and the thermal 
resistance values solved for. Using these values, insulation can be designed and 
thermal coatings to the exterior and interior of the satellite can be looked at. The 
insulation and coatings will keep the internal thermal levels within limits set by the 
other components of the satellite. It is hoped that this will be all that is necessary 
and that no radiators or heaters will be needed. Since ISAT-1 is a small enough 
satellite and does not contain such things that require very strict limits this may be 
possible with the exception of the robotic arm camera which will probably need 
some sort of special thermal control. Figure 5.1 shows the upgraded math model. 
A more detailed report may be found in Appendix C.5. 
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Upgraded Math Model 

Six Sections of Satellite 
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Level 3 
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Level 2 
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Level 1 
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3 mm 

Level 3 

A^ 

R8 
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3 

Level 1 



R = 

Thermal Resistivity 



Figure 5.1 Thermal Model 


5.4 Future Work 

Future work for the thermal-control subsystem depends on the final choices 
for the other subsystems and payload components. The heat outputs and 
temperature limits must be specified before design work can continue. When all of 
these are specified the design work can become more detailed and final decisions 
made. At the moment not all choices for payload and bus designs are final. The 
continued upgrading of the math model and results from this model are to be dealt 
with in future work. Special thermal coatings and insulation are being discussed 
and will be recommended at a point in which all component choices have been 
made final. Also, the final decisions about heaters and radiators will depend on 
these choices. When all choices have been made possible problem areas can be 
considered. 

Another important aspect of the future work will be deciding what thermal 
tests will have to be performed before launch. A testing schedule will be 
developed and testing cycles decided upon. These tests will include component- 
level tests and system-level tests. The testing process is very important to the 
dependability of the thermal-control system. For example, multilayer blankets are 
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affected quite considerably by overlapping at joints and their thermal performance 
degrades quite severely. They are also very sensitive to being handled and installed 
since smashing them could ruin their thermal productivity. Also, test times in 
vacuum chambers must be long enough for the air trapped within the multilayer 
blankets to vent. Otherwise the blanket will not perform as it would in space. 
This is why the testing of the thermal system is so very important. If anything gets 
by the testing the whole satellite project could be jeopardized. 


5.5 References 


Larson, Wiley J., Wertz, James R. 1992.. Space Mission Analysis and Design. 
Second Edition. Kluwer Academic Publishers. 409-430. 
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Chapter 6 Attitude and Control 
6.1 Introduction 

The attitude determination and control subsystem stabilizes the vehicle in 
the correct orientation for the mission. This requires both accurate knowledge of 
the spacecraft attitude and 3-axis control to reject disturbances. 

Sensors use external references such as the sun, the earth's horizon, the 
earth's magnetic field direction and the stars to correctly determine the vehicle's 
attitude. Control techniques may include: gravity-gradient control, momentum 
wheels, reaction wheels, control moment gyros, thrusters, or magnetic torquers. 


6.2 Overview For ISAT-1 

Three scenarios must be considered for the control of ISAT-1. These are: 
despin, normal operation and inverted orientation. Despin includes the damping of 
initial spin rates from orbit insertion. For normal operation, a more precise attitude 
is required and disturbances must be rejected. Inverted orientation is a special 
problem which must be addressed due to the use of the gravity-gradient boom. 

A block diagram for ISAT-1 is shown in Figure 6.1. The disturbance 
torque's which may be encountered include: third body effects, oblateness, hits, 
solar radiation, residual dipole and atmospheric drag. All disturbances except the 
last two are considered negligible for this vehicle. For each of the three scenarios 
mentioned, a controller and feedback gain must be determined. 
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Disturbance Torques: 

1) hits 

2) oblateness 

3) 3rd body effects 

4) solar radiation 

5) residual dipole 

6) atmospheric drag 

For each scenario the controller and the feedback gain must be designed. 

Figure 6.1 Control Block Diagram 

6.3 Gravity-Gradient Control 

For ISAT-1 a passive control technique called gravity-gradient control will 
be utilized. This method requires neither power supply nor attitude sensors. 
Instead, the inertial properties of the vehicle will determine its stability. 

A body subjected to gravitational force naturally aligns its axis of minimum 
inertia along that line of force Therefore, by attaching a boom onto ISAT-1, its 
inertial properties can be altered so that it is inherently stable in an earth-facing 
attitude. 

It has been found in historical data that s x = .98 provides the desired 
stability, s x is defined in Eqn. 6.1. On the other hand, Sy = j (Eqn. 6.2) should be 

avoided because it would cause resonance. 



Eqn. 6.1 
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s y- 


Iy - I7 


Eqn. 6.2 


6.4 Magnetic Torquers 

Since gravity-gradient control causes oscillations and does not stabilize the 
local yaw axis, another method of control is needed. 

Magnetic torquers provide an inexpensive method of three-axis control 
with no moving parts and relatively low power consumption. Magnetic coils on 
rods generate dipoles and produce torque proportional to the Earth's changing 
magnetic field. 


6.5 Sensors 

The attitude sensors chosen for ISAT-1 are a sun sensor and a 
magnetometer. These are light-weight, inexpensive, require little power, and are 
sufficiently accurate for ISAT-l's mission. Sun sensors detect the angle of incident 
light to determine orientation. Magnetometers measure the direction and 
magnitude of the Earth's magnetic field and are used in conjunction with the torque 
rods. 


6.6 Despin Model 

In order to solve the problem of attitude control for ISAT-1 a computer 
simulation has been designed. For the despin scenario, the satellite has been 
modeled as a 50 kg homogeneous mass before boom deployment. The orbit is 
assumed circular with a semi-major axis of 7178 km and an inclination of 50 
degrees. 

A bang-bang control law was used to model the torque rods on ISAT-1. 
The torque (T) applied to the vehicle is equal to the saturation moment (m sa t) 


Chapter 6 Attitude and Control 


70 


created by the torque rod and is applied in the opposite direction of the changing 
db 

magnetic field (“jp. 

T = -m sa t * sign Qp Eqn. 6.3 

A spherical harmonic Earth magnetic field was used and the saturation 
dipole of the torque rod was set a 1.0 A*m 2 this creates a maximum possible 
torque of 4.3(1 0'^) N*m on the vehicle. 

Using rigid body dynamics and two-body orbital mechanics with the one 
primary being Earth, test runs were performed. The program is included in 
Appendix B.6. An initial spin rate was given to each body axis and the despin 
control history was shown graphically in Figure 6.2. For each of these cases, all 
disturbance torque's were neglected and perfect attitude knowledge was assumed. 
As can be seen in Figure 6.2, reasonable initial spin rates can be brought under 
control after only a couple of orbits. 



Figure 6.2 Despin Control 
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6.7 Capture 

After the despin has been initiated, the next important procedure is to 
determine exactly when the gravity gradient boom should be deployed. The object 
is to have the right conditions so that the satellite ends up in the proper orientation 
with respect to the Earth. In order to accomplish this, a relationship between pitch 
rate ( q. ) and pitch attitude (q) needed to be found. Equation Eqn. 6.4 is the 
resulting equation from the original pitch equation and is the relation between q 
and (q). 


,2 = 


— ! K 

3 w 0 2 


Syo 


+ sin 2 (q 0 ) 


Eqn. 6.4 


Assuming that there are no external torque's acting on the body, and if (a) 
is held constant and (q) and (q) are varied from to ^ . From this 

information, a stable range for boom deployment can be found around (q) equal to 
0 and a < 1 . The corresponding (q) quantities can then be determined. 


6.8 Future Work 

Since the control of ISAT-1 is completely dependent upon attitude 
knowledge, the modeling of sensors is very important. A complete understanding 
of the sensors is needed, and a method of filtering the data must be developed. 

The effects of residual dipoles and atmospheric drag must be added to the 
computer simulation and analysis. It must first be determined what will be the 
residual dipoles that are caused by internal electronics in the satellite. This 
required knowledge of all payload systems. 

The scenarios of normal operation and inverted orientation must also be 
analyzed. In normal operation, a method of dampening precession and nutation 
must be found. For inverted orientation, the problem of flipping the spacecraft 
back to a normal orientation must be solved. 
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Chapter 7 Ground Systems and Operations 
7.1 Introduction 


Ground Systems and Operations cover a wide field of study. It deals with 
the planing of all the phases of the Iowa Satellite Project, starting from the 
conceptual design until ISAT-1 plummets into Earth's atmosphere. First, it is 
responsible for the compilation of detailed time lines, both overall and a more 
detailed test schedule. Second, systems for the ground station will be selected. 
Finally, flight operations and procedures must be decided on. 


7.2 Time Lines 

The time lines will provide a basis for an orderly development of ISAT-1. 
Time lines must be made for many aspects of the project to give designers and 
operators a sense of direction. This spring semester, the focus of operations was 
on a revision of the overall and more detailed test time lines. First, the overall time 
line maps the phases of development from preliminary design to sustained 
operations. Second, the test time line was revised to reflect changes in the overall 
time line and more detailed input from design groups. 


7.2.1 Overall Time Line 

Now knowing that it is possible to have a late 1995 launch it was necessary 
to revise the overall schedule (since the Summer 1992 report was based on a late 
1994 launch). 
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7.2. 1.1 Preliminary Design 

This phase has primarily been completed. This includes the development of 
the overall concept of the satellite. The payload packages were selected. The 
communication frequencies have been requested. 


7.2. 1.2 Detailed Design 

This semester has been focussed on this phase. Knowing the conceptual 
design, this phase concerned the design and selection of hardware that will 
compose each subsystem. Specific subsystem configurations will be drawn, down 
to the wiring. Ground system configuration will also be decided . 


7.2. 1.3 Construction 

This phase will be the actual fabrication of the satellite. Components will 
be purchased and made for each subsystem and hardware for the ground station 
will be purchased. 


7.2. 1.4 Testing 

Every subsystem, system and component must be tested to determine if 
that unit will operate in orbit. There will be four levels of testing. 


7.2.1.4.1 Proof of Concept 

The proof of concept test will determine if a certain subsystem will be 
feasible and will provide the data that is desired. For example, the Meteoroid 
Location package could utilize a computer code. This code would be used to 
determine the position of particle impact from acceleration values in certain 
locations. 
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7.2.1.4.2 Acceptance Tests 

From the detailed design, a list of components that comprise each 
subsystem, will have been compiled. This detail design will govern the operating 
parameters of each component, while the acceptance test will determine if these 
parameters are met. Again for the Meteoroid package, the accelerometers should 
be tested to determine if they meet the factory specifications (since they will be 
purchased). 


7.2. 1.4.3 Integrated Tests 

The Integrated tests will consist of connecting all the components of a 
particular subsystem and determine if they work together. A possible test for the 
Meteoroid location package could utilize the mockup of the satellite. A number of 
accelerometers will be attached to the structure and connected to an outside 
computer at which time a simulated meteoroid impact will be conducted (possibly 
a bb gun). Then the computer will record and process the acceleration values and 
determine the location. 


7.2.1.4.4 End-to-End Tests 

Finally, end-to-end testing will be done. This will test the satellite, which is 
fabricated and put together as it would be in orbit, to determine if all the 
subsystems will work together (even the communications link will be used). Here, 
the meteoroid location payload will be operated in tandem with the other 
subsystems in order to determine if this payload will operate adversely in 
conjunction with its neighboring subsystems (such as electromagnetic 
interference). 
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7.2. 1.5 Training 

This phase will provide future operators working knowledge of the 
satellite. This will be a combination of lecture, workbook and hands on. For 
example, the operators must know how to operate the robotic arm without 
affecting the satellite's stability. 


7.2.1.6 Shipment 

This phase will transport the satellite to the launch site. A method of 
transportation must be found that will not damage the satellite. 


7.2.1.7 Pre-launch Tests 

These tests will confirm the operational status of the satellite. This phase 
will determine if damage was incurred to the satellite during transportation or 
some other source. This will also include on pad tests, if at all possible. 


7.2.1.8 Launch 

Here the satellite's health will be observed during launch. 


7.2.1.9 Initial Ops 

Following launch, certain events must be performed to make the satellite 
operational. Also a series of tests must be completed to determine if the satellite is 
operating properly. 
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7.2.1.10 Sustained Ops 

Sustained operations will begin when the satellite has passed all tests. This 
is the scientific data collection phase. Operation time lines, that outline step by 
step procedures, must be written for periods of time, such as a month. 


7.2.2 Detailed Test Time Line 

The entire satellite must be tested. That includes tests for all components, 
down to the last wire. Thus it is necessary to develop a time line that will orderly 
layout the testing periods. This will make it possible to provide certain tests with 
information collected during a previous test. The following list of tests has a 
format consisting of three parts, the begin and completion dates, equipment used 
to conduct the test, and a brief description. 


7.2.2. 1 Structures Group 

7.2.2. 1.1 Buckling Analysis of Gravity-Gradient Boom 

Begin: Jan 93 

Completion: May 93 

Equipment: ANSYS, Workstation 

Determine the required boom sizing for the boom to withstand space environment. 


7.2.2. 1.2. Boom Thermal Bending 


Begin: 

Completion: 

Equipment: 


Jan 93 
May 93 

ANSYS, Workstation 
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Determine the bending properties of the boom and if necessary change the 
boom's design characteristics to operate safely with respect to thermal fluctuations. 


7.2.2. 1.3 Outgassing Properties of Composite Materials 

Begin: Jan 93 

Completion: May 93 

Equipment: Vacuum chamber, scale, MTS, etc... 

This experiment exposes composite material to simulated space 
environment. Changes in weight and material properties should be recorded. 


7.2.2. 1.4 Material Property Tests 

Begin: Jan 93 

Completion: May 93 

Equipment: MTS and computer/software package 

Sufficient data should be taken so that accurate material properties can be 

used in the finite element modeling. 


7.2.2. 1.5 Finite Element Analysis 

Begin: Jan 93 

Completion: May 93 

Equipment: ANSYS, Workstation 

This test will verify that the boom will withstand the forces it will 

experience during flight. This test should also simulate the launch phase. 
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7.2.2. 1.6 Release Mechanism Verification 


Begin: 

Completion: 

Equipment: 


June 93 
Dec 93 

Command equipment 


Verify consistent operation. 


7.2.2. 1.7 Dynamic Testing 

Begin: Sep 93 

Completion: Dec 93 

Equipment: Shaker, Accelerometers, computer, etc... 

Vibrational testing should simulate both launch loads and transportation (to 
launch site) loads. This experiment should be the combined effort of both the 
structures and integration teams. 


7.2.2. 1.8 Non-Destructive Evaluation 

Begin: Sep 93 

Completion. May 94 

Equipment: Dependent on type of material 

Verify that composites and/or aluminum is free of defects. 


7.2.2.1.9 Strain Gauge Check 


Begin: 

Completion: 

Equipment: 


Sep 93 
Dec 93 

Ohm meter, solder, wire, shaker 
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This experiment will determine the status of the strain gauge and how they 
will react to vibrational loads. 


7.2.2.1.10 Integrated Structure Tests 

Begin: Jan 94 

Completion: Aug 94 

Equipment: Shaker, accelerometers, computer, etc 

(foam core components) 

The structure will be tested to withstand launch and transportation loads. 
The structure will include the strain gauges and components represented by foam 
core materials (this will determine if all the components will fit). The test will also 
determine the data collection of the strain gauges. 


12.2.2 Power Group 

7.2.2.2.1 Solar Array Output 

Begin: Sep 93 

Completion: Dec 93 

Equipment: Volt & Amp meters. Simulated sun source 

This experiment will determine the status of each individual solar cell. 


1. 2.2.22 Battery Capacity 

Begin: Sep 93 

Completion: Dec 93 

Equipment: Volt & Amp meters, load 


Verify the battery parameters. 
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1.2.2.2.3 Power Regulation Device 

Begin: Sep 93 

Completion: Dec 93 

Equipment: Volt & Amp meters, power supplier, power storage 

device, and load 

This test will determine the operational status of the power regulation 

device. 

7.2.2.2.4 Internal Meter Check 

Begin: Nov 93 

Completion: Dec 93 

Equipment: Volt & Amp meters, loads, power supply and 

storage devices 

This test will verify the subsystem, that determines the status of the entire 
power system, is working 

7.2.2.2.5 Integrated Power System Tests 

Begin: Jan 94 

Completion: May 94 

Equipment: Volt & Amp meters, load, simulated Sun 

Determine the characteristics of the entire power system. Verify the 
charging of the batteries directly through the solar array and so forth. 
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7.2.2.3 Attitude Control and Determination Subsystem 

7.2.2.3.1 Gravity-Gradient Simulations 

Begin: Jan 93 

Completion: June 93 

Equipment: Workstation/Computer, Code 

The code developed for this simulation should be general, allowing for the 
user to adjust the satellite's properties (distribution of mass). 


7.2.2.3.3 Attitude Sensors 

Begin: Sept 93 

Completion: Oct 93 

Equipment: ? 

Determine if the sensors operate within manufactures' parameters. 


7.2.2.3.4 Torque Rod Test 

Begin: Oct 93 

Completion: Dec 93 

Equipment: ? 


7.2.2.3.5 Integrated ADC Test 

Jan 94 
May 94 

Mockup satellite, 3 dimensional rotating stand, 
simulated forces. 


Begin: 

Completion: 

Equipment: 
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This test will place a mockup version of ISAT-1 in a frame that is allowed 
to rotate in 3 dimensional space. The mockup will have simulated forces applied to 
it in order to determine the status of the torque rods (This test may not be valid, 
since the test is on the Earth's surface). 


7.2.2.4 Payload 

7.2.2.4.1 Accelerometer Data 

Begin: Jan 93 

Completion: May 93 

Equipment: Mock-up, accelerometers, computer, code, bb gun. 

Using computer simulation, determine if the combination of accelerometer 
data can be triangulated to determine the locations of impacts. Develop required 
software. 
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7.2.2.4.2 Laser Reflectors 

Begin: Jan 93 

Completion: Aug 93 

Equipment: ? 

Determine the feasibility of the laser ranging experiment. 


7.2.2.4.3 Plant Experiment 

Begin: Jan 93 

Completion: Aug 93 

Equipment: ? 

Determine the feasibility of the plant experiment. 


7.2.2.4.4 Robotic Arm 

Begin: Jan 93 

Completion: Aug 93 

Equipment: Computer, Simulation Code 

In conjunction with attitude group, determine the stability of the entire 
satellite when the arm is functioning. Also, determine, via computer simulation, 
the envelope of the arm. 


7.2.2.4.5 Radiation Detector 


Begin: 

Completion: 

Equipment: 


Sep 93 
Dec 93 

Radiation source. Reference detector, computer. 
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7.2.2.4.6 Acoustic Sensor/Micrometeoroid Detector 

Begin: 

Sep 93 

Completion: 

Dec 93 

Equipment: 

? 

This is a follow up on the proof of concept test concerning the 

accelerometers. 


7.2.2.4.7 CCD Camera 


Begin: 

Sep 93 

Completion: 

May 94 

Equipment: 

Thermal Chamber, shaker, ? 

Determine the survivability of the CCD Camera. 

7.2.2.4.8 Telescope 


Begin: 

Sep 93 

Completion: 

May 94 

Equipment: 

? 

Determine if the telescope can switch power settings. 

7.2.2.4.9 Robotic Arm Mock-Up 

Begin: 

Sep 93 

Completion: 

Dec 93 

Equipment: 

ISAT-1 mock-up and command software 


Chapter 7 Ground Systems and Operations 


86 


Using command software, can the robotic arm be controlled? Can failures 
be corrected. 

7.2.2.4.10 Video Display 

Begin: Nov 93 

Completion: May 94 

Equipment: Thermal chamber, Volt & Amp meters, video 

camera 

Determine the characteristics of the robotic arm camera. 

7.2.2.4.11 Status Check Subsystem Tests 

Begin: June 94 

Completion: Aug 94 

Equipment: ? 

Tests must be done to determine the operational status of each experiment 
package. Certain failure scenarios should be listed and then determine how each 
scenario can be concluded upon, via sensors. 

7.2.2.4.12 Meteoroid Location System Tests 

Begin: Jan 94 

Completion: May 94 

Equipment: ISAT-1 mockup, accelerometers, computer, code, 

micrometeoroid equivalents. Volt and amp meters. 

Check if the system still functions with the predicted 
momentum transfers. Measure power consumption. 
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7.2.2.4.13 Telescope and CCD Camera Combo 

Begin: Apr 94 

Completion: Aug 94 

Equipment: Telescope, CCD camera, visual reference, computer. 

Determine the compatibility of the two units and the operational reliability 
using commands via a computer. 

7.2.2.5 Thermal Control System 

7.2.2.5.1 Thermal Couple Check 

Begin: Sep 93 

Completion: May 94 

Equipment: Volt meters, thermal chamber 

12 . 2 . 5.2 Thermal Determination Subsystem 

Begin: Sep 93 

Completion: May 94 

Equipment: ? 

Determine the operational status of the thermal detective devices in order 
to determine the thermal environment. 

7.2.2.5.3 Integrated Thermal Control Test 


Begin: 

Completion: 

Equipment: 


June 94 
Aug 94 

Thermal chamber, IR source, etc... 
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7.2.2.6 Communications 

This topic can be found in the Summer 1992 report. This section was not 
revised due to insufficient knowledge of the subject. 


7.3 Data Management 

Data management is concerned with how much data in bits of information 
each subsystem will be allowed to generate. To determine this, it was necessary to 
discuss the overall configuration of the satellite and ground system, and the orbital 
characteristics. 


7.3.1 Satellite & Ground System Configuration 

The satellite configuration has three levels. The payload and health check 
systems are represented by the first level. On the second level are the two flight 
computers and data storage devices. On the final level are the two communication 
packages. These comm units transmit to the ground antenna at Boone which is 
connected to a networking system. This network consists of the primary 
computer, tracking computer, data storage computer and filtering computer. 

The data created by the subsystems, health checks, and monitoring systems 
will be routed to one of the flight computers via a serial junction. The data will 
either be directly sent to the next serial junction or stored in the memory of the 
flight computer to be sent later. The second serial junction will route the data to a 
comm unit for transmission back to Boone. There is an inherent limiting factor in 
these communication packages. They operate at 9600 baud (9600 bits per 
second). 
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7.3.2 Orbital Characteristics 

Since the rate of the transmitter is known, the next step is to determine the 
amount of time we have to transmit. This is inherent in the satellite's orbital 
characteristics. 

A selection of orbits have been selected for analysis. Each orbit has two 
sets of data. First is the line of sight data, when a person at Boone can observe the 
satellite. The second set of data is the communications envelope. Any time the 
satellite is within this envelope, a communications link can be established. This 
data gives the time of pass, pass duration, time since last pass, and total data 
transmitted during the duration. For example, taking an orbit of 800 km in altitude 
and inclination of 50 degrees a pattern can be seen. Generally the passes come in 
sets of 5 and 4 with a hour and three quarters time interval, see Figure A. 7.1 and 
Table A.7.1. Each pass ranges up to 6 minutes in length. Between each set of 
passes there is roughly a 17 hour period. The main item to note, is that the orbit's 
durations are not constant throughout the set of passes. Each pass must be dealt 
with individually. 


7.3.3 Data Rate 

With this overall picture it was possible to determine when and how long 
the satellite is in communications range. The data displays the time for each pass 
over an arbitrary six day week. With the duration of each pass and the rate of the 
transmitter known, it is possible to create a basic data transmitting schedule. 

An example of data transmitting is given below, 


9600 baud 


9600 bits per second 

1 Byte 

= 

8 bits 

9600 baud 

= 

1 200 bytes per second 


C-2- 
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For 8 minutes of transmitting (which corresponds to the average maximum 
duration). 


576,000 bytes transmitted 
(.576 Megabytes) 

This corresponds to about a little over a third of the memory on a 3.5 inch 
High Density floppy (IBM formatted 1.44MB). Thus, it will be necessary to limit 
the data flow of each subsystem. 


7.3.4 Data Rationing 

Having only a limited amount of data that can be transfer it is necessary to 
delegated sizes of data creation to individual subsystems. Each orbit sets an Xi 
amount of data that can be transmitted during each pass i. This amount can be 
calculated from the Orbital Appendix (duration of the pass times the data rate). 
For example, the 45th pass for the 800 km altitude, 60 degree inclined orbit has a 
duration of 3.77 minutes. This corresponds to a total data transmission of 4.35 
Mbits (having two communication units operating at 7600 baud). Each subsystem 
is given a portion of this amount to transmit down to the Boone station. Each 
group must utilize this transmitted data in an useful manner. 


7.3.4. 1 Data Receiving Groups 

The data will be divided into five groups. Command and Control, 
Structures, Attitude Determination and Control, and Power will be portioned a 
minimal (enough to be useful) amount of data. This amount of data will be a set 
value; every pass will consist of this size of data transmission. The Payload group 
will have the left over amount. This amount must be divided among each payload 
subsystem (decided by the payload group and specific event). 



Chapter 7 Ground Systems and Operations 


91 


7.3.4.2 Command and Control 

The data collected for Command and Control is basically confirmation 
requests. During operations, certain events are initiated. If an event is initiated 
prematurely, by accident, etc it may be detrimental to the satellite. Thus, each 
command will be followed by a transmitted confirmation statement. For example, 
an operator types in the commands to extend the boom. The satellite receives this 
uplink and automatically transmits a confirmation statement where the operator 
will reply either yes or no (whatever the conditions are at that time). 


7.3.4.3 Structures 

Multiple strain gauges and other devices will be attached to the structure in 
various locations (locations that are important). These measuring devices will 
create data that will be processed into structural status report. 


7.3.4.4 Attitude Determination and Control 

ADC will be given a certain amount of data. This amount must be 
manipulated in a fashion so that the attitude of the satellite can be determined in 
the least amount of space. This will require analysis in what type of measurements 
are to be made and how often they are made. 


7.3.4.5 Power 

The power group will be delegated an amount of data transmission to 
determine the operational status of the system. Voltage and current checks must 
be done in various areas of the subsystem to determine if each component is 
working properly (solar cells, batteries, regulation device, etc ). 
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7.3.4.6 Payload 

The payload group will receive any left over amount of data that the other 
subsystems do not use. This data must be delegated by the group among the 
various payload packages. 


7.3.4.7 Data numbers 

As stated earlier. Command and Control, Structure, ADC, and Power will 
be given a set amount of data transfer per pass. Each group will receive 10,000 
bits of data per pass. Depending on the duration of the pass, the payload will pick 
up the rest of the data transfer. 


7.4 Personnel 


It will be necessary to hire personnel to operate the satellite. The length of 
time for each set of passes is roughly 7 hours for 5 pass set (set type A) or 5.3 
hours for 4 pass set (set type B). If we allow an hour and a half prior to the initial 
pass and an hour of logging out (including composing a report concerning the pass 
set), operators will work a 10 or 7 hour shift. One factor to note, is the shifts do 
not always come at desirable time (such as a 7 to 5 job). Some shifts may start a 
midnight while others may start in late afternoon. Whatever the case may be, the 
Boone station must be manned 17 hours out of a 48 hour period. 
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7.5 Ground Systems Selection 
7.5.1 Purpose 

The final area focussed on is ground systems selection. This is the 
selection of hardware for the Boone ground station. The hardware must operate 
properly; communicate with and track the satellite, store data, monitor health 
status, and initiate events. 


7.5.2 Ground Station Layout 

The main room contains the control three-bay, educational desks, 
equipment rack and secretary's area. The storage room will be used for storage 
and as a snack room for operators, it will also contain the large screen TV. The 
executive room will consist of the meeting table and the engineers work space. As 
of March 15, 1993 the equipment that will be purchased is shown in Table 7.1. A 
schematic diagram of the communications system is shown in Figure A. 7.2. 


Ferrups UPS RE1 .8KVA 
Citadel Power Conditioner CLT1001RTW 


Laseijet 4 
Deskjet 550C 
DraftPro DXL 
NP2020 Copier 
Scanjet lie 
Laseijet FAX 


80486/50 Control Computer 
20" SVGA Color Monitor 
80386SX/20 Control Computer 
20" SVGA Color Display 
80386/33 Data Archive Computer 
17" SVGA Color Display 
80486/50 Engineer's Computer 
20" SVGA Color Display 
Macintosh 2vx Education Computer 
16" SXVGA Color Display 


Table 7.1 Ground Equipment 
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The furniture will be from Rockwell and AMCO. The equipment rack, 
three-bay and some of the furniture will be purchased from AMCO. Rockwell will 
donate some furniture and possible some hardware. 


7.5.3 Network 

The computers at the Boone Ground Station will be consolidated in a 
networking scheme. This will provide faster processing of data by having the 
individual computers parallel process their designated data. For example, the data 
for tracking the satellite will be directed to the tracking computer where it will be 
processed, while at the same time, the health status data will be sent to its 
computer for processing. Several options have been looked at: UNIX, Novell, 
Windows NT. At the time of this report, the process of choosing has just started. 
The Windows NT, though, looks the most promising. It is however not on the 
market. It is believed that it is in the beta version phase and that a copy will be 
made for ISAT testing usage. 


7.6 Future Work 

Only a portion of Ground systems and Operations has been covered to 
date. The following list is certain areas that need to be determined, 

1 . ) A final detailed test time line 

2. ) A revision of the data management 

— > Need a final orbit 
3 .) Final ground systems selection 

4. ) Event Scenarios 

— > Problem solving procedures 

5. ) Personnel Management 


Study of the above areas will be attempted, time permitting. 
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Chapter 8 Communication System 
8.1 Summary/Overview 

Functionally the communication system of the satellite can be divided into 
two categories. The main comm system which is a part of the satellite bus and the 
experimental comm system which is a part of the satellite payload. For the purpose 
of identification the first one is designated as COMM-1, Figure 8.1, while the other 
as COMM-2, Figure 8.2. 

Practically the two systems are not totally independent. They share some of 
the common system elements and also one of them (COMM-2) acts as a standby 
to the other (COMM-1). 

In the later part of this section, design considerations for COMM-1 and 
COMM-2 are presented. 



Figure 8.1 COMM-1 Block Diagram 



Figure 8.2 COMM-2 Block Diagram 
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8.2 COMM-1 
8.2.1 Purpose 

This would be the main communication link between satellite and ground 
station. All the Command and Control data would be handeled through this link. 
Additionally, following other types of data would also be transferred using this 
link: 


Uplink: 

a. Picture display data. 

b. Onboard computer software update data. 

Downlink: 

a. Bus health status. 

b. On-board experiments data. 

c. Earth imaging data. 

d. Ground sensors data. 

8.2.2 Data Rate / Bandwidth Calculations 

Average number of passes per day 4 

Average contact time per pass 10 min. or 600 sec 

Uplink: 

Data to be uplinked: 

a. Command and Control data 

b. Software Update 

c. Picture display data 

Assumption of data rate less than 10k bps 

Bandwidth requirement 10k Hz 


Downlink: 

Data to be downlinked: 

a. Command and Control data 

b. Bus status report 

c. Earth imaging data 

d. Ground sensors data 


Assumption of data rates: 
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Ground sensors 10k bps 

All other 1 Ok bps 

Total 20kbps 

Bandwidth requirement 20k Hz 


8.2.3 F requency Selection 

The recent World Administrative Radio Conference (WARC) in Spain 
allocated new frequencies for small satellite systems. These are: 

Frequency Band Purpose 

137.000— 137.025 MHz Primary 

137.175— 137.825 MHz 

148.000— 149.900 MHz 
400.150— 401.000 MHz 

137.025— 137.175 MHz Secondary 

137.825— 138.000 MHz 

Although the ultimate selection of frequency is subject to the availability of 
a particular frequency band and approval by the FCC, following are the relative 
merits and demerits of above mentioned frequency bands: 

a. At higher frequencies the available bandwith is larger 

b. At higher frequencies antenna sizes become smaller, thus easier to manage. 

c. Free space loss is higher at higher frequencies. 

d. Doppler shift is more significant at higher frequencies than lower ones. 

Considering all these factors and our system requirement it appears that the 
choice of higher frequency band would be more suitable for our satellite. 


8.2.4 Link Calculations 

Following link equation is used for link calculations: 

(SNR) = (EIRP) + (Gr/Ts) - (FSL) + (228.6) 
where 

(SNR) = signal to noise ratio 

(EIRP) = equivalent isotropic radiated power 

(Gr) = receiver antenna gain 

(Ts) = equivalent receiver noise temperature 

(FSL) = free space loss 
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Link calculations are done at two different frequencies. 

fl = 150 MHz wavelength 1 = wll = 2.0 m 
f2 = 450 MHz wavelength2 = wl2 = .67 m 

Minimun antenna angle above horizen = 15 degrees 

Orbital altitude 550 km 

Maximum distance between satellite and ground station 
R = 550/Sin 15 or 2125 km 

(FSL)l = 142 db 
(FSL)2= 152 db 

Downlink calculations: 


(Ts) 

= 1000 K 

(Gr) 

= 10 

(Gr/Ts) 

= -20 db 

(EIRP) 

= Pt*Gt 

Pt 

= 2 Watts 

Gt 

= 1 

(EIRP) 

= 3 db 


(SNR)l = 3 - 20 - 142 + 228.6 = 70 db 
(SNR)2 = 3 - 20 - 152 + 228.6 = 60 db 

Maximum data rate = (MDR) = (SNR) - (Eb/No) 
where 

(Eb/No) = Min ratio of signal power per bit to noise power density 

= 15 db (typical value corresponding to Bit Error Rate of 10E-6 

(MDR)l = 70 - 15 = 55 db = 300k bps 
(MDR)2 = 60 - 15 = 45 db = 30k bps 


Uplink calculations: 


Pt 

= 2 Watts 

Gt 

= 10 db 

(EIRP) 

= 13 db 

Gr 

= 0 db 

Ts 

= 1000 K 

(Gr/Ts) = 

-30 db 
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(FSL)l = 142 db 
(FSL)2 = 152 db 

(SNR)l = 13 - 30 - 142 + 228.6 = 70 db 
(SNR)2 = 13 - 30 - 152 + 228.6 = 60 db 


8.2.5 Maximum Doppler Shift 

(to be determined) 


8.2.6 Contact Time Calculations 

Radius of earth 
Height of satellite 
Value of u 
Orbital Period 

Ground antenna min elevation angle 
Coverage Sector angle at ground 
Coverage Sector angle at earths center 
Coverage Time 


6378.165 km 
550 km 

3.986032 10E5 
95.65 min 
15 degrees 
75+75 degrees 
37 deg 
10 min 


8.2.7 System Reliability/Redundancy Requirements 

Since this is our main communication link with the satellite a high level of 
accuracy and system reliability is required. Typically a bit error of one in 1000000 
is considered reasonable. To achieve high degree of reliability following 


8.2.8 System Weight/Size 

(to be determined) 

8.2.9 Antenna Design 

Ground Tx/Rx Antenna : 

A directional tracking antenna would be employed for Ground Station in 
order to achieve maximum contact time per pass and reasonably high directivity. 

Satellite Tx/Rx Antenna: 


(to be determined) 
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8.2.10 Power Consumption 

On-board Transmitter 6 watts (only during transmission) 
On-board Receiver 2 watts (operational mode) 

1 watt (standby mode) 


8.2.11 System Cost 

(to be determined) 


8.2.12 Manufacturer Information 

1 . Cynetics Corporation 

P.O. 2422 
3824 Jet Drive 
Rapid City , SD 57709 
Phone (605)394-6430 
Fax (605) 394-6456 

2. Motorola Inc. 

Strategic Electronics Division 
Advanced Programs Development 
2501 S. Price Road 
Chandler , A Z 85248-2899 
Phone (602)732-3015 


8.3 COMM-2 
8.3.1 Purpose 

This experimental communication system aims at establishing a data 
link with the satellite from a small inexpensive and portable ground 
transponder. The idea is to have multiple ground transponders (from 100 to 
200) spread all over the state which could relay some kind of local data 
(like temperature, pressure, humidity, wind velocity etc.) up to the satellite 
whenever the satellite passes over them . 


8.3.2 Data Rate / Bandwidth 

Average contact time per pass 
Average number of passes per day 
Data uplink per transponder per pass 
200 transponders 


8 min or 480 sec. 
4 

24k bits 
4.8M bps 




Chapter 8 Communication System 


103 


System data rate 


4.8M/480 or 10k bps 


System bandwidth 


12kHz 


8.3.3 Frequency Selection 

The recent World Administrative Radio Conference (WARC) in 
Spain allocated new frequencies for small satellite systems . These are: 

Frequency Band Purpose 

137.000— 137.025 MHz Primary 

137.175— 137.825 MHz 

148.000— 149.900 MHz 
400.150— 401.000 MHz 

137.025 — 137.175 MHz Secondary 

137.825— 138.000 MHz 

Although the ultimate selection of frequency is subject to the 
availability of a particular frequency band and approval by the FCC, 
following are the relative merits and demerits of above mentioned 
frequency bands: 

a. At higher frequencies the available bandwidth is larger compared to 

lower frequencies . 

b. At higher frequencies antenna sizes become smaller , thus easier to 

manage . 

c. Free space loss is higher at higher frequencies . 

d. Doppler shift is more significant at higher frequencies than lower ones 

Considering all these factors and our system requirement it appears 
that the choice of higher frequency band would be more suitable for our 
satellite . 

8.3.4 Multiple Access Technique 

Various possible multiple access techniques were studied and the 
relative merits and demerits were evaluated. 


a. Frequency Division Multiple Access: 
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1. Each of the 200 Ground Transponders will have a distinct carrier 
frequency. Data will be uplinked at a rate less than 100bps using BPSK. 
Inter-channel spacing will be extremely narrow, ie. of the order of few 
hundred hertz. 

2. The receiver on-board satellite (Rx-2) will receive all the 200 
channels, downconvert them to an IF frequency and then separate them 
out. Data of each channel can be stored in the memory or directly sent to 
the downlink transmitter (Tx-2). Channel separation for so close channels 
could be very complex. 

3. Each Ground Transponder will have a small very simple receiver 
that would detect the presence of signal coming from the TT&C 
transmitter (Tx-1) on-board satellite. 

4. When the detected signal crosses certain threshold, it would be 
an indication that the satellite is within the range of the Ground 
Transponder, so it will start transmitting the stored data and will continue 
transmission until the satellite gets out of range. 

5. It is expected that each Ground Transponder will barely manage 
to uplink all the data that was required to be transferred. 


b. FDMA(Bent-pipe scheme): 


1. Each of the 200 Ground Transponders will have a distinct carrier 
frequency. Data will be uplinked at a rate less than 100bps using BPSK. 
Inter-channel spacing will be extremely narrow, ie. of the order of few 
hundred hertz. 

2. The receiver on-board satellite will receive all the 200 channels 
together as a single band and will convert the whole block to an IF level. 
The IF signal will be fed to the downlink transmitter Tx-2 for re- 
transmission to the Ground-Station. 

3. At the Ground-Station the receiver will split the individual 
channels apart and detect the data in each channel. 


4. Each Ground Transponder will have a small very simple receiver 
that would detect the presence of signal coming from the TT&C 
transmitter (Tx-1) on-board satellite 
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5. When the detected signal crosses certain threshold, it would be 
an indication that the satellite is within the range of the Ground 
Transponder, so it will start transmitting the stored data and will continue 
transmission until the satellite gets out of range. 


c. Packet Data Transfer Scheme: 


1. All the 200 Ground Transponder will have identical transmitter 
that share a common frequency band of around 50Khz. 

2. Each Ground-Transponder will have a receiver to detect the 
presence of satellite and to decode various messages from the satellite. 

3. When a Ground-Transponder detects the presence of satellite, it 
transmits its ID and a "data transfer request" message. 

4. The satellite receiver Rx-2, receives and decodes the ID and the 
message and transmits a "request acknowledge" message along with the ID 
of the specific Ground-Transponder, using another transmitter Tx-3. 

5. The Ground-Transponder receives the message and starts 
transmitting its data packet. 

6. When the satellite receiver Rx-2, receives the end-of-data header, 
it transmits a "receipt acknowledge" message. 

7. The "request acknowledge" signal from the satellite serves as an 
indication to the other Ground-Transponders that the channel is busy while 
the "receipt acknowledge" message indicates that the channel has become 
idle. So another Ground-Transponder can "ask" for data transfer. 

8. The only possibility of Ground-Transponder data clash is during 
the initial data transfer request phase. If that happens, the satellite receiver 
won't be able to decode the message so the satellite will transmit a "data 
clash" message. 

9. In such a situation, the Ground-Transponders will retry the "data 
transfer request" after some random delay to avoid a clash again. 

10. The transmission of various messages from the satellite to the 
Ground-Transponders will be quite frequent, ie. for each Ground- 
Transponder to transfer data, at least two messages needs to be transmitted 
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by the satellite. This amounts to an average of about 400 messages per pass 
of about 8 minutes. 

Comparison of Various Schemes: 

1 . FDMA Scheme: 

Advantages: 

1. No additional transmitter on-board satellite is required. 

2. Ground-Transponder receiver can be very simple. 
Disadvantages: 

1 . Satellite receiver design for the separation of closely 
spaced channels can be very complex. 

2. Ground-Transponder transmitter design would be tough, 
since the demand on frequency stability would be less than 
100 hz at around 200 Mhz. 

3. Doppler frequency shift can be a serious problem. 

2. Bent-Pipe Scheme: 

Advantages: 

1 . No additional transmitter on-board satellite is required. 

2. Ground-Transponder receiver can be very simple. 

3. No on-board data storage is required. 

4. Satellite receiver would be very simple. 

Disadvantages: 

1 . Ground-Transponder Transmitter design would be tough, 
since the demand on frequency stability would be less than 
100 hz at around 200 Mhz. 

2. Doppler frequency shift can be a serious problem. 

3. Ground-Transponder needs to transmit higher power to 
maintain the same bit-error-rate. 
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3. Packet Data Transfer Scheme: 

Advantages: 

1 . Simple satellite receiver design. 

2. Much simpler Ground-Transponder transmitter design. 

3. Doppler frequency shift won't be a serious problem. 
Disadvantages: 

1 . Additional transmitter on-board satellite (Tx-3) is required. 

2. Ground-Transponder needs a relatively complex 
receiver to decode various messages from the satellite. 


8.3.5 Link Calculations 

Following link equation is used for link calculations: 


(SNR) = 

(EIRP) + (Gr/Ts) - (FSL) + (228.6) 
where 

(SNR) 

= signal to noise ratio 

(EIRP) • 

= equivalent isotropic radiated power 

(Gr) 

= receiver antenna gain 

(Ts) 

= equivalent receiver noise temperature 

(FSL) 

= free space loss 


Link calculations are done at two different frequencies . 
fl = 150 MHz wavelength 1 = wll = 2.0 m 
f2 = 450 MHz wavelength2 = wl2 = .67 m 

Minimun antenna angle above horizen = 30 degrees 

Orbital altitude = 550 km 

Maximum distance between satellite and ground station 
R = 550/Sin 30 or 1100 km 

(FSL) = (4*22*R/7*wl)**2 


(FSL)l = 137 db 
(FSL)2 = 146 db 
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Uplink calculations: 


Pt 

2 Watts 

Gt 

10 db 

(EIRP) 

13 db 

Gr 

0 db 

Ts 

1000 K 

(Gr/Ts) 

-30 db 

(FSL)l 

137 db 

(FSL)2 

146 db 


Maximum data rate = (MDR) = (SNR) - (Eb/No) 
where 

(Eb/No) = Min ratio of signal power per bit to noise power density 

= 1 5 db (typical value corresponding to Bit Error Rate of 1 OE-6 

(MDR)l = 70 - 15 = 55 db = 300k bps 
(MDR)2 = 60 - 15 = 45 db = 30k bps 

(SNR)l = 13 - 30 - 142 + 228.6 = 70 db 
(SNR)2 = 13 - 30 - 152 + 228 6 = 60 db 


8.3.6 Maximum Doppler Shift 

(to be determined) 


8.3.7 Contact Time Calculations 


Radius of earth 
Height of satellite 
Value of u 
Orbital Period 

Ground antenna min elevation angle 
Coverage Sector angle at ground 
Coverage Sector angle at earths center 
Coverage Time 


6378.165 km 
550 km 

3.986032 10E5 
95.65 min 
30 deg degrees 
60+60 deg 
30 deg 
8 min 
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8.3.8 System Reliability / Redundancy 

Since this would be an experimental system its reliability constraints 
both in terms of bit-error-rate and equipment redundancy would be much 
lesser than COMM- 1 . A BER of 1 0E-5 is considered suitable for the type 
of data to be normally handeled by this system . 


8.3.9 System Weight / Size 

(to be determined) 


8.3.10 Antenna Design 

(For Satellite Receiver) 

It is assumed that the frequency bands for all the four links of 
COMM-1 and COMM-2 (ie. up and down links of COMM-1 and up and 
down links of COMM-2) are adjacent to each other. Therefore a common 
antenna can be used for both COMM-1 and COMM-2. A circulator would 
be required to connect the onboard Transmitter(s) and Receivers) to a 
common antenna. The uplink signals of COMM-1 and COMM-2 would be 
separated by means of suitable filters, while the downlink signals of 
COMM- 1 and COMM-2 can be combined through a matchingh unit . 


8.3.11 Power Consumption 

On-board Receiver 

operational mode 4 watts 

standby mode 2 watts 


8.3.12 System Cost 

(to be determined) 

8.3.13 Perspective Suppliers 

a. Cynetics Corporation 
P.O. Box 2422 
3824 Rapid City , SD 37709 
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Figure A. 1.4 Batteries 
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Figure A. 1.6 Flight Communications 
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Figure A. 1.7 Camera Module 
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Figure A. 1.1 2 Packing Arrangement 
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Table A.2.1 Mass Budget 
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ISAT Structural Maee Estimate 


| ISAT Structural Mass Estimate 






















note : 

- same structural conf iguration for Al & Gr/Ep | 





- fasteners not included 







- bus sensors/wiring/etc. 

not included 





- Aluminum density taken from €061 -T6 

: 2710 kg/m~3 

... 

r 


- Graphite Epoxy density: 

1600 kg/m ^3 




























[SIDE PANEL: 

each panel 

is 0.170 by 

0.643 






area : 

0. 10931 










' 



ItOP/ BOTTOM PANEL: 

area of hexagon: ( 1/2 ) *apoth.em*perimeter 





apothem: 

0 . 17 







perimeter 

1 . 02 







area 

0.0867 













STRINGER: 

a = outer flange length 






t = thickness of stringer 






area: t*a 

- 0. 5 (t**2) (tan 30) 












SHELVES : 

each shelf 

is a hexagon 

L 






area : 

0.0867 


r 



































SIDE PANELS: 
















Panel 

Total 


Mass 

Mass 

Panel 

^ ... 


Thickness 

Area 

Volume 

Al 

Gr/Ep 

Thickness 



(m) 

(nT2) 

<m~3) 

(kg) 

(kg) 

( in) 



















0.001016 

0.65586 

0. 000666 

1 .80582 

1 . 06617 

0.040 



0. 001366 

0.65586 

0. 000896 

2 . 42790 

1.43345 

0 . 054 



0. 001716 

0.65586 

0. 001125 

3 . 04999 

1.80073 

0.068 



0.002066 

0. 65586 

0.001355 

3 . 67207 

2. 16801 

0.081 



0. 002416 

0.65586 

0.001585 

4 . 29415 

2. 53529 

0 . 095 



0.002766 

0.65586 

0.001814 

4 . 91623 

2.90257 

0.109 



0. 003116 

0.65586 

0 . 002044 

5 . 53832 

3.26986 

0.123 


_ 

0. 003466 

0.65586 

0.002273 

6 . 16040 

3 . 63714 

0.136 







. - j 



TOP/ BOTTOM PANEL: 




. 









__ j 



Panel 

Total 


Mass 

Mass 

Panel 



Thickness 

Area 

Volume 

Al 

Gr/Ep 

Thickness i 



(m) 

(m-2) 

(itT3) 

(kg) 

(kg) 

(in) 



















0 . 001016 

0.1734 

0.00018 

0.47743 

0.28188 

0.040 



0.001366 

0.1734 

0.00024 

0.64190 

0.37898 

0.054 





0.001716 

0.1734 

0.00030 

0.80637 

0.47609 

0.068 



0.002066 

0.1734 

0.00036 

0.97084 

0.57319 

0.081 



0.002416 

0 . 1734 

0 . 00042 

1.13531 

0.67030 

0.095 



0.002766 

0.1734 

0.00048 

1 . 29970 

0.76740 

0.109 



0.003116 

0. 1734 

0.00054 

1.46425 

0.86450 

0.123 



0.003466 

0 . 1734 

0.00060 

1.62872 

0.96161 

0.136 



0.003816 

0 . 1734 

0.00066 

1.79319 

1.05871 

0.150 





0.004166 

0.1734 

0.00072 

1.95766 

1.15582 

0.164 



0.004516 

0.1734 

0.00078 

2.12213 

1.25292 

0.178 



0.004866 

0.1734 

0.00084 

2.28660 

1.35002 

0.192 



0.005216 

0.1734 

0.00090 

2.45107 

1.44713 

0.205 



0.005566 

0.1734 

0 . 00097 

2.61554 

1.54423 

0.219 



0.005916 

0 . 1734 

0.00103 

2.78001 

1.64134 

0.233 



0. 006266 

0.1734 

0.00109 

2.94448 

1.73844 

0.247 



t 







SHELVES: 
















Shelf 

Total 


Mass 

Mass 

Shelf 



Thickness 

Area 

Volume 

A1 

Gr/Ep 

Thickness 



(m) 

(m^) 

{m“3 ) 

(kg) 

(kg) 

(in) 




0.002032 

0.2601 

0. 00053 

1.43230 

0.84564 

0.080 

0.002382 

0.2601 

0. 00062 

1 . 67900 

0 . 99129 

0 . 094 

0.002732 

0.2601 

0.00071 

1.92571 

1.13695 

0.108 

0.003082 

0.2601 

0.00080 

2 . 17241 

1.28261 

0.121 

0 . 003432 

0.2601 

0 . 00089 

2 .41912 

1.42826 

0.135 

0.003782 

0.2601 

0 . 00098 

2.66582 

1.57392 

0.149 

0. 004132 

0 . 2601 

0.00107 

2.91253 

1 . 71957 

0.163 

0.004482 

0.2601 

0 . 00117 

3.15923 

1.86523 

0.176 

0.004832 

0.2601 

0.00126 

3.40594 

2 . 01089 

0.190 

0.005182 

0.2601 

0.00135 

3 . 65264 

2 . 15654 

0.204 

0.005532 

0.2601 

0.00144 

3.89935 

2.30220 

0.218 

0.005882 

0.2601 

0.00153 

4.14605 

2.44785 

0.232 

0 . 006232 

0 . 2601 

0 . 00162 

4 . 39276 

2.59351 

0.245 

0. 006582 

0.2601 

0.00171 

4 . 63946 

2 . 73917 

0.259 

0 . 006932 

0.2601 

0.00180 

4 . 88617 

2.88482 

0.273 

0.007282 

0.2601 

0.00189 

5.13287 

3 . 03048 

0.287 

0 . 007632 

0.2601 

0. 00199 

5.37958 

3 . 17613 

0.300 

0.007982 

0.2601 

0. 00208 

5.62628 

3.32179 

0.314 

0.008332 

0.2601 

0. 00217 

5.87299 

3.46745 

0.328 

0.008682 

0.2601 

0.00226 

6.11969 

3.61310 

0.342 

0.009032 

0.2601 

0.00235 

6.36639 

3.75876 

0.356 

0.009382 

0.2601 

0.00244 

6.61310 

3.90441 

0.369 

0.009732 

0.2601 

0.00253 

6 . 85980 

4 . 05007 

0.383 


STRINGERS : 


Flange 

Length 


Flange 

Thickness 


Mae s Flange 

Gr/Ep Thickness 




(m) 

(m) 

(m-2) 

(m~3) 

(kg) 

(kg) 

(in) 









10% 

0.0170 

0. 001016 

1.697E-05 

1 . 09E-05 

0.01746 

0.01746 

0. 040 


0.0170 

0. 001366 

2 . 268E- 05 

1.46E-05 

0.02334 

0.02334 

0. 054 


0.0170 

0.001716 

2.832E-05 

1.82E-05 

0.02914 

0.02914 

0.068 


0.0170 

0.002066 

3.389E-05 

2 . 18E-05 

0.03487 

0.03487 

0.081 


0.0170 

0.002416 

3 .939E-05 

2 . 53E-05 

0 . 04052 

0.04052 

0.095 


0. 0170 

0.002766 

4 .481E-05 

2 . 88E-05 

0. 04610 

0.04610 

0. 109 


0.0170 

0.003116 

5 . 017E-05 

3.23E-05 

0.05161 

0.05161 

0 . 123 









15% 

0.0255 

0. 001016 

2. 561E-05 

1.65E-05 

0.02635 

0.02635 

0.040 


0. 0255 

0. 001366 

— 

3 .429E-05 

2 . 21E- 05 

0. 03528 

0.03528 

0. 054 


0. 0255 

0.001716 

4 . 2 9 IE- 05 

2 . 76E- 05 

0.04414 

0.04414 

0.068 


0.0255 

0.002066 

5.145E-05 

3 . 3 IE- 05 

0.05293 

0.05293 

0.081 


0.0255 

0.002416 

5. 992E-05 

3 . 85E-05 

0.06165 

0.06165 

0.095 


0.0255 

0. 002766 

6. 832E-05 

4 . 39E- 05 

0.07029 

0.07029 

0. 109 


0.0255 

0.003116 

7. 666E-05 

4 . 93E-05 

0. 07886 

0.07886 

0.123 









20% 

0.0340 

0. 001016 

3 . 425E 05 

2 . 2E - 05 

0. 03523 

0.03523 

0.04 0 


0. 0340 

0.001366 

4 -591E-05 

2 . 95E - 05 

0. 04723 

0.04723 

0.054 


0. 0340 

0.001716 

5.749E-05 

3 . 7E- 05 

0. 05915 

0.05915 

0.068 


0. 0340 

0.002066 

6 . 901E- 05 

4.44E-05 

0.07100 

0.07100 

0.081 


0.0340 

0.002416 

8 . 046E-05 

5. 17E-05 

0.08278 

0.08278 

0.095 


0.0340 

0.002766 

9.184E-05 

5. 91E- 05 

0. 09448 

0.09448 

0.109 


0.0340 

0.003116 

0. 0001031 

6 . 63E-05 

0.10611 

0.10611 

0.123 


Table A.2.2 Structural Mass Estimate 




Appendix A.3 
Payload 






















Appendix A.7 

Ground Systems and Operations 




Figure A.7.1 Orbital Pass History 
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r ^ ^ 9 t ^ ^ 9 CO O N ^ N 1 ^ r- 1 - ^ t- |\ |\ f- ^ r- i- ^ r- N CO ^ O CO (D 
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CO 
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00 
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CO 
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CO 
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CO 
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1.78 

LL l 

1.76 

00 
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CD 
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I 
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1 
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CM 
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CM 
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00 
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00 

00 
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CO 

05 

05 
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00 

05 

05 

05 
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W r 1 1 ten :ov 

4 3 
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4. -? 

t na nr 
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in r _ t to r. ; 

4 4 


4 4 


4 4 

e i s u o n i : 
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: 
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assumes ai i st nge r s same 
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oj ■ 
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n u Hi d e r o f . i .a i . a is- o * ! 1 1 s 


N 0 T E . th.-c oaraiTicte^ listed be lew contain the arrsc- 

d_ mens i ora ncs parameters used t-hr oughc-u t the d “ooi 
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so . i o ) 


C 0 it M o N / 

INPUT 

NSTRNG , 5 

TRH 

OR , o \< i 

n 

COriMQN / 

LOADS 
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coordinates of 7 he section *- ot mas: 
coordinates of tne section * pt mar 


tens number of data points to prose 
the number of stringers 
tne n u m oe r of s p a r c a o s 


IhLIZATION 
d j. r / v ar c c! . f 
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This subroutine creates a unifies. data matrix of the 
stringers and skins 

VARIABLE LISTING 

R - real 

DP - double precision 
1 - integer 

STRNGR(DP) - data array for the stringer 
skin(DP) - data array for the wing skin 
ALLDAT ( DP )- the unified data array 
ALlDAT array structure 
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we i g h t / m 2 of stringer 
mas s of s t r i nge r ( spa r ) 
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"b" dimension of stringer 
skin thickness for panel 
weight/ in of skin 
G "of skin 

mass of skin added to strings 1 ' 
new value of point mass equivalent 


note: skin data is controlled by the first index of aildat 

e . g , alldat(i,S) is the thickness of panel 1 , NOT panel 

1(1)- iooo control variable 
NSTRNGfl i - the number of stringers 

NUMPT£( I i *■ the number of data points in the current arrav 
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call ar ea (a.b. t.ar ; 
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CONTINUE 


ADD SKIN DATA TO ALL STRINGER VALUES 

note: for version i all skin panels are the same - 

a method needs to be determined for getting the proper 
skin properties with mid-side stringers when that happens 
use the first panel ent^y for now 

do 5l 0 i = i . ns t r ng 

ALLDAT ( I ,8 ) -skin( 1 , 1 5 
ALLDATt I , 9 )=skin( 1,2) 

ALLDAT ( I , 10 ) = skint 1,3) 

SET ADDED SKIN MASS AND POINT MASS VALUE TO ZERO 
ALLDAT ( I . 11 ) = 0 . 0D0 
ALLDAT ( I , 1 2 ) - 0 , 0D0 
con ti nue 


Sort the stringer and spar caps into clock-wise fashion, 
beginning with the front upper* stringer 

note. this sort worked once - hopefully it still works 1 DGC 4/3/93) 

Split, upper and lower stringers for sorting 
NUPR. ‘0 
N LWR = 0 
NC0L = 12 


numpts- ns t. r ng+ neaps 

DO 300 1=1, NUMPTS 
Z- ALLDATt 1,2) 

IT ( z . GT . 0 . ODO ) TH EN 
NUPR = NUPR+ 1 
DO 310 J=1 , NCOL 

UPR ( NUPR . J ) = ALLDAT ( I „ J ) 
CONTINUE 
else 

N L WR. - N L WR+ 1 
DO S20 J = 1 .. NCOL 

LWR \ N LWR . J ) = ALLDAT ( I , J ) 
CONTINUE 
END IF 
CONTINUE 

SORT UFR INTO ASENDING ORDER 
CALL SORTUPl UPR. NUPR, 12,1 .SRT1 ) 
SORT LWR INTO DESCENDING ORDER 
CALL SORTDN « LWR . N LWR ,12. 1 ,SRT2 ) 


p f C‘ C CT ivf D I jT 



/ /f Iyer/ hd/ug rad/36 1/cooper/aer o464 .dir 


Page 


840 

830 


860 

850 

* 


870 

* 


995 

990 


+ 

* 


DO 830 1=1 ,NUPR 

DO 840 J=1 .NCOL 

ALLDAT ( I, J)=SRT1(I. J) 

CONTINUE 

CONTINUE 

LIMIT = NUPR+NLWR 

DO 850 I=NUPR+1 , LIMIT 
J = I-NUPR 
DO 860 K= 1 , NCOL 

ALLDAT ( I , K ) =SRT2( J , K ) 

CONTINUE 

CONTINUE 

ADD 1ST ELEMENT TO END TO COMPLETE “CIRCLE" 

DO 870 J = 1 , NCOL 

ALLDAT (LIMIT+i.J) = ALLDAT ( 1 . J ) 

CONTINUE 

CHECK THE DATA AT THIS POINT - OUTPUT TO FILE ALLDAT 1 . DAT 

OPEN ( UN IT = 3, FILE= ’ . ,/ou tpu t. dir/all da tl . dat ’ , status: ’ unknown ’ ) 

write! 3 , * ) ’ This is the data at the end of UN1FER’ 
write(3,* ) ' it should have all original data sorted’ 
write! 3, * ) ’ in clockwise order’ 

do 995 i=l,numpts 

wri te( 3 , 990 )(alldat(i,j),j=l,l2) 
continue 

format! lx , 12dl6 . 5 ) 
close(3) 


retu r n 
END 


NOTE: SPAR CAP COORDINATES ARE THE INPUT COORDINATES NOT 

THE CENTROIDS OF THE SPAR CAPS 



ANSYS Finite Element Model 


/inter, no 
/prep7 

/title, ISAT Structure - Model 2 
et, 1, 63 

r, 1, 0 .001016, 0.001016,0.001016,0.001016 

mp, ex, 1, 100e6 

mp, dens, 1,2. 6265e-4 

et,2, 63 

r, 2, 0.002032, 0.002032, 0.002032, 0.002032 

np,ex,2,100e6 

mp, dens, 2,2 .6265e-4 


n, 1, 0 .170, 0, 0 
n, 2, 0.15725, -0.0220836,0 
n, 4, 0.09775, -0.1251164,0 
fill 

n, 5, 0.085, -0.1472 
n, 6, 0.0595,-0.1472 
n, 8, -0.0595, -0.1472 
fill 

n, 9,-0.085,-0.1472 
n, 10, -0.09775, -0.1257764 
n, 12, -0.15725, -0.0220836 
fill 

n, 13,-0.170,0 
n, 14,-0.15725,0.0220836 
n, 16,-0.09775,0.1257764 
fill 

n, 17, -0.085, 0.1472 
n,18, -0.0595, 0.1472 
n, 20, 0.0595, 0.1472 
fill 

n, 21, 0.085, 0.1472 
n, 22, 0.09775, 0.1251164 
n, 24, 0.15725, 0.0220836 
fill 

n, 25, 0, 0, 0 

n, 26, 0 . 085, 0 

n, 27, 0.0425, -0.0736 

n, 28, -0.0425, -0.0736 

n, 29, -0.085,0 

n, 30, -0.0425, 0.0736 

n,31, 0.0425, 0.0736 


ngen, 3, 31, 1,3 1,1, 0,0, 0.063 9 
ngen, 4,31,63,93,1,0,0,0.0762 
ngen, 4, 31, 156,186,1,0,0,0.063866667 
ngen, 2, 31, 249, 279, 1,0, 0,0. 095 

type, 2 
real, 2 

en, 1,1,32,33,2 
egen, 9 , 31 ,1,1,1 



en, 10,4,35,36,5 
en, 11,5,36, 37, 6 
egen, 9,31,10,11,1 

en, 30, 8, 39, 40, 9 
en,31, 9,40,41, 10 
egen, 9,31,30,31, 1 

en, 50, 12, 43, 44, 13 
en, 51, 13,44, 45, 14 
egen, 9,31,50,51,1 

en, 70, 16, 47, 48, 17 
en, 71, 17, 48, 49, 18 
egen, 9,31,70,71, 1 

en, 90, 20, 51, 52, 21 
en, 91,21, 52, 53, 22 
egen, 9,31, 90, 91,1 

en, 110, 24, 55, 32,1 
egen, 9,31,110,110,1 

type, 1 
real, 1 

en, 120, 2, 33, 34, 3 
en, 121, 3, 34, 35, 4 
egen, 9,31, 120 , 121,1 

en, 140, 6, 37, 38, 7 
en, 141, 7, 38, 39, 8 
egen, 9, 31, 140, 141,1 

en, 160, 10,41,42, 11 
en, 161, 11, 42, 43, 12 
egen, 9,31,160,161,1 

en, 180, 14, 45, 46, 15 
en, 181, 15, 46, 47, 16 
egen, 9, 31, 180, 181,1 

en, 200, 18, 49, 50,19 
en, 201, 19, 50, 51, 20 
egen, 9, 3 1,2 00, 2 01, 1 

en, 220, 22, 53, 54, 23 
en, 221, 23, 54, 55, 24 
egen, 9,31,220,221, 1 

type, 2 
real, 2 

en, 240, 31, 25, 26, 26 
en, 241, 26, 25, 27, 27 
en, 242, 27, 25, 28, 28 
en, 243, 28, 25, 29, 29 
en, 244, 29, 25, 30, 30 
en, 245, 30, 25, 31, 31 
en, 246, 22, 31, 26, 24 
en, 247, 24, 26, 1,1 



en, 248, 1,26, 2, 2 
en, 249, 2, 26, 27, 4, 4 
en, 250,4,27,5,5 
en, 251, 5,27,6,6 
en, 252, 6, 27, 28, 8 
en, 253, 8,28, 9, 9 
en, 254, 9, 28, 10, 10 
en, 255, 10, 28, 29, 12 
en, 256, 12, 29, 13, 13 
en, 257, 13, 29, 14, 14 
en, 258, 14, 29, 30, 16 
en, 259, 16, 30, 17, 17 
en, 260, 17, 30, 18, 18 
en, 261, 18, 30, 31, 20 
en, 262, 20, 31, 21, 21 
en, 263, 21, 31, 22, 22 

egen, 2,62, 240 , 263,1 
egen, 2, 93, 264, 287 , 1 
egen, 2, 93,288, 311, 1 
egen, 2, 31, 312, 335,1 

d, 80, all, 0 
d, 82 , all, 0 
d, 173, all, 0 
d, 175, all, 0 

f , 254, fz, -9 
f , 275 , f z, -9 
f , 276, fz, -15 
f , 256 , f z , -7 . 5 

f ,258,fz,-7.5 
f ,277, fz, -10.5 
f ,260,fz,-9 

f , 264, fz, -1 .5 
f ,278, fz, -10.5 
f , 265, f z, -3 
f , 266, f z, -9 

f , 279, fz, -15 
f , 268, f z, -7 . 5 
f ,270,fz,-7.5 
f ,274,fz,-10.5 
f , 249, fz, -3 
f , 250 , f z , -1 . 5 

f ,252,fz,-1.5 
f , 253, fz, -3 

f , 173, fz, -2 .375 
f , 175, fz, -2 .375 
f,185,fz,-13.375 
f,186,fz,-13.375 

f ,184,fz,-16 
f , 183, fz, -11.2 
f , 182, fz, -11.2 
f , 181, fz, -16 



f , 180, fz, -38 .2 


f / 68, f z, - 
f , 70, fz, - 
f , 72, fz, - 
f , 74, f z, - 
f , 76, fz, - 
f,78,fz,~ 
f,79,fz,~ 
f , 80, f z, - 
f , 82 , fz, - 
f , 83, fz, - 
f , 84, fz, - 
f , 86, f z, - 
f ,64,fz,- 
f , 66, fz, - 

f,89,fz,- 
f , 90, f z, - 
f , 91, fz, - 
f,92,fz,~ 
f , 93 , f z, - 
f , 88, f z, - 
f , 87, f z, - 


iter, 1, 1, 
af write 
finish 


/eof 


3 

3 

1.1 

1.1 

1.5 

2.5 
2 

3.75 

3.75 
2 

2.5 

1.5 

2.75 

2.75 

9.076 

5.76 
4.2 
9.576 
9.576 
7.756 
6 


1 



Appendix B.6 
Attitude and Control 



Dynamic Equations 

For a rigid spacecraft the dynamic equations of motion are given by Euler’s equation 

—-H = -Jx/? + fj + T e 

at 

H = Iu 


where 


H = 

angular momentum of the spacecraft 

I = 

moment of inertia tensor of the spacecraft 

w = 

spacecraft angular velocity 

Td = 

disturbance torques 

T c = 

control torques 


In the principal axis coordinate system the dynamic equations in component form are 

= T d s d” Tc t + (/jr - It ) 

= U, + T Cp + (It ~ I,)u M W t 
= U, + r c , + (/, - I f )u x w t 


du x 

u ~dT 

du v 

lv ~dt 

du, 

Il ~dT 


Kinematic Equations 


The kinematic equations using Euier symmetric parameters can be expressed as follows 


where 


and 


dq 

dt 



n = 


“ 0 

W, 

— w y 

Wc 


0 

Wf 

Wy 

** 

-Wg 0 

w. 

m — w* 


f — 

0 

w s £3 

+ x 

' o ’ 

-—Wo 

0 



where A is the attitude matrix and w 0 is the orbital rate. 


The angular rates are related as follows: 


and conversely 



cos B 

sin 0 tan ^ 
— sin $ sec ^ 




Wy 

= 


. 


cos 6 
0 

sin 6 


0 

1 

0 


0 sin# 1 T w, *1 

1 — cos 0 tan ^ I I w y I 

0 cos ^ sec J [ w, J 

— sin 0 cos <f> 1 T 4 1 

sin^ 0 

cos 0 cos $ J l ij> J 
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Attitude Matrix 


The attitude of the spacecraft can be described as a transformation from the nominal reference frame to 
the body frame. Here, the reference frame is centered at the spacecraft’s center-of-mass and is composed of 
the y-axis parallel to the negative orbit normal, the z-axis toward nadir, and the x-axis normal to both the 
y-axis and z-axis. For circular orbits, the x-axis is parallel to the spacecraft’s velocity vector. The relation 
between the nominal reference frame, E , and the body reference frame, e } is 


Sr 


' E t ‘ 

ty 

. Cl . 

= A 

E, 

E, 


where A is the attitude rotation matrix. Using the yaw-roll-pitch (3-1-2) Euler rotation sequence the attitude 
matrix, A, can be expressed as 


A = 


cos 0 cos ip — sin 0 sin <j> sin ^ cos 0 sin ip + sin 0 sin <f> cos %!> — sin $ cos ^ 

— cos <f> sin cos <j> cos xp sin <f> 

sin B cos rp + cos 0 sin ^ sin ^ sin 0 sin — cos 6 sin ^ cos ^ cos $ cos <f> 


where ^ = roll, 0 = pitch, and = yaw. 

The Euler angles can be obtained from the elements of the matrix: 

<t> = arcsin ( 023 ) - 90° < <f> < 90° 

6 = arctan 0 < 0 < 360° 

V> = arctan 


f-— ) 

v <»«/ 


0 < i> < 360 # 


Alternatively, the attitude matrix can be expressed using Euler symmetric parameters ql, q2, q3, q4 as 
follows 

fi “ 92 “ 9i + 9? 2 (9192 + 9394 ) 2 (9193 — 9294) 

2 (9192 ~ 9394) -9? + 92 - 93 + 94 2 (9293 + 9194) 

2(9193 + 9294) 2(fl293 -9l94) -9? -92 + 93 + 94 J 


A = 


and conversely, the Euler symmetric parameters can be found from the attitude matrix 

94 

92 


^( 023 - 032 ) 


4^ (a 31-°13) 


93 = 


^(012-021) 
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Gravity-gradient Torque 

The gravity-gradient torque T lt about the center of mass is given by 

T »=% [*• x ('■"•)] 

where 

/i = Earth’s gravitational constant 

R» = geocentric radius vector to the spacecraft’s center-of-mass in the body reference frame 
R$ = magnitude of R, 

R$ = unit vector of R t 
I = moment-of-inertia tensor 


Solar Radiation Torque 


The solar radiation torque on the spacecraft, r to i ar is the vector sum of all the torques on the individual 
surfaces generated by a solar radiation force on those surfaces, and is given by 


n 

T, oiar = ^2 Ri X F{ 
1=1 


where Ri is the vector from the spacecraft center-of-mass to the center-of-pressure of the sth element. 
The force generated by solar radiation on a fiat plane is given by 


where 


P = -PAaosB [(1 - C,)5+ 2 (c, cos* + ivj 


P 

A 

S 

N 

0 

C t 

Cd 


mean momentum flux acting on a surface normal to the solar radiation 
surface area 

unit vector from the spacecraft to the sun 
unit outward normal vector of the surface 
angle between S and N 
coefficient of specular reflection 
coefficient of diffuse reflection 


Aerodynamic Torque 

The aerodynamic torque, f atro , is calculated the same way as the torque due to solar radiation. The 
aerodynamic torque is calculated by summing the vectors of the torques on individual surfaces generated by 
the force due to aerodynamic drag, i.e. 

w 

T.'ro = y^ Rj X Fj 

•sl 

where again Ri is the vector from the spacecraft center-of-mass to the center-of-pressure of the tth element. 
The aerodynamic force on a flat plane with area, A , is given by 

P = -\c d pV 2 A (ti v)v 

and the force on a sphere with radius, R> is 

P = -\CdpV 2 *&V 

m 


1 



where 


Cd = coefficient of drag 

p = atmospheric density 

V = velocity vector of spacecraft relative to the atmosphere 

V = unit vector of V 

V = magnitude of V* 

N = unit outward normal vector of plane 

Magnetic Torque 

The torque due to the spacecraft’s magnetic moment is given by 

f m = rh x B 

where 

T m = magnetic torque acting on the spacecraft (N • m 2 ) 
m = spacecraft magnetic moment (A m 2 ) 

B = Earth’s magnetic flux density ( Wb/m 2 ) 


Earth’s Magnetic Field 

The Earth’s magnetic field, J3, can be represented as the gradient of a scalar potential function, V 

B = -VV 

and V can be represented by a spherical harmonic model 

V (r, 6, 4) = a^2 (- )" +1 £ (s"- m cos + A"- m sin nufi) P" m (0) 

nxsl msO 

where 


r = 

geocentric distance 

9 

coelevation 

* 

east longitude from Greenwich 

a = 

equatorial radius of the Earth (6371.2 km) 

0 and A = 

coefficients of the IGRF model 

m = 

order of model 

n — 

degree of model 

P = 

Gaussian normalised Legendre function 


The field in topocentric coordinates is then 

. - /a\ ^ 

B r = ]T(-) + 


n=l 


msO 


b, = 


•4^ /0\"+ J aP"- m 

" E (;) E (*"’"* coe m* + /»"-"* sin m*) W 

n=l msO 


B 4 = ~ E (r)" +J E n»(— y B,m sin m^ + h n,m cos m^)P n,m (0) 

nsl m=0 
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where 

B r = radial component , outward positive 
B$ = coelevation component, South positive 
B+ = azimuthal component, East positive 
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Various Coordinate System Transformations 

The inertial cartesian coordinate system is related to the spherical local tangent coordinate system by 



where R is the orthogonal rotation matrix 

sin <j> cos 0 cos cos 0 — sin 0 

R = sin ^ sin 0 cos <f> sin 6 cos 0 

cos <f> — sin 4> 0 

where 

^ = coelevation 

0 = right ascension 


The polar coordinate system is related to the perifocal coordinate system by 



where R is the orthogonal rotation matrix 

cos v —sin v 0 * 

R = sin v cos v 0 

0 0 1 _ 

where v is the true anomaly 

The inertial cartesian coordinate system is related to the perifocal system by 



where R is the rotation matrix 



cos ft cos u; — sin Osina; cost 

— cos SI sin u> — sin ft cos w cos* 

sin 0 sin « 


R = sin (l cos w + cos ft sin w cost 

— sin ft sin w + cos ft cos w cost 

— cos ftsint 


sin us sin t 

cos u) sin t 

cost 

where 




Q = 

longitude of the ascending node 



U) = 

arguement of periapsis 




t = inclination 
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DE SIGN, FABRICATION, AND TESTING OF A 
GRAPHITE-EPOXY COMPOSITE GRAVITY-GRADIENT BOOM 
FOR A SMALL SATELLITE 
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Dr. T. J. McDaniel 3 and Dr. Leveme K. Seversilce 4 , Advisors 
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Ames, Iowa 50011 


Abstract 

Passive stabilization methods for satellites have 
undergone extensive research and development 

Recently the number of small satellites (satellites less 

than 100 kg.) has increased dramatically. This has 

lead to increased use of passive stabilization methods, 

such as gravity-gradient The core of a 
gravity-gradient stabilization system is a deployable 
boom with a damping mechanism. Traditionally, this 
boom is constructed from metal alloys. Uneven 
heating and cooling occurs when these alloys are 
exposed to varying solar radiation conditions. This 
can induce thermal vibrations which can lead to 
undesired satellite attitude inversions. Graphite-epoxy 
composites can be fabricated to minimize thermal 
expansion. This will be beneficial when applied to 
gravity-gradient booms. The goal of this project is to 
demonstrate the use of graphite -epoxy composites in 
gravity-gradient booms. This project encompasses: 
the use of a sate llit e attitude simulation program for 
boom sizing and determination of gravity-gradient 
boom loading, development of joint-locking 
mechanisms for boom deployment, and selection and 
testing of appropriate fabrication methods. 

Gravity-Gradient Stability 

The Earth’s gravitational field provides a ready source 

of stability for satellites in low Earth orbit. The 
inverse square nature of the Earth’s gravitational fi»Vl 

causes a torque to act on a spacecraft in orbit This 

torque causes the spacecraft’s principal axis of 
minimum moment of inertia to align itself with the 
local vertical. 1 


The magnitude of the torque is a function of the 
distance from the center of the earth (orbital radius) 
and die ratios of the mass moments of inertia. A 
cartesian coordinate system is introduced with the 
origin at the spacecraft center of mass, with the x-axis 
along the orbital velocity vector, y-axis perp endicular 
to the orbital plane, and the z-axis nadir p ointing and 
completing a right-handed system. 



Figure 1. Spacecraft Axis System 2 
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Figure l shows a model satellite in Earth ortrit far 
reference. The spacecraft’s stability in a gravity-field 
can be measured by the introduction of three stability 
criteria 1 : 



(la) 


(lb) 

©.-a,-!,)/!. 

(lc) 


These criteria represent the stability of a satellit e as 
the ratios of mass moments vary. This is shown 
graphically in Figure 2. Most spacecraft are 
symmetric about the z-axis so that 1,=^ When 1,*^ 
the region of stability lies along the 6, axis, and the 
only criterion of interest becomes 6,. With an 
axisymmetric spacecraft, equation la above reduces 
to: 


The determining factor for stability then becomes the 
ratio of 1, to t Note that as I, increases, ©, 
approaches its limiting value of unity. Equation 2 
suggests that shapes with an 1JI* ratio near zero are 
the most stable. This corresponds geometrically to 
long, slender objects. 

Spacecraft design considerations such as 
arrangements, thermal control, and available launch 
envelopes do not always allow for the use of long 
slender bodies. However, the mam moments of 
inertia of spacecraft may be modified once in orbit to 
achieve the desired stability. While many possibilities 
exist to achieve a favorable mass moment ratio, QJIJ 
one of the most common methods is the deployment 
of a long boom with a tip mass; more often known as 
a gravity-gradient boom. 


e.«l-(I,/U (2) 



Figure 2. Stability Held 
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Gravity-Gradient Booms 

Considerable research effort has been directed toward 
development and testing of gravity-gradient booms. 
These efforts include the development of several types 
of booms. A bi-metallic boom was most frequently 
used. This consists of two thin strips of n»*nHir 
•Boy approximately 0.005 centimeters. (0.002 in.) 
thick. Beryllium-aluminum and beryllium -copper 
woe commonly used in this application. Thes e two 
•hips were carried into orbit in coils, much as 
common tape measures are stored. A motorized unit 
would unwind the two tapes. The tapes would then 
buckle together, inter-locking to form a closed 
section. This type of boom continues to be used 
today. Figure 3 shows some representative of this 
type of boom. 



Figure 3. Tape boom examples 2 

Several difficulties arose with this design. Despite 
the closed section, the deployed boom is wed: in 
torsion. The design has also been observed to 
experience thermal flutter. Thermal flutter is the 
bending of the boom out of the desired plane due to 
uneven heating. This bending changes the 
spacecraft’s inertia properties and causes undesirable 
attitude behavior. The behavior has been severe 
enough to cause satellite a t titude inversions as 


documented in the HILAT and Polar BEAR 4 - 3 
spacecraft. 

As more small sa t ell i t es are constructed, complexity 
becomes a planning consideration. The support 
system for the tape-design booms includes a 
motorized deployment system. This system 

introduces m ec h a n ical uncertainty and power 
requirements into the design of a spacecraft 
Simplification of this system would reduce spacecraft 
complexity and remove potential sources of failure. 

Composite Mate rial Application 

The two most significant drawbacks of the currently 
available gravity-gradient booms are their complexity 
and undesirable thermal behavior. Only limited 
success has been achieved with the n* of coatings to 
control the thermal properties of metal alloys. The 
mechanisms for deployment of metal alloy booms 
also remain complex. The application of composite 

material technology can be used to address these two 

considerations. 

Composite materials can be fabricated with favorable 
thermal expansion characteristics. By varying ply 
orientations and sequences, it is possible to create a 
material with a near zero coefficient of thermal 
expansion. Composite materials designed with low 
thermal expansion coefficients could be applied to 
current design configurations. This should reduce the 
effects of thermal flutter and its aMrafaifrvi satellit e 
a ttitu d e disturbances. 

The mechanical properties of composite materials also 

allow for several possible design change , The 
favorable strength to weight ratio of composite 
materials allows for smaller closed section designs. 
This facilitates the use of a telescoping design which 
can be carried internally as a closed section. 
Telescoping closed section designs do not require a 
drive motor for deployment. A closed section may be 
deployed by a spring system or by use of a gas 
charge to pressurize the internal volume of the boom. 
Either of these deployment systems would require less 
volume and would be less complex than a motorized 
unit 


3 



Boom Section Development 

Woric began to study various boom section 
configurations to determine their feasibility. At this 
point* no consideration was given to final boom size 
(length, diameter, tip mass, etc.) or loading 
conditi on s. Rather, efforts were focused on the cross- 
sectional properties and methods of deployment. 

Figure 4 shows the first open-section concept. This 
de si gn would consist of two "C" sections with flanges 
attached. To store this structure, the smaller iwtins 
"C" section would have to buckle into the larger "C 
section. 



Two test specimens of this design were fabricated. 
The first test specimen consisted of three plys of 
uniaxial graphite-epoxy composite [0°^0°,0°]. The 
specimen was oven -cured under vacuum. Two "C" 
sections were fabricated in succession, with the first, 
smaller section serving as a mold for the second, 
larger section. These two portions were then epoxied 
together to form a closed section. This specimen was 
evaluated qualitatively to determine its buckling 
ability. It was found that the smaller section would 
not buckle without loss of matrix integrity. 

The second specimen consisted of two plys of [0,90] 
woven graphite-epoxy composite. This specimen was 
constructed using the same method as the first and 
was also qualitatively evaluated. It was found that 


the woven material was still too "stiff" to buckle 
without significant cracks developing in the composite 
matrix. 

Figure 5 shows the second open-section concept 
This design would consist of concentric, interlocking 
"C" sections. The inner, smaller radius section, 
would slide in rails provided by die outer 
This design would provide for a telescoping boom. 



Figures. Open-Section Concept 2 

Two test specimens of this design were fahnctH 
The first, or inner section, was fabricated as a "C" 
section and cured in an oven. The second, or outer 
section, which provides the slide rails, was then 
formed over the inner section and cured. After the 
second cure, it was discovered that the two sections 
had fused. The fabrication process used did not 
a deq u a t ely separate the two sections to prevent epoxy 
from joining the two sections during the second cure. 
It was not possible lo evaluate this specimen. The 
second specimen was constructed in the same manner , 
but with added release material to prevent the fusing 
of the inner and outer "C" sections. This specimen 
also used a third outer section. This specimen was 
again qualitatively evaluated. It proved difficult lo 
fabricate a section with uniform guide rails. This 
created points along the structure that would bind 
during deployment or retraction. This difficulty arises 
out of the fabrication technique and the mandrel Hat 
was used. It is possible that with a different mandre l 
this problem could be corrected. In wtHirinn to 
deployment and retraction problems, this section 
proved very weak in torsion because of the open- 
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section design. The open-section also exhibits 
■•otwaufimn bending stiffhess, as the test specimen 
was found to bend very easily in the open direction. 

Rgore 6 shows the closed-section concept This 
design consists of concentric circular sections. This 
concept would be used much like a te l escoping radio 
antenna with the coocentric sections being slid along 
one another. 



Figure 6. Closed-Section Concept 


Deplovment/Retraction Methods 

Deployment schemes were considered in the 
development erf the three concepts outlined above. 
The double "C" section (Figure 4) was considered to 
be deployable and retractable by buckling of the 
t ra iler radius "C" section into the larger section, and 
then "rolling' these two sections on a spool (much 
like a retractable tape measure). The interlocking "C 
section (Figure 5) would be deployed and retracted by 
telescoping, as would the closed-section. 

Several mechanisms were considered to perform these 
tasks. A servomechanism would be required to 
deploy and retract the double ”C section by «nn»Hng 
a spool. This would add weight and complexity to 
the system. A boom of this type would be locked 
into place by the un-buckling of the maiw section . 

The interlocking "C section could be deployed by a 
spring. Retraction of this design would be more 


difficult One possibility was to attach a cord to the 
outer section and then reel this cord in to retract the 
boom. It is still uncertain what type of inrHng 

mechanism could be employed to lock the boom in a 

deployed configuration. 

The closed-section could be deployed by a 
charge. This was the simplest mechanism considered. 
The use of a cord and reel system could provide 
retraction capability. It may also be possible to 
design the magnetic portion of the satellite control 
system to avoid the need to retract the boom once it 
is deployed. 

Boom Section Selection 

The closed-section design was selected for further 
development This section was chosen far its 
torsional properties and its uniform bending stififoess. 
The simplicity of the deployment mechanism was also 
a consideration. The current assumption is that the 
m a gn etic portion of die *atrii«t* mtimd ** cont rol 
system will be able to invert the satellite with the 
boom deployed, so the boom will only be deployed 
once. 

Figure 7 and Figure 8 show the derails 0 f the joint 
design. Figure 7 is a cut-away of the outer wall of a 
section. Figure 8 is a view of an inner section, with 
the outer section not shown. This design provides for 
a forward stop collar located on the inner wall of the 
outer boom section, and a aeries of forward 
segmented stops located on the outer wall of the inner 
boom section. A segmented stop collar is also 
located on the inner wall of the outer section. A rear 
stop collar is also located on the inner section. 

When pressurized, the boom extends until the forward 

tiop collar on the outer section comes into «*iiyt 
with the segmented collar on the inner section. The 
segments on the inner section slide through the gqa 
in the collar on the outer section, and are then rotated 
to lock the boom in its deployed state. The joint 
design also provides for the placement of three 
P das: one each at the base of the inner section, the 

top and base of the rear stop collar to provide a gas 
seal for the pressurized deployment of the boom. 
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Figure 7. Joint Detail A 



Figure 8. Joint Detail B 


Boom Sizing 

The operational task of a gravity-gradient boom is to 
modify the inertia properties of a spacecraft to 
increase its stability. The degree of stability increase 
desired is the driving factor in sizing the boom’s 
length and tip mass to a particular satellite. Iowa 
State University is currently undertaking an effort to 
design, construct and launch a small satellite, ISAT-1 
(Iowa Satellite One). The preliminary design of this 



Figure 9. ISAT geometry 

satellite was used to size the boom length and imm 
for this project The current design configuration of 
ISAT is shown in Figure 9. 

Current efforts of the Iowa Satellite Project provided 
much of the needed information about the 
altitude dynamics. Mainly, the attitude deter mination 
and control group provided the mass moment inertias 
of die spacecraft body. These inertias assumed 
uniform mass distribution within the A 

design sizing code with boom length and tip mass as 
control variables was written to determine an 
optimum length and tip mass for the The 

code allows the user to select the desired stability in 
terms of the stability criterion ©,. The user also 
selects the desired range of boom lengths and dp 
masses. The code then iterates through these two 
variables, ca l culati n g the new mass moment of inertia 
I,, and uses the new inertia to calculate the value of 
© r The calculated value of 6, is subtracted from the 
target value, and the absolute value of the difference 
is written to a data file along with the boom length 
and mass. This information can then be plotted on a 
contour plot 
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Figure 10. Boom Sizing Contour Plot 


Figure 10 shows the contour plot for a stability 
criterion (8J of 0.98. Several candidate designs ran 
be seen to meet the desired stability. The final boom 
sue was s e lecte d to be 3 meters in length with a 6 
kilogram tip mass. This re pre s ents almost a fifteen 
fold increase in ^ 

Stability considerations determined the boom hmg th 
and tip mass. The diameter and thickness of each 
section would be determined according to the Wfe 
acting upon each section. acting on the 

gravid-gradient boom were modeled as the Earth's 
gravitational force acting on the tip mass. A section 
sizing code was written to size the boom sections 
based on the known loads, the material p rop erti es of 
graphite-epoxy composites 6 , and a selected factor of 
safety. It was not possible to accurately determine 


dynamic loads on the gravity-gradient boom due to 
insufficient source data from the BAT project, so the 
factor of safety used in the sizing analysis was 
increased. 

The code was set up to allow the user to choose the 
desired tip mass and boom length for section sizing 
analysis. The forces and moment on the boom were 
then determined based on the user's input of the ti p 
mass and the maximum deflection of the 
from the vertical axis. Once the user had com plet ed 
the inputs of number of sections, minim um section 
thickne ss, gap between sections, and initial b ase 
section outer diameter, the code would size the boom 
sections. Hus entailed determining the necessary 
outer diameters, section thicknesses, and inner 
diameters conesponding to each boom section. If the 
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lection thickness needed was found by the code to be 
smaller than the minimum input section thirimeas , the 
code would substitute in the minimum section 
t hickness and continue the analysis. This measure 
was i n c l ude d in order to allow for more reasonable 
composite layup thicknesses in the fabrication 
process. 

Results of section sizing yielded necessary section 
thi ck n esses for a number of base section outer 
diameters. A base section having an outer diameter 
of 3.81 centimeters was chosen based on the limited 
volume fraction of ISAT-1 which would be required. 
Table 1 shows the section sizes in terms of outer 
diameters, thicknesses, and corresponding inner 
diameters. It should be noted that in all cases the 
thicknesses used are the minimum section thicknes s, 
and the gap between sections is constant for ease of 
fabrication. 


Boom 

Outer 

Section 

Inner 

Section 

Dia. 

Thickness 

Dia. 

(#) 

(cm) 

(cm) 

(cm) 

1 

331 

0.102 

3.71 

2 

331 

0.102 

331 

3 

3.40 

0.102 

330 

4 

330 

0.102 

3.10 

5 

3.00 

0.102 

2.90 

6 

2.79 

0.102 

2.69 

Table 1. 

Boom Section Sizes 


Fabrication Methods 

Fabrication methods were developed to produce the 
necessary sections to conform to the design. This 
included the construction of suitable mandrels and 
development of layup procedures. This also included 
developing a fabrication technique for the section 
joints. 

The mandrels used for curing of the boom sections 
consisted of two sections of steel pipe, each 0.61 
meters in length. The two pipes had inside diameters 

of were 3.43 centimeters and 3.94 centimeters. These 
pipes were split along their length at the diameter. 
These were complemented by three rubber hose 
sections of 0.76 meters in length. The hoses had 
outside diameters 2.62 centimeters, 3.1S cenitmeters, 
and 3.61 meters. This provided the capability to 


fa b ri ca t e sections of varying diameter and 
The hoses were fitted internally with steel rods along 
their entire length. Both ends of the hoses were 
sealed, one co n ta in i ng a valve, to allow pressurization 
during the curing process. 

Each section was fabricated from three plys of woven 
[0° ,90°] graphite-epoxy composite. The composite 
material was cut to a length of 0-56 meters and a 
width equal to the circumference of the desired 
section. The three plys were then debullced with the 
edges staggered prior to being placed in the mandrel 
The debullced composite was then placed in the 
mandrel and wrapped around the hose. Finally, the 
pipe section was clamped around it Bleeder cloth 
was inserted between the composite material and the 
outer pipe as necessary to form each section. 

With the mandrel clamped together, the inner hose 
was pressurized to 275,7903 Pa (40 psi,), and the 
mandrel was placed in an oven to cure for 3 1/2 hours 
at 176° C. Upon completion of the cure, the mandrel 
was disassembled, and the completed boom section 
was removed. 

To avoid the added complexity and difficultly 
involved with co-curing the joints along with the 
respective boom sections, joints from separate cure 
cycles were epoxied in at a later time. A suitable gas 
seal was added to each section, and the boom was 
assembled. 

Testing 

Once the boom was complete, testing was c on duct ed 
to ensure that the mechanical performance of the 
gravity-gradient boom design was satisfactory. The 
first test to be run was an analytical test involving 
ANSYS finite-element analysis. The boom was 
modeled on ANSYS using plate elements. A static 
analysis was performed. Figure 11 gives a 
representation of the deflection resulting in a 
maximum load condition applied to the boom and its 
tip mass. The boom was found to have a maximum 
deflection at the tip of 0.05 meters (1.969 inches). 
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Figure 11. ANSYS Model 


Recommendations 

Further effort in this subject area is needed to enhance 
the understanding of the application of graphite -epoxy 
composites to gravity-gradient booms. The feasibility 
of the ap pl i cati on has been demonstrated. However 
additional development work and testing is needed to 
provide the validation necessary to prepare hardware 
for use in a space application. 

Additional testing of the boom section joint to verily 

die gas seal design and deployability is needed. The 

" Wgnjtude of gas charge needed for deployment in 
orbit is yet to be determined. Also, laboratory testing 
of the joint mechanics would serve to validate the 
length of the joint in respect to the boom sections. 

Additional mechanical testing would serve to 
characterize die composite material properties 
obtained under the cure conditions used in this 
project The repeatability of mechanical p mp w iir^ 
would build confidence in the production of flight 
hardware. 


This project has demonstrated the feasibility for use 
of graphite-epoxy composites in gravity-gradient 
booms for sm a ll satellites . Given the mechanical mi 
thermal properties of graphite-epoxy c omposites , a 
boom was designed and fabricated. The boom 
showed good attitude stability chararteri wi™ and 
encouraging mechanical behavior. 
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Introduction 

While in orbit, ISAT-1 will be in an environment, filled with meteoroids and orbital debris. Within 
2000 km of Earth's surface there is about 200 kg of meteoroid material, most of it only 0.1 mm in 
diameter. [1] In addition to this naturally occurring material there is estimated to be 1.5 to 3 million 
kg of man-made material as of mid-1988. Nearly 5000 objects with altitudes below 2000 km are 
continuously tracked by the U.S. Space Command Space Surveillance Center. Of these 5000 objects, 
only 6% are operational satellites and over 40% are fragments caused by breakups and collisions. [2] 
U.S Space Command can only track objects that are 10 cm in diameter or larger; therefore, any 
estimates of debris population smaller than this comes from impact and fragmentation modeling and 
from studying satellites that have been retrieved from orbit such as the Long Duration Exposure 
Facility (LDEF). There is an estimated 1000 kg of debris with diameters less than 1 cm and about 300 
kg of debris with diameters less than 1 mm. The greatest threat occurrs at altitudes lower than 1000 
km, or lower low Earth orbit (LEOl), the region where ISAT-1 will be operating. The population in 
LEOl is growing by about 120 objects per year as of 1989. The population is expected to increase 
substantially in the next decade with the international space station, proliferation of micro-satellites, 
and new collisions and breakups. Due to an increase in the use of higher inclination orbits over the 
last decade, the average relative velocity of debris has also increased substantially and is expected 
to continue. [3] 

Number of Impacts 

The number of impacts on a spacecraft depends on the particle flux and the spacecraft's size, 
shape, and orientation. Flux is the number of objects impacting on a randomly tumbling surface per 
unit time and area averaged over time. Because the particle flux is usually calculated for a tumbling 
surface, the k-faclor method is used to evaluate the expected number of impacts on a surface with 
a fixed orientation. The constant k is a function of the spacecraft's shape, as well as its zenith and 
azimuth angles. The zenith angle, 6. is the angle between the Earth vertical axis and the outward 
surface normal and the azimuth angle. <f > , is the angle between the spacecraft’s velocity vector and 
the outward surface normal. Values for the k-factor can be found from Figure 1. 

The number of impacts on a spacecraft, N , during a period of time is given by the following 
expression 


where 


N 



kiFrAidt 


n, — number of surfaces 

k = k-factor 

F t = cumulative flux on a randomly tumbling surface 

A = surface area 


( 1 ) 


Once N has been determined, the probability of exactly n impacts occuriug in the corresponding 
time interval is found from Poisson statistics 


P„ = 



( 2 ) 


1 
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Zenilh Angle 

THETA = 90 

- THETA = 75 AND 105 

— — THETA =60 AND 1 20 

— - - THETA = 45 AND 135 

THETA = 30 AND 150 

THETA = 15 AND 165 



Figure 1: K-factor For Values of Zenith and Azimuth Angles 



Meteoroids 


Meteoroids that stay with their natural parent body and create periods of high flux are called 
streams. Random fluxes are called sporadic. The meteoroid environment that was used in this 
analysis was comprised of the average sporadic meteoroids and a yearly average of the stream 
meteoroids. 

Flux 

The cumulative meteoroid flux, F r , for meteoroids of mass, m, is given by the following expression 

Fr = SjCbF^ p ( 3 ) 

where F* p is the interplanetary flux at 1 A.U. and is given by 

Fl v (m) = c& J(cim° 306 + Co) + c-j(m -f c^m 2 + csm 4 ) -0 36 + c 6 (»n+ c 7 m 2 )~° 85 j (4) 

where 


Co 

= 

3.156 x 10 7 

Cl 

= 

2.2 x 10 3 

Co 

= 

15 

C3 

= 

1.3 x l(T y 

Ca 

= 

10" 

C5 

= 

10 27 

Cf> 

= 

1.3 x 10~ 16 

cr 


10 6 


Sf is the Earth shielding factor and Ge is the focusing factor caused by the Earth’s graviational 
field. They can be calculated from the following expressions 



c _ (1 + cosrj) 
} 2 

(5) 


+ 

t-H 

II 

CD 

(6) 

and 




Re 



sin rj = — 

(Re + H) 

(7) 

where 



Re = 

Earth’s radius + 100 km atmosphere (0478 km) 


H = 

height above atmosphere 



t = orbit radius 
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Mass Density 


The estimated mass density for meteoroids, p, of mass, m, is given by the following expression 


p (gm/cm 3 ) 


Velocity Distribution 


2 m < 10“ 6 gm 
1 10“ 6 < m < 0.01 gm 

0.5 m > 0.01 gm 


( 8 ) 


Because of the precession of a satellite’s orbit and the tilt of the Earth $ equatorial plane with 
respect to the ecliptic plane, the meteoroid environment can be assumed to be omnidirectional 
relative to the Earth. The normalized velocity distribution with respect to the Earth is given by 


{ O.llJ 
3.32* 
1.69* 


112 11.1 < v < 16.3 km/s 

328 x 10 5 r” 5 34 1 6.3 < t; < 55 km/s (9) 

695 x 10~ 4 55 < v < 72.2 km/s 

The normalized velocity distribution is the fraction of impacts in a 1 km/s interval and is given 
by the following expression 


/'(«>) = 


/(*) 


( 10 ) 


FH*)dv 

The velocity relative to the spacecraft is just the vector sum of the meteoroid's velocity relative 
to the Earth and the spacecraft's velocity relative to the Earth. 


Orbital Debris 

Flux 

The cumulative flux, F rt of debris of diameter d and larger on a randomly tumbling surface is 
given by the following equations for altitudes below 2000 km. 


where 


F r = H<t>y {Figi + Fig?) 


H 

4>\ 




10 

1q( WB-lfe -1 5 ) 


, _ Ol 

* (*! + l) 

Fi = 1.22 x 10- s <r 2 5 

F 2 = 8.1 x 10 ,o (rf + 700) -6 

/ (1 + ,)"- 19881 for l < 201 1 

9 ' \ (1 + ?)- 3 (l + g')" -J011 ' for t > 2011 

92 = l + (p(‘~ 1988)] 


( 11 ) 

(12a) 

(12b) 

(12c) 

(12d) 

(12e) 

(12f) 

(12g) 
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where 


d 

t 

h 

S 

V 

q and q ' 


debris diameter, cm 
year 

altitude, km 

solar radio flux, F 10 . 7 . for t-1 year, 10 4 Jy 

assumed annual growth rate of intact objects in orbit, (p = 0.05) 
estimated growth rate of fragments in orbit, (q = 0.02, q ' = 0.04) 


and ¥ is the inclination function given by Table 1. 


Table 1: Inclination Function Values 


Inclination, deg 

¥ 

28.5 

0.91 

30.0 

0.92 

40.0 

0.96 

50.0 

1.02 

60.0 

1.09 

70.0 

1.26 

80.0 

1.71 

90.0 

1.37 

100.0 

1.78 

120.0 

1.18 


Mass Density 

Knowledge of debris shape and density is very scant. For the purposes of this analysis, the 
objects were assumed to be spherical. The average density, p, for an object with diameter, d, can 
be taken as 


P 



4 

2.8<r 074 


for d < 0.62 cm 
for d > 0.62 cm 


(13) 


Velocity Distribution 

The collision velocity distribution is the number of impacts with velocities between v and v + dv 
relative to the spacecraft. Averaged over all altitude's, the distribution is given by 


f(v) = {2vv 0 - v 2 } jcexp “ +Fex p ~ j + HC (4tn\, - v 7 ) (14) 
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where 


A = 
B = 

C = 

D = 
E = 

F = 


a = 

h = 

v„ = 


2.5 


f 0.5 

*<60 

< 0.5 - 0.01 (*' - 60) 

60 < * < 80 

[ 0.3 

i > 80 

f 0.0125 

i < 100 

\ 0.0125+ 0.00125 (t 

- 100) i > 100 

1.3 - 0.01 (i - 30) 


0.55 + 0.005 (* - 30) 



I 

( 


0.3 + 0.0008 (» - 50) 2 
0.3 — 0.01 (/ — 50) 

0.0 

18.7 

18.7 + 0.0280 (/ — 90 ) 3 
250.0 


i < 50 
50 < i < 80 
/ > 80 

i < 00 
60 < i < 80 
i > 80 


1 .0-0.0000757 (/- 60)-' 
f 7.25 + 0.015(1-30) *<60 
\ 7.7 » > 60 


(15) 


and i is the orbit inclination in degrees. When / (t;) is negative, the function is to be set to zero. 
The normalized velocity distribution is the fraction of impacts over a 1 Icm/s interval and is given 
by the following expression 


/» 


fjv) 

r f( v ) dv 


(16) 


Results 

It was assumed for this study that ISAT-1 will be orbiting with a perfect attitude, i.e. non- 
rotating and the gravity-gradient boom will be pointing exactly towards zenith. With this assump- 
tion, the zenith angle will be 90 degrees. Figure 2 shows the azimuth angle for the six lateral surfaces 
on ISAT-1. Using Figure 1, the k-factor for ISAT-1 was found to be 7.80. The surface area was 
0.10931m 2 for each of ISAT-l’s six lateral sides. It was also assumed that ISAT-1 will be in an orbit 
with an altitude between 500 and 1000 km and an inclination of 60 degrees. 

Figure 3 shows the expected number of impacts on ISAT-1 from meteoroids of varying mass 
during five years in orbit, it can be seen that only one impact from meteoroids of 0.1 mg and 
larger is expected. At the same time, meteoroids of 1 pg and larger are expected to impact ISAT-1 
about 5.5 times per day. Figure 4 shows that most meteoroids have velocities between 12 and 16 
km/s relative to the Earth. A spacecraft in a 600 km circular orbit will be traveling at 7.6 km/s. 
inclination and altitude have little effect on the meteoroid flux or velocity distribution. 

Figures 5 and 6 show the number of impacts from orbital debris and the debris's velocity dis- 
tribution, respectively. It can be seen from Figure 5 that only one impact is expected from debris 
with diameters 0.3 mm and larger. Using equation 13 the mass for a 0.3 mm piece of debris will 
be about 0.45 mg. It is also expected that debris with diameters of 10 pm and larger (17 ng and 
larger) will impact ISAT-1 about once per day. Figure 6 shows that most orbital debris will impact 
ISAT-1 with velocities between 12 and 15 km/s for an inclination of 60 degrees. 


6 


/ 




Figure 2: Azimuth Angle Geometry of 1SAT-1 
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Log 10 of Meteoroid Mass, gm 


Figure 3: Expected Number of Impacts on ISAT-1 From Meteorids During 5 Years 
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Figure 5: Expected Number of Impacts on ISAT-l From Orbital Debris During 5 Years 
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Debris Velocity Relative to Spacecraft, km/s 


Figure 6: Normalized Orbit Debris Velocity Distribution 
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Introduction 

ISAT-1 will obtain the necessary power to run the bus and payload subsytems by solar arrays 
body-mounted on the lateral faces of the bus. This analysis was done to find out how much power can 
be generated as a function of orbit altitude, sun elevation above the orbit plane, ISAT-l's position 
in its orbit about the Earth, and the efficiency of the solar array. 

Coordinate System 

The analysis was done using Earth-referenced spacecraft-centered celestial coordinates as shown 
in Figures 1 and 2. The orbit plane is the equator of the globe plot and the direction of the Earth, 
the yaw direction, is held fixed. The Earth's disk is the shaded circle in Figure 2. Because the yaw 
axis is always facing nadir, this coordinate frame rotates once per orbit in inertial space about the 
orbit pole; therefore, any objects that are essentially fixed in inertial space during one orbit, such 
as the sun or moon, will appear to rotate about the orbit pole once per orbit. An eclipse will occur 
when the sun passes behind the disk of the Earth. As the Earth travels in its orbit about the sun 
and the longitude of the ascending node of the spacecraft's orbit rotates about the Earth's pole 
due to perturbations, the sun will travel up and down on the globe plot from a maximum of the 
inclination of the spacecraft's orbit plus the obliquity of the ecliptic above the orbit plane to this 
maximum below the orbit plane and back up to this maximum above the orbit plane during one 
year. For example, if the spacecraft’s orbit has an inclination of 60 degrees, the sun will travel from 
a maxi m u m elevation of 83.4 degrees above the orbit plane to 83.4 degrees below the orbit plane 
and then back up to 83.4 degrees above the orbit plane over the course of one year. 

In this coordinate system the normal to a given face is represented by point F on the celestial 
sphere shown in Figure 3. The elevation of the sun above the orbit plane is represented by 0s- 
and the elevation of point F above the orbit plane is represented by 0p. Axs is the azimuth of the 
sun with respect to nadir and Azp is the azimuth of point F with respect to nadir. The rotation 
angle A Az is the difference in azimuths between the sun and point F and varies uniformly from 0 
to 360 degrees for a circular orbit. The angle between point F and the sun is represented by 0. The 
incident solar energy on face F is then 


I = AKcj cos0 (1) 

where 

A = area of face F 

K = solar radiation constant (1358 W/m 7 ) 
tj = overall efficiency of solar panel 


and from the law of cosines for spherical trigonometry 

co6 0 = cos 0p cos 0s + sin 0p sin 0s cos ( A Ar) (2) 

When 0 > 90 degrees the sun is on the "back side” of face F and will not cause any power generation. 
Eclipses occur at sun azimuths given by 


where 


Azotic = arccos 



sin p = 


Rb 

Re + H 


(3) 

(4) 


1 



where 

p — angular radius of a spherical Earth 

Rb = radius of the Earth (6378.14 km) 
H = spacecraft altitude 



Mlnition of a Sptwrlcal Coordinate System on th* Unit Sphere. The point P is at 

■n azimuth of 50 dag and elevation of 35 deg, notmaSy written as (50°. 36*). 


Figure 1: 


Results 

ISAT-1 will be a hexagonal cylinder with the solar panels mounted on the six lateral sides. Five 
•ides will have 0.1024 m 3 of solar cells while the side with the launch vehicle adapter will have 
0.512 m 3 of solar cells. Because not all six sides have the same solar cell area an average power was 
calculated at each point in the orbit. This was done by "rotating” the model about the longitudinal 
axis through 360 degrees by steps of 5 degrees and averaging the power generated. This technique 
was used for generating Figures 4, 5, and 6. 

Figure 4 shows the average amount of power generated in each orbit as a function of sun elevation 
and solar array efficiency at an altitude of 600 km. The average power at each sun elevation was 
determined by calculating the power generated at 30 second intervals during one orbit, including 
eclipse times, and then taking the average of the power generated. Altitudes between 400 km and 
1500 km affect the average power generated very minimally. The linear decrease in power for lower 
efficiencies is from the linear relation in Equation 1. 

Figure 5 shows the power generation profile during one orbit for an efficiency of 0.18 and at 
an altitude of 600 km. As can be seen from this figure, sun elevations near 0 degrees cause great 
flucuations in the power generated at different points in the orbit. This comes from the fact that 
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Figure 2: 


the sun is shining mainly on the nadir and zenith ends of the spacecraft at azimuths near 0 and 180 
degrees. When the elevation of the sun nears 90 degrees (or -90 degrees) the ends of the spacecraft 
are never in sunlight at any time during the orbit and the sunlight is concentrated on three of the 
lateral sides. This causes the power generation curve to become more flat. Also at higher sun 
elevations, the eclipse times become shorter. This is also reflected in Figure 7. 

Figure 6 shows the average amount of energy generated during one orbit. This was calculated 
by numerically integrating the power generated during one orbit with time. As can be seen there 
is more energy generated at higher altitudes. This occurs because the eclipse times are shorter at 
higher altitudes. The curves in Figure 6 closely match the curves in Figure 4. 

Figure 7 shows the time a satellite spends in sunlight and in eclipse as a function of altitude and 
sun elevation. 
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Figure 3: 
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Power, watts 



Figure 4: Average Power Generated During One Orbit at an Altitude of 600 km 
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Time During One Orbit, minutes 


Figure 5: Power Generation Profile For One Orbit at an Efficiency of 0.18 and Altitude of 600 km 


6 






Figure 7: Time Spent in Eclipse and Sunlight 
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ABSTRACT 


ISAT-1 is a satellite being designed, built, and tested by students from the state 
universities of Iowa. The objective of this project was to design a method for determining 
the thermal control needs of the ISAT-1 Satellite. Space is a harsh environment with a 
temperature of four degrees kelvin, and the only method for heat transfer in space is 
through radiation. A satellite in space will be subjected to four main types of radiation: 
direct solar radiation, albedo radiation or solar radiation reflected from the earth, earth 
infrared radiation, and emission from die satellite itself. As the satellite orbits die earth, 
these different radiations will fluctuate causing fluctuations in temperature of the satellite. 
The fact that not all equipment on a satellite can operate at any temperature necessitates a 
method for controlling the temperatures seen by the satellite. This method is called thermal 
control. 


In order to solve die problem of determining the thermal control needs of the ISA 
1 Satellite, a computer program was generated. This computer program was designed to 
used for any satellite configuration by having the user input die properties of the satellite 
be analyzed. The program works by splitting up the satellite into nodes which are 
representative of different parts of the satellite and pieces of equipment on the satellite. By 
inputting orbit information of the satellite, space environment information, and the thermal 
properties of the nodes of the satellite, the program simulates the temperatures the satellite 
will see during it's orbit These te mp eratures are output in the form of high, low, and 
mean temperature for each node during the orbit Output data on the transient response for 
each node can also be obtained from the program to be used for graphing. The thermal 
control system in this program is represented by thermal resistances between each node 
which are input to the program. By inputting different values for the thermal resistances 
(R-values), the temperature of the die nodes on the satellite can be controlled. Once the R- 
values that keep each node within it's respective temperature tolerance are found, they can 
be used to choose materials, such as insulators or conductors, that will be used for the 
thermal control system of the satellite. 


Not enough specific information about the ISAT-1 Satellite was known to perform 
and in depth analysis using this program. However, with the general configuration of the 
satellite, properties were assumed, and the program was ran to determine if it functioned 
property. The results showed that the program worked, and it is suggested, that as more 
information about the ISAT-1 Satellite is obtained, the program is used to find the thermal 
control needs of the satellite. 
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RECOMMENDATIONS 

Engineers use mathematics to understand and predict behavior in the physical 
II & mathematical model created fails to do so, then engineers are no longer 


eraineers — they are simply mathematicians. Therefore, it is desirable to verify a model in 


oner to insure both its relevance to and accurate description of die system for which it was 


created. An improper model will fail both to fully describe the chosen system and to predict 
l behavior, thus, die general rationale for model verification is obvious. 


tt&b 

I The particular rationale for verification of the Thermal Math Model (TMM) for the 
ISAT-1 can be summarized as follows. It is desired to demonstrate that the principal heat 
transfer paths and pertinent heat transfer mechanisms have been correctly anticipate and 


understood. Model verification is also necessary in order to check die validity of the 
assumptions made concerning surface properties, thermal interfaces, and insulation 
effectiveness. Furthermore, it is important to verify the accuracy of the predictions that the 
model mak es as these will be used to determine pe rfor mance reliability 

Three principal methods exist for TMM verification: verification by independent 
analysis, verification with flight data, and verification by test Verification by independent 
analysis involves a separate design team conducting a parallel analysis of the system. In 
this particular case, this method is not an option because this team is the only team which 
will be involved with this phase of the design for ISAT-1. In large scale satellite 
operations, often several, subsequent vehicles will be launched which are similar to one 
another in design and/or orbit This allows for the es tablishme nt of a performance histor y 
constructed from the knowledge gained through flight data. ISAT-1 is the first of its 
particular design to be launched; therefore, there is no previously acc umulated flight data. 
Due to this fact verification with flight data is not a possibility. The remaining method is 
verification by test This has become the chosen method for verifying die TMM of ISAT-1. 

Two chief test types have been selected by this team to be recommended to the 
ISAT-1 committee. They are namely, a Thermal Balance Test and a Thermal Vacuum Test 
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Detailed testing procedures have as yet to be decided. However, the following basic 
information concerning die tests has been acquired. 

In effect, the thermal balance test is the qualification of die Thermal Control System 
(TCS) as it verifies the analytical thermal model This test will be used to demonstrate the 
functional capability of die on-board thermal control hardware. Thermal balance testing is 
vital in establishing the ability of the TCS to maintain components, subsystems, and space 
vehicles within specified temperature tolerances for critical mission phases. 

A thermal vacuum test is not the same as a thermal balance test While the balance 
test is designed to simulate actual mission operation, the vacuum test is merely a hot / cold 
- soaking cycle which is used to subject the satellite to its anticipated environmental 
extremes in a worst-case scenario. 

At present very little information is available concerning die payload for die satellite. 
In order far the design, or the computer program to be fully implemented, more specific 
data must be defined concerning ISAT-1. Complete information needs to be compiled 
which must include detailed information for all the intended on-board components (size, 
location, material composition, mass, power requirements, thermal tolerances, and duty 
cycles). 

Once this information is known, the existing nodal network can be expanded to 
accommodate the true configuration for ISAT-1. The design described herein can then be 
used to generate the solution for the Thermal Control System of the satellite. The solution 
will contain all the TCS requirements, specifically the R-values, which can be manipulated 
through material / device selection to satisfy the thermal criteria for the payload 
components. 

CONCLUSIONS 

To date the design team for the TCS of ISAT-1 has succeeded in producing a 
design to determine satellite thermal control needs. Any configuration can be 
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ac^ommod 


ated, and die design has the capability to determine the TCS requirements for the 


particular configuration chosen. The requirements are in the form of R-values far thermal 
distance. 

I In this preliminary analysis, a hypothetical configuration based upon the limited 
component information currently available has been used. The primary intent of this 
arjalysis was to prove that the present design is capable of producing accurate and 
meaningful results. The R-values are indeed meaningful as they represent design elements 
which can be effectively generated and modified to control and regulate thermal criteria for 
the components. The accuracy of the R-values lies in the fact that they have been derived 
from the appropriate mathematical models in which all the relevant physical phenomena 
were considered. 

A major accomplishment of the design and analysis is that the need for a TCS for 


ISAT-1 has been confirmed. The analysis run for the hypothetical configuration has 
in d i cate d the existence of potential problem areas with regard to thermal overload. The 
computer program has shown its ability to provide a solution to this problem by generating 
die R-values necessary to control the thermal overloads through appropriate manipulation 
fix., material / device selection and component location). 


INTRODUCTION 

The Iowa Satellite Project, ISAT-1 will not just be another co mmuni cations 
satellite, it will be an educational tooL The ISAT- 1 will be designed, built, operated, and 
used by students, and the general public throughout Iowa. Although all three of the Iowa 
Universities will be involved in the project, only die University of Iowa and Iowa State 
University will design the satellite. The University of Northern Iowa's main function will 
be to provide educational material and projects from the satellite information for use by 
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grade school age students and up. Iowa State will design orbit and the space systems part 
of the project, while University of Iowa will be designing the thermal control system or 
TCS for the satellite. More background about the satellite, including the equipment and 
experiments on board the satellite can be found in the Appendix A. 

The objective of this project was to analyze the temperatures die ISAT-1 satellite 
will be subjected to, and then design die method for determining die thermal control needs 
of die satellite. 

The thermal control system of a satellite is needed in order to maintain structural and 
equipment integrity over long periods of time. The fundamental requirement of die thermal 
control system is to allow optimum performance of components by controlling the 
temperature in which they operate. If it were feasible to operate the equipment and 
experiments at any temperature, there would be no need for a thermal control system. 
Reliable long term performance of most spacecraft components takes place at or near room 
temperature (293 K). Most of die equipment on board die satellite have certain temperature 
ranges that they can operate within, without thermal failure. (Fisher, 1993) However, as 
the satellite orbits the earth, the temperature in space is 4 K. This does not mean the 
satellite will be at 4 K because it will be subjected to different heat fluxes. 

The only means of heat transfer in space is through radiation and conduction. 
There are four types of radiation that a satellite will be subjected to in space: radiation from 
the sun, albedo radiation (reflected sunlight off the earth), earth shine (infrared radiation 
from earth), and satellite e mis sion. It can be seen from Figure 1, that the radiation from the 
sun and the albedo radiation will vary as the satellite passes in and out of the shadow of the 
earth. Therefore, the satellite's temperature will fluctuate with time. These fluctuations in 
temperature as well as the temperature in space make up the need for a thermal control 
system of a satellite. 



5 



Figure 1: Radiation in the Space Environment 
Design of the thermal control system for die ISAT-1 can be divided into five parts: 
define the spacecraft environment, define thermal control options, select and analyze 
thermal configuration, and select components.(Iowa, 1992) In order to perform a full 
analysis, concurrent engineering has had to take place between die state universities, and, 
as information is lacking due to this aspect of the project, personal goals for the purpose of 
the class have been set Some information on the satellite's structure and components has 
been obtained, and it was decided to use this information to begin an analysis of the 
temperatures the satellite would see in space. To perform this, a computer program has 
been written to simulate the satellite in orbit Based on the satellite's material properties the 
program analyses the temperatures of the components of die satellite as it orbits the earth. 
The thermal control needs of the satellite are found by using this program. This paper will 
describe the analysis of the satellite performed and results found to dat e , 

THEORY: 

The spacecraft thermal balance between cold space (4K) and solar, planetary, and 
equipment heat sources is the means by which the desired range of equipment and 
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structural temperatures arc obtained. In this analysis, the thermal balance was achieved in 

the form of a FORTRAN program, that determined the satellite temperature as a function of 
time. 

Before the formulation of the thermal balance program, die im p or t a nt characteri stics 
of die space environment in conjunction with the mathematical techniques to solve the 
thermal balance had to be considered. In the following discussion, the characteristics of the 
space environment will be limited to those encountered by a smyJifo in low earth orbit. 

The thermal conductivity of the Earth's atmosphere below an altitude of 90 km is a 
function of temperature change and is independent of variations in pressure. Above 90 Wn 
the molecular mean free path, or the distance a particle may travel before colliding with 
another, becomes comparable to the distance in which a temperature change significantly 
varies, and therefore thermal conductivity becomes pressure dependent At an altitude of 
300 km, the molecular mean free path increases to where thermal conductivity is negligible. 
Therefore, convective heat transfer may be neglected above 300 km.(Fortescue, 1991) 

Hence, heat transfer between the satellite and die space environment above 300 km 
is limited to radiation heat transfer only. Sources of radiation that may effect the physical 
temperature of the spacecraft are direct solar radiation, solar radiatio n reflected from nearby 
planetsfalbedo radiation), thermal energy radiated from nearby planetsfinffarcd radiation), 
and radiation out to deep space from the spacecraft itself.(Fischer, 1993) 

Once the sources of radiation have been determined, the radia tion ne eded to be 
converted into heat transfer. When converting the radiation incident upon the satellite into 
heat transfer, several parameters must be understood. These parameters consist of the 
radiation intensity, degree of collimation, and die enrissivity and absorptance of the c«t»mtr 
skin. It should be noted that these parameters vary for each type of radiation and are 
therefore discussed individually. 

Solar radiation shall be considered first because it has the highest intensity. This 
solar radiation intensity, Sq, from a point source, for any distance D from the sun may be 


1 
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expressed as shown in Equation 1, where the total power output of the Sun, P, is 
3.8x1025 W. 

So= ^ (1 > 

The radiation from die Sun may be considered to be collimated, or in other words, 
entirely made up of parallel waves. This collimation may be assumed due to the radiation 
originating from a very large distancc(D=1.50xl0^ 1 m) from the earth. The resulting value 
of Solar radiation is 1430 W/m^ for the closest approach to the Sun, which represents the 
worst case scenario. (Fischer 1993) 

In determining the solar radiation heat transfer, 0^ . the product of the solar 
intensity, area the radiation is incident upon, A^i^, and absorptance, ct, as shown in 
Equation 2, must be found. 

Qsolar = CtAsolarSo (2) 

Another type of external heat transfer that must be considered is albedo radiation 
heat transfer. The albedo radiation intensity may be calculated as a fraction of solar 
radiation, this fraction varies with earth surface conditions. This fraction may be 
determined by a weighting factor that varies 0.8 to 0.03 for a cloud shrouded or vegetation 
covered earth, respectively.(Fischer, 1993) 

Another consideration in determining the intensity of this type of radiation is the 
noncollimation caused by the reflection properties of the earth. Not c ollimated radiation 
may be taken into account by the use of a view factor. A discussion of the use of a view 
factor is imperative in determining the albedo radiation incident upon the satellite and 
therefore follows. 

A view factor represents the variation in radiation inc ident on a surface du e to 
noncollimation of radiation. In order to illustrate the use of a view factor, there are two 
limiting cases that need to be considered. For the first case, where the satellite would 
receive the most radiation, the surface of die earth may be thought of as a flat plate 
extending infinitely in all directions. This infinite flat plate would reflect a fraction of the 
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solar radiation incident upon it, which has been previously discussed by the use of a 
weighting factor. The reflected radiation would be non collimated or more specifically 
reflected in all directions. With the satellite above this infinite flat plate, the radiation 
incident on a surface of the satellite parallel to the plate, would not only come from the area 
directly below the satellite, but also from areas far outside. In other words, the radiation 
from the flat plate is emitted in a disperse pattern, that causes an increase in the radiation 
flux of the parallel satellite surface. Furthermore, radiation incident on the sides of the 
satellite due to this disperse pattern, would also add to the radiation flux. 

For the second limiting case, the radiation emitted from the infinite flat plate would 
be collimated, therefore only the parallel surface of the satellite, facing the infinite plate, 
receives radiation. This second limiting case represents the least amount of radiation 
incident upon the satellite. The actual case is between these two extremes, the reflected 
radiation is non collimated and the surface of the earth is neither flat or infinite. By using 
the view factor, this variation caused by non collimated radiation may be handled. 

The intensity of radiation, J a , received by the satellite therefore depends on the 

product of view factor, F, solar radiation intensity and weighting factor, a, as shown by 
Equation 3. 

Ja = SoaF (3) 

At last, the albedo radiation heat transfer, Qaibedo> may be found by the product of 
the albedo radiation intensity, Ja, the area of the satellite which has radiation incident upon 
it, Aaibedo» and the absorptance of the satellite skin, as shown in Equation 4. 

Qdbedo * aAdbedoJ* (4) 

The third external heat transfer due to infrared radiation that must be considered, is 
emitted from the earth acting as a black body radiator. Intensity of the infrared radiation, 
Jp, is determined by the Stefan-Boltzmann fourth power law, which relates die product of 
temperature, T, to the fourth power, and the Stefan-Boltzmann constant, o, with the 
radiation intensity, as shown by Equation S. 
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Jp-tfT* (5) 

Therefore, to determine the infrared heat transfer, the product of the infrared radiation 
intensity, the area of the satellite which has radiation incident upon it and the emittance of 
the satellite skin, e, as shown in Equation 6, must be found. 


Qinfra s CAmfraJp (6) 

The final heat transfer with die external space environment is from the satellite out 
to deep space. This heat transfer occurs because the satellite, acting as a black body 
radiator, emits radiation out to deep space. As previously discussed, a black body 
radiator’s radiation intensity may be found by the Stefan-Boltzmann fourth power law. In 
determining the heat transfer to space, Qintemal. the product of the intensity of the black 
body radiation, the total emitting area of the satellite, A ^nim. die Stefan-Boltzmann 
constant, and the emittance of the satellite skin, as shown in Equation 7, must be found. 

Qextemal * ^saldhlet^gatelliie-T^apace) (7) 

Once the important characteristics of die space environment have been investigated, 
a thermal balance of the satellite is determined by setting the heat transfer to deep space 
from the spacecraft, equal to the solar, albedo and infrared heat transfer. The resulting 
polynomial is shown in Equation 8, 

EAradiantafT 4 satellite — T 4 space) = CtAsolarSo 4- CtAalbedoS albedo ^ 

+ eAinfraredS infrared 

where the term on the left hand side is die heat transfer to deep space, and the terms on the 
right hand side are the solar, albedo, and infrared heat transfer, listed first, second and 
third, respectively. The results acquired by this thermal balance may be considered a 
steady state solution to the energy balance. 
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The radiations from the different sources, incident on the satellite, vary with the 
satellites position as it orbits the earth and also as the earth orbits the sun. Hence, a 
transient analysis was applied to account for the effects of die changing incident radiation 
on the satellite, with respect to time. This transient analysis was performed by equating die 
product of die satellite bus mass, the specific heat of die bus, and the change in t emp er a ture 
with respect to time to the sum radiation heat transfers and internal heat generation, 
Qintemal, shown in Equation 9. It should be noted, that die internal heat generation will 
initially be assumed to be constant 

dT 

mcp— = aAaolarSo + aAa lb edoSalbedo + EAinfraredS infrared 

dl <9) 

+ Qintemal — EAradiantoCT 4 satellite — T 4 space) 


The resulting equation is a first order nonlinear differential equation, whose 
solution would be very intensive if solved by hand. Using a numerical solution, such as 
the Fourth Order Runge-Kutta method, the solution may be simplified. With utilization of 
a computer program the solution may be found for many different times for many different 
satellite characteristics, in much less time than one hand calculations. 

In order to employ the Runge-Kutta method. Equation 10 had to be transformed 
into the form 

T « F(T,t) (10) 

where T is the time derivative of temperature and F is a function of tune and temperature. 
The general solution of this method takes the form shown in Equation 1 1. 


yi+i 


yi+ 


XWjKi 

XWj 


( 11 ) 


As may be seat from the general form die value of the variable y is calculated for each time 
increment with the previous y value plus die summation of the weighted gradients, WjKj, 
divided by the summation of the weighting factors, Wj. Due to the need for a previous 
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value y, an initial condition must be specified. It should be noted that as j increases, the 
order of the Runge-Kutta method increases, increasing die accuracy but also increasing die 
solution time. For this analysis, the fourth order Runge-Kutta method, which is the most 
commonly used order, providing adequate accuracy and efficiency, was used. 

Now that an external analysis has been completed, an internal analysis may be 
performed. To begin this analysis, the internal components and equipment may be 
modeled as nodes with individual temperatures. The heat transfer between these nodes is 
related by a thermal resistance. To illustrate the relation between two components. Figure 2 
has been included. In this figure, R12 is the thermal resistance between component 1 and 2 


R 

Node 1 Node 2 

Figure 2 : Thermal Relation Between Two Internal Components 
represented by node 1 and 2, respectively. To determine the internal heat transfer between 
these two component, the quotient of temperature differential an thermal resistance, as 
shown in Equation 12, must be found. 


HAArH 


T 2 


Q12 = 


T1-T2 

R12 


( 12 ) 


Now that the fundamental concepts of heat transfer between the external 
environment and satellite, and heat transfer between the internal components have been 
determined. These external and internal analysis may be combined into one analysis, 
shown by Figure 3, that takes into account all heat transfers. By relating the external heat 
transfers from space to the satellite sldn and the satellite sldn to the internal components of 
the satellite, a combined analysis may be performed to determine the temperatures of 
internal components. 
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Figure 3: Combined External and Internal Analysis 


Computer Program 

The analysis outlined in the previous section, leaves die problem of determining a 
method to solve far the temperature of each node on a satellite over a range of time. It was 
d e cided that this problem would be handled easiest by using a computer program to solve 
for die temperatures on the satellite. 

The nodal analysis as outlined in die previous figure, shows, that the space 
environment affects only the nodes that are representative of the skin of the satellite. The 
satellite sldn is represented by three nodes, one for the top, bottom, and sides of the 
satellite. The temperatures that the external nodes exhibit, in turn, will be transmitted to the 
internal components of the satellite depending on the thermal contact between the nodes. 
Also, the internal nodes can affect die temperature of the outer shell nodes of the satellite 
through die internal nodes respective heat generation terms. Due to these many interactions 
between different factors which determine the temperature of each node of the satellite, it 
seems that a computer program could get rather complex. However, since the satellite is 
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broken up into nodes, arrays can be used to describe the properties of each of these nodes, 
and how they interact with the space environment and each other. The array style of the 
problem leads to a computer program which only needs to solve one function (the 
temperature function). This temperature function can now be represented in an array 
format as shown in Equation 13. 

F = liT = ^^^(QtafO+QsotO+QalbCO+QinfO)* 

I 2 - ( ^ ' ~ T(J) -e(i)A(i)g(T(i) 4 -To 4 ) 
i j ROj) 

In this equation, "i" represents the node number being evaluated. The other 
variables represent either properties of the satellite that would have to be input to the 
program, or values for heat flux terms that the program must find for each node. To 
understand this better, each term in the equation will be described. The "T" and "t" 
represent temperature and time respectively. These are the values that the p rogram will try 
to find for each node. The mass (m), specific heat (Cp), and internal heat generation (Qint) 
arc values that describe properties of nodes, and must be input to the program. The next 
terms (Qsol, Qalb, Qinf) represent the different types of radiation from the space 
environment that will be incident on the external skin G f the satellite. For this reason, these 
variables arc set to zero for internal nodes. Also, since the values of these radiations vary 
with the particular time position of the satellite in orbit, the values will be solved for by the 
program dependent on the time. The next term describes the heat flux due to the thermal 
contact between nodes. This term is very important because it is a function of the thermal 
resistance between nodes, which represents the thermal control system of the satellite. It 
should be noted that this term has to be calculated for the interaction between every node on 
the satellite and the node being evaluated. The last term describes the Stefan-Boltzmann 
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fourth power law for radiation emission from the satellite skin. This term is set to zero for 
any of the internal nodes. 

The temperature function equation and the information outlined in the theory section 
of this paper were used to formulate a computer program that would solve for the 
temperature of each node on the satellite as a function of the time in it’s orbit The name of 
the program is "isaLftn," and for a detailed description cm die operation and background of 
die program see Appendix B. This appendix contains a program handbook far "isaLftn." 
The general form of the program is a Runge-Kutta numerical differential equation solver 
which was modified to handle arrays. The differential equation solver evaluates the 
temperature of the satellite in second time steps. At each time step, every node is evaluated 
by calling different subroutines to find the value of the terms in the temperature function 
equation. The temperatures of each node of the satellite are then found, and the rime step is 
incremented. The program was initially formulated to analyze the ISAT-1 satellite, but has 
been modified to accept any satellite configuration that falls under the assumptions used in 
generating the program (See the Program Handbook). Also, the program could be 
modified in any way needed to analyze any satellite configuration. 

The program that was generated to evaluate the temperature function for each node, 
needs certain input parameters. These input parameters are entered in die form of an input 
file which contains orbit information, space environment information, and the thermal 
properties of the satellite. An example of the input file format can be found in the Program 
Handbook in the Appendix of this paper. 

The orbit information needed is information about the period of the orbit and the 
amount of time the satellite will be eclipsed by the earth. The space environment 
information needed describes the intensity or amount of the different types of radiation that 
will be incident upon the satellite. The first radiation is the solar radiation which typically is 
constant in space near the earth. The second type of radiation is the albedo radiation, or die 
radiation from the sun that is reflected back off of the earth. The input for this radiation is 
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represented by a decimal fraction of the solar radiation. The thud type of radiation is the 
earth shine, or the infrared radiation that the earth emits due to the fact that is a thermal 
emitter and follows the S tefan-Boltzmann fourth power law. 

The first satellite properties that must be input are the length (or height) and 
effective diameter of the satellite. Next, the initial time to begin die analysis and die initial 
temperature of the satellite must be input For ease, these input values are usually just 
assumed to be an initial time of zero, and an initial temperature of 273 K. These values are 
needed due to the fact that die program uses an ordinary differential solver and must have 
initial conditions to begin the analysis. The next satellite properties for the input file deal 
with the nodal analysis of the satellite. First, the number of nodes to be simulated must be 
input By allowing the program user to input the number of nodes, many different satellite 
configurations can be handled by the program, and as a thermal control analysis becomes 
more developed, the program can handle this expansion. Then, the mass, specific heat 
and internal heat generation of each node must be input The next values describe the 
properties of the skin or satellite covering. The values needed are the absorptivity (alpha) 
and emissivity (epsilon) of the different sides of the satellite (top, bottom, and sides). 
Also, two view factor values are input These view factors describe the view factor 
between the sides of the satellite to earth, and the view factor between the bottom of the 
satellite to the earth. These view factors are needed in determining the amount of albedo 
and earth shine radiation incident upon the satellite. 

The last properties needed to describe the satellite in the input file, are the thermal 
resistance, or R-values between the components and between the components and sides of 
die satellite. The units for the R-values is in Kelvin's per watt An R-value is needed for 
each interaction between each node. This means that if two components are not near each 
other, or are thermally isolated from each other, a large R-value such as 109 can be input 
Conversely, if two components are in perfect thermal contact such as between the sides 
and top and bottom of the satellite, the R-value should be very small. It is suggested that 
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some kind of preliminary calculation for these R-values be pe rfor med because die program 
will not handle extremely small R-values or values less than 10 ~ 2 . The values in this 
situation can be calculated by equation 14.(Chapman, 1984) 

R -r X <“»> 

where: Dx * conduction distance, m 

k = thermal conductivity of the material, WAn-K 
A * the area of conduction (or touching), m^ 

Due to die set up of the program, an R-value is also needed for die thermal contact between 
a node and itself, these R-values can just be set to one in the input file. 

As was mentioned, these R-values are important because they represent die thermal 
control system of the satellite. It should be noted that these R-values represent only a 
passive thermal control system (one without moving parts or a working fluid). For this 
reason, this program will only evaluate satellite temperatures for passive thermal control 
systems. The program, however, can be used in two different ways. The first way the 
program can be used is to determine if a previously designed thermal control system for a 
satellite is correct For this type of analysis, the R-values of the materials for the thermal 
control system can be input into the program, and the high* low, and mean temperature of 
each node of the satellite will be output These temperatures can then be compared to the 
temperature tolerances set for each node (or piece of equipment) to see if the thermal control 
design is correct The second way the program can be used is to design a thermal control 
system. This is performed by assuming R-values for the satellite, running the program, 
and then modifying the R-values until the temperatures of each node fall within their 
respective temperature tolerances. If there is not a combination of R-values that satisfies 
die thermal criteria of the satellite under this situation, this shows that there may be a need 
for sense kind of active thermal control system. 
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The program was tested by performing an analysis of the ISAT-1 Satellite. At the 
time the program was developed, not enough specific information about all die satellit e and 
components were known to perform an in depth analysis. However, the main 
configuration of the satellite was known, and by assuming the properties for the 
configuration of the satellite, the program was run and tested to see if it would produce 
results. 

RESULTS AND DISCUSSION 

For the hypothetical configuration (Fig. 4), the the rmal resistance matrix (6 x 6) 
contains (36) elements. Symmetry considerations (Le., Rjj = Rj0 reduce this number to 
(21). Further reduction of values is achieved by the geometry and known mechanical 
behavior of the system. For example, the R-value between the top and the bottom of the 
structure is defined as infinity due to their thermal isolation from one another. Also, the R- 
value between the top and the sides is defined as B 0 due to their identical mater ial 
composition and direct contact which results in complete the rmal c ommunic a tion . S imilar 
definitions exist which reduce the number of final R-valnes to be designed for to only (7). 
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Figure 4: Nodal Configuration 

In order to accommodate the actual, potential temperature responses of the nodes 
during orbit, two worst-case scenarios were investigated. The first of which was to 
si mul a t e a situation in which all the nodes were in complete thermal isolation. This was 
achieved by setting the design R-values to 00 . The temperature responses for this 
configuration after a twelve-orbit simulation are displayed in Table 3. The design R-values 
are also indicated. 


Table 3: Node temperatures after 12 orbits - 
Complete thermal isolation 


Satellite 

Component 


KBIISI 

Mean 
Temp. (K) 

Sides 

304.8 

280.7 

292.7 

Top 

301.0 

279.7 

290.3 

Bottom 

302.0 

280.4 

291.4 

Level 1 

360.9 

277.5 

349.2 

Level 2 

342.2 

323.7 

332.9 

Level 3 

289.1 

284.8 

286.9 


®14 = ®15 s ®ir R 2r R 34 s *4S s ®St * 


There are several conclusions which can be drawn from this result The first is that 
when the nodes are in thermal isolation, no heat dissipation is allowed; therefore, as the 
number of orbits increases, the mean temperatures will all steadily increase. This has the 
potential to cause thermal overloads in the nodes (depending on their individual temperature 
tolerances). Although very little detailed information is available on component temperature 
tolerances, information about the batteries is known. Specifically, batteries need to be 
maintained at around 273 K. The batteries are contained within level 3 which indicates a 
mean temperature of 349.2 K. Therefore, a second conclusion to be drawn is that the 
batteries are experiencing thermal overload. 
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The results after 12 orbits for the opposite worst-case are displayed in Table 4. This 
case s imul a t es a situation in which the nodes are in complete thermal c ommnnirarinn with 
one another. This scenario was achieved by setting the R-values to *0. 


Table 4: Node temperatures after 12 orbits - 
Complete thermal communication 


Satellite 

Component 

WBS-atSM 

■KESiH 

mBSSiSm 


Sides 

307.2 

301.4 

304.3 

Top 

306.9 

301.4 

304.1 

Bottom 

307.0 

301.5 

304.2 

Level 1 

307.1 

301.7 

304.4 

Level 2 

307.1 

301.7 

304.4 

Level 3 

307.0 

301.6 

304.3 


®ir^ij = ®ir®2r®34 s ®4j s ®« " 0 


The main conclusion here is that when the nodes are allowed to fully dissipate heat 
to their surroundings (and surrounding nodes) the satellite will eventually reach the rmal 
equilibrium. The potential problem with respect to component temperature tolerances is that 
temperature differentiation among nodes is prevented which can also result in potential 
thermal overloads. This phenomenon is evidenced by the fact that the batteries are srili in a 
state of thermal overload at 304.4 K. 

Obviously, the two aforementioned cases are unrealistic. In reality one would 
expect the proper R-values to lie between 0 and °°. Owing to the fact that no detailed 
information is available to direct the convergence upon the proper R-values, a trial and error 
method was employed with die goal to drive the temperature of the batteries as close to 273 
K as possible. The best result possible was achieved with the R-values displayed in Table 
5. 


Table 5: Node temperatures after 12 orbits 
Trial & error best case 
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Satellite 

ComDonent 

High 

Temn. flO 

Low 

Temn f K) 

Mean 
Temo. fKl 

Sides 

309.5 

295.7 

302.6 i 

Top 

306.2 

298.0 

302.1 

Bottom 

307.3 

297.2 

302.2 

Level 1 

307.6 

299.7 

303.7 

Level 2 

307.5 

304.9 

306.2 

Level 3 

306.4 

302.9 

304.7 


R , 4 =0.1 R, 5 =lE +6 R, 4 =lE +6 

R24 = 0.1 R 24 = 0.1 R 45 =lE +6 

R s< =0.1 

As indicated, the lowest temperature possible, with this configuration, was 303.7. 
Since this is still in violation of the tolerance for die batteries, it is concluded that a TCS is 
necessary for ISAT-1. Manipulation of component location and materi als / devices 
used to affect the heat transfer among the components will be necessary in order to direct 
the convergence upon the proper R-values necessary to maintain proper temperature 
requirements. This is the next extension of the TCS design. 


I 
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1.0 Introduction 

1.1 Program Background 

isatftn is a fortran program which simulates the temperature fluctuations a satellite 
will be subjected to as it orbits the earth in space. The program is set up so that the 
properties of the sate l li t e, as well as orbital information, arc input, and die temp erature 
response of different parts of the satellite is output The progra m can be usedtn two 
general ways. First the program can be used to determine the thermal control needs of a 
satellite. Or, secondly, it can be used to determine if the designed thermal control s ystem 
of a s a t e ll i te will keep the components of the satellite in their respective temperature 
tolerances. 

As a satellite orbits the earth, it will pass in and out of the shadow of die earth. 

This causes the different types of radiation incident upon the satellite to vary with Him* 
Since radiation is the only method of heat transfer in space, these fluctuations in the 
inc id e nt ra di a tio n will cause die temperatures seen by die skin, or covering, of the 
to flu c tu a t e. The skin temperature fluctuations will he t ranslated imn fluc tuations in 
temperature of die internal components of the satellite. Thermal control of a is 

n eeded in order to keep all of die different components of a satellite within temperature 
tolerances for the safe operation of the equipment Since, a satellite is generally comprised 
of many different components, many different temperature needs may exist on one satellite 
This brings about the need for an easy method to analyze die temperature of the 
components of a satellite as it orbits die earth. ISATJTN was originally developed to 
determine the thermal control needs of the ISAT-1 Satellite, but, this program is 
for use with other satellites. 

In order to use this program for analyzing the temperatures of a satellite, the satellite 
must be able to fall under some of the assumptions used in generating die prog ram These 
assumptions could be modified by making changes to the program, but, without changing 
die program the satellite must fall within these specifications, first, die satellite must have 
an orbit height of at least 300 km. This is due to the fact that this program does not take 
into account the effects of aerodynamic heating, which are existent for orbits below 300 
km. Second, it was assu m ed that the satellite must be able to be nyxlel ed as a cylinder with 
an effective diameter. Most satellites have die general shape of a cylinder (such as 
hexagonal or other sided geometries) and spin about their center axis. Third, the 
must orbit with it's bottom side pointing toward the earth. The pro gram was set up this 
way due to the orbit configuration of the ISAT-1 satellite. This assumption may limit the 
use of the program for other satellites, but die progra m can be mnHifiiyt fairly easy by 
following the outline of the program shown later in this manual Other assumptions for 
this program will deal with the input to the program, or die modeling of the gaming 
properties, which the program needs, to perform calculations. These assumptions will 
become apparent in the next sections of die manual which describe die program use and 
show an example of die program capabilities. 

Once an input rile is prepared for the ISATJTN program, the program will run, 
and two different options for output are available. The first option, is a display of the high, 
low, and mean temperatures in degrees kelvin, that each satellite component will see 
during its orbit The second option is a data output of the temperature response (kelvin) 
with time of any component on the satellite. This data can be used to make a plot of die 
temperature response of different components on the satellite during orbit to study carefully 
die pattern of temperature change of particular components. 
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12 Program Method 

The program is set up to perform a nodal heat transfer analysis of the satellite. 

Each component of the satellite represents a node. This means that there will always be at 
least three nodes in the analysis because die sides, top, and bottom erf die satellite each 
count as one node. The program simulates the temperature interaction of the space 
environment with theses three nodes and then uses a finite difference methyl to determine 
die interaction of the outer sltinof the satellite with the internal components, and the internal 
components with each other. The thermal control system, in this program, is represented 
by thermal resistance values (R-values) between each internal node and between internal 
nodes and the external nodes or sldn. These R-values are representative of the material 
properties of either an insulation or some kind of conductor which is to be used in the 
thermal control system. In this sense, this program can only be used to simulate a passive 
thermal control system, or one in which there are no mechanical or flowing parts. It non 
however, point out problem areas that may need active thermal control to maintain the 
correct temperatures. 

As was mentioned, the program can be used in two ways. The first method is to 

determine die thermal control needs of a satellite. This can be accomplished inputting 
certain R-values into the program, running it, and then checking the output to determine if 
the input R-values are sufficient to keep each node within it's temperature tolerance. If the 
R-values arc correct, then materials can be chosen to use on the satellite that satisfy the 

p articul ar R-values. The second method is used to check to see if a thermal control design 
is correct. In this method, the properties of die satellite and the R-valucs that are 
characteristic of the thermal control system design can be input into the program. The 
temperature output will then show whether or not the R-values (or thermal control design) 
will keep the satellite components at their respective temperatures. 

The next section describes how to run the program, the input that is and the 

input format, and the units of all of the input characteristics. 

2.0 Program Operation 

The program operation begins by compiling the program. This step is completed 
by typing "f77 isatftn" and pressing return (notice isatftn is in lower case letters). Now, 
the compiler has produced a file called "a.out", and by typing H a.out" and pressing return 
the program will run. But, before the program can be ran, an input file has to be prepared 
which contains orbit information, space environment information, and die thermal 
properties of the satellite. 

The orbit information needed is information about the period of the orbit, and the 
amount of time the satellite will be eclipsed by the earth. The units for die input time is in 
minutes and the eclipse time input is represented by entering die times that the satellite will 
enter the shadow and exit die shadow of the earth. The space environment information 
needed describes the intensity or amount of die different types of rartinrirm that will be 
incident upon the satellite. TTie first radiation is the sola - radiation which typically is 
constant in space near the earth. TTie second type of radiation is the albedo radiation, or the 
radiation from the sun that is reflected back off of the earth. Hie input for this radiation is 
represented by a decimal fraction of the solar radiation. The third type of radiatio n is the 
earth shine, or the infrared radiation that the earth emits due to the fact that is a thermal 
emitter and follows the Stefan-Boltzmann fourth power law. The units for the solar and 
earth shine radiations must both be in watts per meter squared. 



The first satellite properties that must be input are the length (or height) and 
rafec&ve diameter of the satellite. The units on these properti es must be in meters. Next, 
the initial time to begin die analysis and the initial temperature of die satellite must be input 
For ease, these input values are usually just assumed to be an initial time of zero, and an 
initial temperaturc of 273 Kelvin. These values arc needed due to the fact that the program 
uses an ordinary differential solver and must have initial conditions to begin the analysis. 
The next satellite properties for the input file deal with the nodal analysis of the satellite. 
First, the number of nodes to be simulated must be input Then, the mass, specific heat, 
and internal heat generation of each node must be input The units on these values are 
kilograms, joules per kilogram degree kelvin, and watts, respectively. The next values 
describe die properties of the skin or satellite covering. The values need e d are the 
absorptivity (alpha) and emmissivity (epsilon) of the different sides of die satellite (top, 
bottom, and sides). Also, two view factor values arc input These view factors describe 
die view factor between the sides of the sa t elli t e to earth, and the view factor between die 
bottom of the satellite to the earth. These view factors are needed in dete rminin g the 
a mo u nt of a lb edo and earth shine radiation incident upon the 

The last properties needed to describe the satellite are the thenmul resistance, or R- 
values between die components and between the components and sides of die The 

units for the R-values is in kelvins per watt. An R-value is needed for each infraction 
between each node. This means that if two components arc not near each other, or arc 
thermally isolated from each other, a large R-value such as 10 6 should be input 
Conversely, if two components are in perfect thermal contact such as between the 
and top and bottom of the satellite, die R-value should be very small, it is sug gested th at 
some land of preliminary calculation for these R-values be performed because the program 

will not handle extremely small R-values or values less than Hh 2 . The values in this 
situation can be calculated by equation l.[l] 



where: Ax = conduction distance, m 

k = thermal conductivity of the material, W/m-K 
A = the area of conduction (or touching), m 2 

Due to the set up of the program, an R-value is also needed for the thermal contact between 

a node and itself, these R-values can just be set to one in the input file. 

The format of the input file is set up so that each number has a field of thir t^n 
characters. The format for the input file is shown on die next two pages for a nodal 
analysis of six nodes. The numbers that are shown on this page need to be included in the 
input file, and should not be changed. Also, the input values should hr entered ns decimals 
or exponential form. 
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Input File Format 
The input file should start at the top of the page of die file: 

Orbit Time (min) Shadow In (min) Shadow Out (min) 

Solar Radiation (w/m^) Albedo Factor Infrared Radiation (wAn^) 

Diameter (m) Length (m) Initial Time = 0 

Initial Temperature (K) 

Number of Nodes (Integer value including the different sides of the satellite) 

Mass Sides (Node 1) (kg) 

Mass Top (Node 2) 

Mass Bottom (Node 3) 

Mass Node 4 
Mass Node 5 
Mass Node 6 

Specific Heat (cp) Sides (j/kg K) 

CpTop 
Cp Bottom 
CpNode4 
Q) Node 5 
Q>Node6 

0 Internal Heat Generation Sides (Qint) (W) 

0 Qint Top 
0 Qint Bottom 
Qint Node 1 
Qint Node 2 
Qint Node 3 
Alpha Sides 
Alpha Top 
Alpha Bottom 
Alpha Node 4 = 0 
Alpha Node5 = 0 
Alpha Node 6 = 0 
Epsilon Sides 
Epsilon Top 
Epsilon Bottom 
Epsilon Node 4 = 0 
Epsilon Node 5 = 0 
Epsilon Node 6 = 0 
View Factor Sides 
View Factor Top = 0 
View Factor Bottom 
View Factor Node 4 = 0 
View Factor Node 5 = 0 
View Factor Node 6 = 0 
R 11 = 1 

R 12 = small value 

R 13 = small value 

R14 

R15 

R16 

R22 = 1 

R 23 = 10 6 
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\ R24 = 10 6 

R25 = 10 6 
R26 
R33= 1 
R34 

R 35 = 10 6 
R36= 10 6 
R44 = 1 
R45 

R46 = 10® 
R55=l 
R56 
R 66= 1 
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Once this input file is prepared the following steps are needed to complete a run of 
the program: 

%a.out <retum> 

Input filename: 
filename <retum> 

Then die output of the temperatures of die different sides of die satellite and the nodes of 
die satellite will be output after the program is finished running. 

3.0 Thermal Control Example 

3.1 Example Background 

This example was first generated with preliminary information cm the ISAT-1 
satellite. The configuration of die satellite is a hexagonal cylinder with an effective diameter 
of 0317 meters and a height of 0.643 meters. For die preliminary analysis of this satellite, 
it was known that the satellite would be broken up into three main shelves: a battery shelf, 
an experiments shelf, and a communications eqmpment shelf. This configuration is shown 
in Figure 1, and thus leads to a six nodal analysis. 


,Top 


Level 3 

Communications 


Level 2 
Experiments 


Level 1 
Batteries 


.Sides 


Bottom 


Top 


Sides 


Bottom 


R3 


r/IXl 


R7 


Level 3 


Level 2 

M3M 

Level 1 


R* Thermal Resistivity 


Figure 1: ISAT-1 Preliminary Satellite Configuration 

With this configuration, information on each of the six nodes, as outlined in the 
previous section, was needed in order to run the program. For this particular example, the 
thermal control system of die satellite was being designed, so, none of die R-values 
between the nodes were known. Values were assumed for these R-values, and the output 
for these values will be shown. 


3.2 Input Hie 

The mbit for die ISAT-1 satellite was specified as being 95 minutes king with a 
shadow time of approximately 35 minutes. The height of the mint was 550 kilometers. 
These values are reflected in the input file shown on the next two pages. 
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Information on the different nodes was not completely known at the time of 
generating this example. However, it was known that the outer shell of the satellite was to 
be made of 15 kg of Aluminum. The sides of the satellite were to be covered with solar 
cells, and the top and bottom of the satellite could be covered with whatever is needed for 
die thermal control of the satellite. For this example, the bottom and top of the satellite 
were assumed to be covered with white paint (alpha = 0.25, epsilon = 0.85). It should 
also be noted, that the view factors used in this example were calculated from an average 
between a view factor for a flat plate above the earth, and a small sphere above the earth. [2] 
The input file reflects the values that were assumed for the other different properties of the 
satellite. The nodes for the input file were denoted so that node 4 was die battery shelf, 
node 5 was the experiments shelf, and node 6 was the communications shelf. This input 
file was named isatd, and is shown on the next two pages. 
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isat.d Input File 


95.0 20.0 55.0 

1430.0 0.3 230.0 

0.317 0.643 0.0 

273.0 
6 

10.0 
2.5 
2.5 

25.0 

25.0 

25.0 

900.0 

900.0 

900.0 

700.0 

890.0 

850.0 
0.0 
0.0 
0.0 

18.0 

18.0 

4.0 
0.72 
0.25 
0.25 
0.0 
0.0 
0.0 
0.71 
0.85 
0.85 
0.0 
0.0 
0.0 
0.727 
0.0 
0.727 
0.0 
0.0 
0.0 

1.0 
1.0E-1 
1.0E-1 
1.0E-1 
1.0E9 
1.0E9 

1.0 
1.0E9 
1.0E9 
1.0E9 



l.OE-1 

1.0 

l.OE-1 

1.0E9 

1.0E9 

1.0 

1.0E9 

1.0E9 

1.0 

l.OE-1 

1.0 



3.3 Output 


After the input file was prepared, isaLftn was compiled and ran. The procedure 
with output is shown below: 

$ f77 isaLftn 
isaLftn: 

$a .out 

ENTER INPUT FILENAME: 


SIDES 

THIGH TLOW TMEAN 

308.44610 295.68741 302.06675 

TOP 

THIGH TLOW TMEAN 

304.89307 297.99475 301.44391 

BOTTOM 

THIGH TLOW TMEAN 

306.18924 297.21547 301.70235 

NODE 4 

THIGH TLOW TMEAN 

306.51317 299.74645 303.12981 

NODE 5 

THIGH TLOW TMEAN 

305.81584 304.90375 305.35980 

NODE 6 

THIGH TLOW TMEAN 

304.79974 302.94559 303.87267 


This output shows the temperatures for each node for the R-values shown in the input file. 
Notice, that these temperatures are in degrees kelvin. This output shows that with these R- 
values, there could be some nodes which are at too high of a temperature. By changing the 
R-values in the input file and running the program, die corresponding change in node 
temperatures can be found. 

Another form of output that can be obtained from the p rogram by adding sim ple 
write statements to die program, is the temperature output at each time step far any n od c s , 
This information can be used to make temperature verses time plots for any node of the 
satellite. Figure 2 shows a graph of this type, for the sides (Node 1) of the satellite 
described by this example. 
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Figure 2: Temperature Response of Sides of Satellite 


It can be seen from this graph, that the t emp erature of the sides starts out at the initial 
temperature of 273 Kelvin, and then as the satellite continues going through orbits, die 
temperature readies an equilibrium fluctuation of temperature. The program is designed to 
take this equilibrium condition into account before reading final hi g h , low, and mean 
temperatures for each node of the satellite. 

4.0 Program Flow Chart 

The program described by this handbook was developed to study the thermal 
control needs of a satellite. A satellite in space, will interact with die space environment 
which will cause the temperature of die satellite skin to change. These changes is skin 
temperature will cause changes in die temperature of different com ponents erf die satellit e 
depending on die thermal contact between the components and die skin and the components 
with each other. The equation that describes the transient thermal response of a body in 
space is a first order nonlinear differential equation. This program was developed so that 
the temperature of die satellite would be evaluated at each second time step in its orbit Due 
to this fact a numerical differential equation solver was imp lemente d. The differential 
equation solver and temperature function will be discussed further in the next section of this 

mannal. 

Since a satellite is comprised of many components, the temperature of each 
component of the satellite has to be evaluated. In order to perform this, the satellite can be 
broken up into nodes, and the temperature of each node can be evaluated using the 
differential equation solver. The whole analysis results in a program that functions by 
using arrays to describe the different properties of each node of the satellite. The problem 
of finding the temperature at each node is handled by inc rementing die node numbers, and 
solving the temperature function for each node at each time step. Obviously, only the first 
three nodes (the sldn of the satellite) will be affected in different ways by thic external space 
environment, and will radiate heat to space. While, with the other nodes (internal 
\ components) the temperature function will remain the same. These different types of heat 

flux on a node were handled by calling different subroutines to evaluate all of the heat flux 
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cm a node i.e. for internal components, heat flux to and from the space environment was set 
to zero. 


A flow chart of the program is shown in Figure 3. This flow chart shows the 
general method that the program follows to solve for the temperature of different parts of a 
satellite at each time step in it’s orbit It should be noted that the program is set to simulate 
a certain number of orbits in order to reach an equilibrium fluctuation of temperature during 
an orbit The output temperatures of the satellite are not evaluated in the program until the 
equilibrium condition has been reached. 



Figure 3: isat.ftn Program Flow Chart 


5.0 Main Program 

5.1 Temperature Function 

The main point of this program is to solve a temperature function for each node of 
the satellite. This temperature function is a first order nonlinear differential equation that is 
solved for each node at each time step in the orbit of the satellite. The temperature function 
is shown below in equation 2. 
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dT (j) 


dt m (J) Cp(i) 


(Qiat(i) + Qs(D + Qa(J) + Qiif(i) + 


- £(^(J>g(TtD 4 -To 4 )) 

R(jf j 


i 3 


Eq. 2 


From equation 2, it can be seen that there are many different factors which affect the 
nodes of the satellite configuration. The first terms are the solar, albedo, and infrared 
radiation heat fluxes to the external skin of the satellite. The functions for these heat fluxes 
are described in more detail in die next section of this manual. For all internal nodes, these 
values are set to zero. Another term, is the internal heat generation of each node. These 
values are input to the p rogra m and are typically zero for die external nodes of the satellite. 
The next term deals with the heat flux due to thermal contact between nodes. This term 
takes the temperature difference between each node, divides this by the thermal resistance 
between the nodes, and determines the heat flux into or out of the node in question. The 
last term in this equation, is die satellite radiation to space. This radiation is governed by 
the Stcfan-Boltzmann fourth power law, and only is taken into account for the external 
nodes of the satellite. 

It is seen from equation 2, that the temperature function is written in terms of 
arrays. This allows the program to iterate to evaluate each node at each time step in the 
solving of the temperature of the satellite. The format of this equation also makes it easy to 
change the factors that are in the temperature function without changing the function itself, 
if it is desired to modify the program. 

5.2 Ordinary Differential Equation Solver 

The ordinary differential equation solver used in this program is the Runge-Kutta 
numerical differential equation solving scheme.[3] The form of the differential equation 
solver was modified to adapt to the arrays used in the program. The general form of the 
Runge-Kutta method is: 

7 * = F(x,y) 


I Wj Kj 

where: y = Temperature 
x = Time 

W = Gradient factors 
K = Time step 


In the program, the differential equation solver is set up to solve for the temperature 
of each node in terms of second increments in time. This time increment results in a 
program running time that is very long. Without sacrificing accuracy, a time step of every 
ten seconds was implemented in the program to decrease die running time. This time step 
is represented by the variable "H" in the listing of the program found in the Appendix of 
this manual It is recommended that this time step not be increased any more than thirty 



15 

seconds, because there is a trade off with die accuracy of the program due to die fact thflt 
die R-values input in the program cannot be as large eras small as with a small rime step. 
Reference 4 explains this "stability criteria" in more detaiL[4] 

6.0 Heat Transfer Algorithms 

The purpose of this section is to describe the method used to describe die different 
types of heat transfers that are in the temperature function. The first three parts to this 
section describe die radiation heat transfer between die space environment and the outside 
of the satellite. The general equation which governs these heat transfers, or heat fluxes, is 
shown in equation 3. 

Qrad = J • Aj * Alpha * View Factor • Area Factor Eq. 3 

where: J = Radiation intensity (W/m 2 ) 

At = Satellite area affected by the radiation (m 2 ) 

Alpha = Absorptivity of the satellite skin 
View Factor = 1 for Solar, and other values for albedo and 
infrared 

Area Factor = Factor which describes the amount of 
projected area of the satellite surface, 
depending on orientation 

From this equation, the intensity and alpha values, and view factors are input into the 
program, and the areas are calculated from die dimensions input into the program. The l«gt 
variable, the area factor, is what comprises most of the differences for each type of 
radiation cm die different sides of the satellite. These factors are different for each type of 
radiation and each section of die outside of the satellite. The factors were developed for the 
ISAT-1 satellite which has an orbit orientation described by Figure 4. From figure 4, it can 
be seen that the radiation in c i de n t on die different sides of the satellite . will vary depending 
on die time of the orbit and die orientation of the satellite at that time. 



Figure 4: Satellite Orientation During Orbit 
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The other heat transfer factors in the temperature function are fairly self 
explanatory, but, the method for handling die heat flux due to thermal contact between two 
nodes will also be discussed briefly in this section. For more information on the different 
types of radiation explained in these sections and in die temperature function of equation 3, 
see Reference S.[S] 

6.1 Solar 

The solar radiation will only be incident on the top of the satellite and the sides of 
the satellite as can be seen in Figure 4. For this reason, the area factor for the bottom of the 
satellite was set to zero. The top of the satellite will receive 100 percent solar radiation 
when it is perpendicular to the sun and receive 0 percent of the radiation in the shadow of 
the earth. In between these two, the incident radiation was assumed to decrease and 
increase in a linear fashion. This is depicted in Figure 5, which shows the top area factor 
for solar radiation for the period of one orbit This particular orbit is 95 minutes long with 
a 35 minute shadow time, but the program will handle calculating the area factor for other 
orbit parameters. The equations that govern all of the different area factors used in this 
program can be found in the program listing in the Appendix of this manual. The different 
area factor equations are found in their respective radiation subroutines in the listing. 



Figure 5: Solar Area Factor for Top of Satellite 

The area factor for the side of the satellite was handled in the same way as the top 
factor, but the sides of the satellite receive 100 percent radiation just before entrance into the 
shadow of die earth, and receive 0 percent radiation perpendicukur to die sun. The side 
solar area factor is depicted by Figure 6. 
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Figure 6: Solar Area Factor for Sides of Satellite 


62 Albedo 

The albedo radiation is solar radiation reflected off of the earth back into space. For 
this reason, the albedo radiation will only be incident upon the bottom of die satelli te and 
die sides of the satellite. The area factor for albedo for die top of the satellite was thus set 
to zero. The area factors for the other two sides was Handled a little differently than die 
solar factor. As the satellite comes close to entering the shadow of die earth, die albedo 
will start to decrease due to the fact that die shadow will not reflect to die satellite. Also, 
after die satellite has entered the shadow, their will still be a small amount of albedo 
affecting the sa t elli t e. For this reason, die albedo area factor was decreased starting ten 
m i n u t es before entering the shadow of die earth, and does not reach zero until five minutes 
after entry into the shadow. This is depicted by Figure 7. It should be mentioned, that in 
this program, the albedo area factor carried the same value for both die bottom and side s of 
the satellite. 
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Figure 7: Albedo Area Factor for Sides and Bottom of Satellite 

6.3 Infrared Radiation 

The infrared radiation on incident on a satellite remains constant due to die fact that 
die earth is always emitting this radiation throughout the entire orbit The infrared 
radiation, however, will only affect the bottom and sides of the satellite. Because a view 
factor is already used for determining die heat flux due to this radiation cm the satellite, and 
because die radiation is a constant throughout the orbit, die infrared radiation area factor 
was set to one for die bottom and sides of the satellite. This area factor was set to zero for 
die top of the satellite because it does not receive any of die radiation. 

6.4 Thermal Contact Flux 

It is important to understand the method used for determining the thermal contact 
heat flux between nodes for this program, because it depends on the R-values, which are 
die thermal control system. The general equation for finding die heat transfer between two 
nodes and different temperatures separated by a thermal resistor is seen in equation 4. 


Ql2 = 


Tl-T 2 

Rl2 


where: T = Temperature of the node (K) 

R = Thermal resistance between nodes (K/W) 
Q = Heat Transfer between nodes (W) 


Eq. 4 


In the program, the heat flux between nodes was evaluated for each node at each time step. 
This procedure can be seen in the temperature function subroutine of the program listing 
found in the Appendix of this manual. For more information on the method for performing 
a nodal analysis as explained by equation 4, see Reference 6.[6] 


6 -/ 
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8.0 Appendix: Program Listing 

PR OGRAM THERMALCCNTROL 

C THERMAL CONTROL DESIGN PROGRAM FOR THE I SAT-1 SATELLITE 
C 

C THERMAL CONTROL GROUP 
C 

C UNIVERSITY OF IONA 
C 

C GROUP MEMBERS: JOE CLAY, DAVE PRALL, MATT SCHNEIDER, 

C MIKE SKRBICH 

C 

C 

C SPRING 1993 
C 

IMPLICIT REAL* 8 (A-H,0-Z) 

DIMENSION SM(50) , CP (50), Q(50), T(50), ALPHA (50), EPS (50) 
DIMENSION A (50), TH(50), TL(50), TM(50), VF(50), R(50,50) 
DIMENSION TSUB(50) 

CCM40N/ ORPAR/ ORB, TORBCO, SHADEIN, SHADEOUT, TIMEORB 
COM4QN/OTHPAR/TCON, APROJ, SIGMA, SO, RINFRA, ALB, TC2 
COMON/VAR/NV 


C READ INPUT FILE 
C 

CHARACTER* 8, INNAME 

PRINT *, ’ENTER INPUT FILENAME:' 

READ ' (A) ' , INNAME 

OPEN (UNIT - 1, FILE - INNAME, STATUS - 'OLD') 

READ (1,2) TIMEORB, SHADEIN, SHADEOUT 

READ (1,2) SO, ALB, RINFRA 

READ (1,2) DIA, SLEN, TIME 

READ (1,3) T(l) 

DO 7 I- 2,6 
T(I) - T (1) 

7 CONTINUE 

i READ (1,900) NV 

900 FORMAT (13) 

DO 5 IO - 1,6 
DO 6 IP « 1,NV 

IF (IO.EQ.1) THEN 
READ (1,3) SM(IP) 

END IF 

IF (IO.EQ.2) THEN 
READ (1,3) CP (IP) 

ENDIF 

IF (IO.EQ.3) THEN 
READ (1,3) Q(IP) 

ENDIF 

IF (IO.EQ.4) THEN 
READ (1,3) ALPHA (IP) 

ENDIF 

IF (IO.EQ.5) THEN 
READ (1,3) EPS (IP) 

ENDIF 

IF (IO.EQ.6) THEN 
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READ (1,3) VF (IP) 

ENDIF 
CONTINUE 
CONTINUE 
K « 0 
L - 1 

IF (K.EQ.NV) THEN 
L - L + 1 
K - L - 1 
ENDIF 
K - K + 1 
READ (1,3) R(L,K) 

R(K,L) - R(L,K) 

IF (L.LE.NV-1) THEN 
GO TO 4 
ENDIF 

FORMAT (3F13.5) 

FORMAT (F13.5) 

CLOSE (1) 

CONSTANTS 

TCON CONVERTS TIME TO MINUTES IN THE TEMPERATURE FUNCTION 
SIGMA IS THE STEFAN BOLTZMAN CONSTANT 

ORB AND TORBCO ARE USED TO COUNT ORBITS IN THE PROGRAM 
DELTA TIME - 1 SEC 

TC2 - 60.0D0 

TIMEORB - TIMEORB * TC2 

SHADEIN - SHADEIN * TC2 

SHADEOUT « SHADEOUT * TC2 

ORBITS - 12.0D0 

SIGMA - 5.6703D-8 

ORB - 1.0D0 

TORBCO - TIMEORB 

TCON - 1.0D0 

MAIN 

AREAS CALCULATES DIFFERENT AREAS FROM INPUT DIMENSIONS 
CALL AREAS (DIA, SLEN, A, APROJ) 

ODE SOLVER: RUNGE-KUTTA METHOD 

H - 10 

NSTEPS - (TIMEORB* ORBITS) /H 
KPRINTT - 300 
DO 10 I - 1,NV 
TH(I) - 0.0D0 
TL(I) - 1.0D3 
CONTINUE 

DO 20 Y - 1, NSTEPS 

IF (KPRINTT. EQ. 300) THEN 
WRITE (6,101) T(4) 

KPRINTT - 10 

FT-SF. 

KPRINTT - KPRINTT + 10 
ENDIF 

CALL PRLONG (PRINTT, T) 
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IF (ORB.GE.12.ODO) THEN 
DO 30 ICOUN - 1,NV 

IF (T (ICOUN) .GT.TH (ICOUN) ) THEN 
TH (ICOUN) - T (ICOUN) 

END IF 

IF (T (ICOUN) .LT.TL (ICOUN)) THEN 
TL (ICOUN) - T (ICOUN) 

END IF 
30 CONTINUE 

END IF 

DO 40 I - 1,NV 

CALL VECADD (T, 0, TSUB) 

FKO - H*F (TIME, TSUB, I, R, Q, EPS, A, CP, SM, ALPHA, VF) 

CALL RES (TSUB) 

CALL VECADD (T,0.5D0*RK0, TSUB) 

RKL - H*F(TIME+0.5D0*H,TSUB,I,R,Q,EPS,A,CP,SM,ALPHA,VF) 

CALL RES (TSUB) 

CALL VECADD (T, 0 . 5D0*RK1, TSUB) 

RK2 - H*F (TIME+0. 5D0*H,TSUB,I,R, Q,EPS,A,CP,SM, ALPHA, VF) 

CALL RES (TSUB) 

CALL VECADD (T, RK2*1 . OdO, TSUB) 

RK3 - H*F (TIME+H, TSUB, I, R, Q, EPS/ A, CP, SM, ALPHA, VF) 

T(I) - T(I) + (RK0+2.0d0*RKL+2.Od0*RK2+RK3)/6.Od0 
40 CONTINUE 

TIME « TIME + H 
20 CONTINUE 

DO 50 IZ - 1,NV 

TM(IZ) - (TH(IZ)+TL(IZ))/2.0D0 
WRITE (6,*) 

IF (IZ.EQ.l) WRITE (*,*) • 

IF (IZ.EQ.2) WRITE (*, *) ' 

IF (IZ.EQ.3) WRITE (*,*) • 

IF (IZ.GT.3) WRITE (*,*) ' 

WRITE (6,*) (' THIGH TLOW 

WRITE (6,101) TH(IZ) ,TL(IZ) ,TM(IZ) 

50 CONTINUE 

100 FORMAT (3E13.5) 

101 FORMAT (3F13.5) 

END 

C 

C SUBROUTINES FOR SAT. HEAT TRANSFERS 
C 

C TEMPERATURE FUNCTION 

FUNCTION F (TIME, T, I, R, Q, EPS, A, CP, SM, ALPHA, VF) 

IMPLICIT REAL*8 (A-H,0-Z) 

C0M40N/ORPAR/' ORB, TORBCO, SHADEIN, SHADEOUT, TIMEORB 
COMMON/ OTHPAR/TCON, APROJ, SIGMA, SO, RINFRA, ALB, TC2 
COM40N/VAR/NV 

DIMENSION R(50, 50) , Q (50) , EPS (50) , A (50) , SM(50) , CP (50) , T (50) , VF (50) 
(ALPHA (50) 

C INTEGER I 

CALL ORBCOUNT (TIME, ORB, TORBCO, TIMEORB, TIME2) 

CALL SOLAR (I, TIMEORB, SO, TIME2, SHADEIN, SHADEOUT, A, ALPHA, 

&APROJ, QSQLAR) 

CALL ALBEDO (I, ALB, SO, TIME2, SHADEIN, SHADEOUT, ALPHA, VF, A, TC2, 


SIDES' 

TOP' 

BOTTOM' 

NODE ',IZ 

TMEAN') 
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CQALBEDO) 

CALL INFRARED (I, EPS, A, VF,RINFRA,QINFRA) 

TDR - 0.0D0 
DO 60 J - 1,NV 

TDR - TDR + ( (T (J) —T (I) ) /R(I, J) ) 

60 CONTINUE 
119 FORMAT (F13. 5) 

F “ ( (Q(I) + QSQLAR+ QALBEDO+ QINFRA+ TDR - SIGMA * EPS (I) * 
*A(I) * T(I)**4.0D0) * TCON) / (SM(I)*CP(I) ) 

12 FORMAT (6F13. 5) 

RETURN 

END 

C LAUNCH 
C ORBCOUNT 

SUBROUTINE ORBCOUNT (TIME, ORB, TORBCO, TIMEORB, TIME2) 

IMPLICIT REAL* 8 (A-H, O-Z) 

IF (TORBCO.LT. TIME) THEN 
ORB - ORB + 1.0D0 
TORBCO - ORB * TIMEORB 
ENDIF 

TIME2 - TIME - TIMEORB* (ORB-1 . 0D0) 

RETURN 

END 

C AREAS 

SUBROUTINE AREAS (DIA, SLEN, A, APROJ) 

IMPLICIT REAL* 8 (A-H, O-Z) 

DIMENSION A (50) 

CCM40N/VAR/NV 

PI - 3.14159D0 

A(l) - PI * DIA * SLEN 

A (2) - PI * ( (DIA**2 . 0D0) /4 . 0D0) 

A(3) -A(2) 

DO 301 IAO-4,NV 
A(IAO) - 0.0D0 
301 CONTINUE 

APROJ - DIA * SLEN 

RETURN 

END 

C SOLAR RADIATION 

C ONLY EFFECTS THE SHELL TOP AND SIDES I.E. I - 1,2 

SUBROUTINE SOLAR (I , TIMEORB, SO, TIME2, SHADEIN, SHADEOUT, 
tA, ALPHA, APROJ, QSOLAR) 

IMPLICIT REAL*8 (A-H, 0-2) 

DIMENSION A (50), ALPHA (50) 


IF (I.GE.3) THEN 
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QSOLAR - 0.0D0 
END IF 

IF (I.EQ.l) THEN 

IF (TIME2 . LE . SHADEIN) THEN 
FAC - TIME2 /SHADEIN 
END IF 

IF (TIME2 . GE . SHADEOUT) THEN 

FAC - - (TIME2/ (TIMEORB-SHADEOUT) ) +1 . ODO+ (SHADEOUT/ 

« (TIMEORB-SHADEOOT) ) 

END IF 

IF ( (TIME2.GT. SHADEIN) .AND. (TIME2.LT. SHADEOUT) ) THEN 
FAC - 0.0D0 
END IF 

QSOLAR - ALPHA(l) * SO * APROJ * FAC 
END IF 

IF (I.EQ.2) THEN 

IF (TIME2.LE. SHADEIN) THEN 

FAC - - (TIME2/SHADEIN) + 1.0D0 
END IF 

IF (TIME2 . GE . SHADEOUT) THEN 

FAC - (TIME2/ (TIMEORB-SHADEOUT) ) +1 . 0D0- (TIMEORB/ 
i (TIMEORB-SHADEOUT) ) 

END IF 

IF ( (TIME2 . GT . SHADEIN) . AND . (TIME2 . LT. SHADEOUT) ) THEN 
FAC - 0.0D0 
END IF 

QSOLAR - ALPHA (2) * SO * A (2) * FAC 
END IF 

RETURN 

END 

C INFRARED RADIATION 

C ONLY EFFECTS THE SIDES AND BOTTOM I.E. I - 1,3 
SUBROUTINE INFRARED (I, EPS, A, VF, RINFRA, QINFRA) 

IMPLICIT REAL*8 (A-H, O-Z ) 

DIMENSION A(50) ,VF(50) , EPS (50) 

QINFRA - EPS (I) * A (I) * VF(I) * RINFRA 

RETURN 

END 

C ALBEDO RADIATION 

C ONLY EFFECTS THE SIDES AND BOTTOM I.E. I - 1,3 

SUBROUTINE ALBEDO ( I , ALB, SO, TIME2 , SHADEIN, SHADEOUT, ALPHA, VF 
i ,A,TC2,QALBEDO) 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION A(50) ,VF(50) , ALPHA (50) 

S1IN - SHADEIN - 10.0D0 * TC2 
S2IN * SHADEIN + 5.0D0 * TC2 
S10UT - SHADEOUT - 5.0D0 * TC2 
S20UT - SHADEOUT + 10.0D0 * TC2 


FAC2 - 1.0D0 
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IF ((TIME2.GT.S1IN) .AND. (TIME2.LT. S2IN) ) THEN 
FAC2 - — (TIME2/ (15.0D0*TC2) ) + l.ODO + 

£ ( (S1IN) / (15. ODO*TC2) ) 

END IF 

IF ((TIME2.GT.S10UT) .AND. (TIME2.LT. S20UT) ) THEN 
FAC2 - (TIME2/ (15. 0D0*TC2) ) •+ l.ODO - 
£ ( (S20UT) / (15 . 0D0*TC2) ) 

END IF 

IF ( (TIME2.GE.S2IN) .AND. (TIME2.LE.S10UT) ) THEN 
FAC2 - 0.0D0 
END IF 

QALBEDO - SO * ALB * ALPHA (I) * A(I) * VF(I) * FAC2 

RETURN 

END 

C 

SUBROUTINE VECADD (T, SCAL, TSUB) 

IMPLICIT REAL*8 (A-H, O-Z) 

DIMENSION T (50), TSUB (50) 

CCM40N/VAR/NV 

DO 1 KP - 1,NV 
1 TSUB(KP) - T (KP) + SCAL 

RETURN 

END 

C 


SUBROUTINE RES (TSUB) 
IMPLICIT REAL* 8 (A-H, O-Z) 
DIMENSION TSUB (50) 
COM40N/VAR/NV 

DO 133 KL - 1,NV 
133 TSUB(KL) - 0.0D0 

RETURN 

END 



